>> ##-- non central Chi^2 :
>> xB <- c(2000,1e6,1e50,Inf)
>> for(df in c(0.1, 1, 10))
> + for(ncp in c(0, 1, 10, 100)) stopifnot(pchisq(xB, df=df, ncp=ncp) ==1)
> Error: pchisq(xB, df = df, ncp = ncp) == 1 is not all TRUE
> Execution halted
>
> Ok, thanks;
> so, if we want to learn more, we need
> the output of something like
>
> xB <- c(2000,1e6,1e50,Inf)
> for(df in c(0.1, 1, 10))
> for(ncp in c(0, 1, 10, 100))
> print(pchisq(xB, df=df, ncp=ncp), digits == 15)
Here is the results:
> xB <- c(2000,1e6,1e50,Inf)
> for(df in c(0.1, 1, 10))
+ for(ncp in c(0, 1, 10, 100))
+ print(pchisq(xB, df=df, ncp=ncp), digits == 15)
Error in print.default(pchisq(xB, df = df, ncp = ncp), digits == 15) :
object "digits" not found
Thanks again...
adam
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
