[R-br] Seria um caso de R inferno?

Jose Claudio Faria joseclaudio.faria em gmail.com
Segunda Novembro 3 08:03:01 BRST 2014


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
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141103/c5cb3ea7/attachment.html>


Mais detalhes sobre a lista de discussão R-br