<div dir="ltr">Olá,<div><br></div><div>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?</div><div><br></div><div><div># Função de densidade: chisq</div><div>fchisq <- function(x, n)</div><div>  (1/(2^(n/2)*gamma(n/2)))*x^(n/2-1)*exp(-x/2)</div><div><br></div><div>par(mfrow=c(1, 2))</div><div>plot(fchisq(x=1:100, </div><div>            n=10))</div><div>curve(dchisq(x, </div><div>             10),</div><div>      col='red',</div><div>      lw=2,</div><div>      add=TRUE,</div><div>      n=1e3)</div><div>## OK: a função fchisq está correta!</div><div><br></div><div>plot(fchisq(x=seq(0, </div><div>                  50, </div><div>                  len=100), </div><div>            n=10))</div><div>curve(dchisq(x, </div><div>             10),</div><div>      col='red',</div><div>      lw=2,</div><div>      add=TRUE,</div><div>      n=1e3)</div><div>## ?</div><div><br></div><div>Gastei um bom tempo para perceber que o uso de "seq" não é correto nesses casos.</div><div><br></div><div>Ab,</div><div><div class="gmail_signature">///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>Jose Claudio Faria<br>Estatistica<br>UESC/DCET/Brasil<br>joseclaudio.faria at <a href="http://gmail.com" target="_blank">gmail.com</a><br>Telefones:<br>55(73)3680.5545 - UESC<br>55(73)9966.9100 - VIVO<br>55(73)9100.7351 - TIM<br>55(73)8817.6159 - OI<br>///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\</div></div>
</div></div>