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

Jose Claudio Faria joseclaudio.faria em gmail.com
Segunda Novembro 3 19:39:39 BRST 2014


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 . <walmeszeviani em gmail.com> 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
> 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/972d4191/attachment.html>


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