Caros Téo e Walmes, Desculpe a demora na resposta: estava no campo sem web. Muito obrigado: ambas as soluções são satisfatórias!
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 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ On Mon, Nov 3, 2014 at 1:32 PM, walmes . <[email protected]> wrote: > Agora que pude testar o código, acho que encontrei a causa de falha. > > ## 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) > plot(fchisq(x=seq(0, 50, len=100), n=10)) > curve(dchisq(x, df=10), col='red', lw=2, add=TRUE, n=1e3) > > *## Estava sendo usado o INDEX e não o valor de x!!!* > > x <- seq(0, 50, length.out=101) > plot(fchisq(x=x, n=10)~x) > curve(dchisq(x, df=10), col='red', lw=2, add=TRUE, n=1e3) > > À disposição. > Walmes. > > > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. >
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
