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

Téo Calvo teocalvo2 em gmail.com
Segunda Novembro 3 11:31:32 BRST 2014


Esqueci de mencionar:

x=1:100
y=seq( from=1, to=100, by=1)
y==x

Talvez esclareça um pouco mais.

Att Téo Calvo.

On 03-11-2014 08:03, Jose Claudio Faria wrote:
> 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 <http://gmail.com>
> Telefones:
> 55(73)3680.5545 - UESC
> 55(73)9966.9100 - VIVO
> 55(73)9100.7351 - TIM
> 55(73)8817.6159 - OI
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141103/94043aca/attachment.html>


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