
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 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\