Olá,
Estava preparando um material didático para aulas de
introdução à estatística, na parte das variáveis aleatórias
contínuas e me deparei com um problema que ainda não encontrei
resposta. Seria um caso de R inferno?
# Função de densidade: chisq
fchisq <- function(x, n)
(1/(2^(n/2)*gamma(n/2)))*x^(n/2-1)*exp(-x/2)
par(mfrow=c(1, 2))
plot(fchisq(x=1:100,
n=10))
curve(dchisq(x,
10),
col='red',
lw=2,
add=TRUE,
n=1e3)
## OK: a função fchisq está correta!
plot(fchisq(x=seq(0,
50,
len=100),
n=10))
curve(dchisq(x,
10),
col='red',
lw=2,
add=TRUE,
n=1e3)
## ?
Gastei um bom tempo para perceber que o uso de "seq" não
é correto nesses casos.
Ab,
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
Jose Claudio Faria
Estatistica
UESC/DCET/Brasil
joseclaudio.faria at
gmail.com
Telefones:
55(73)3680.5545 - UESC
55(73)9966.9100 - VIVO
55(73)9100.7351 - TIM
55(73)8817.6159 - OI
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\