>>>>> "Adam" == Adam Witney <[EMAIL PROTECTED]>
>>>>> on Fri, 25 May 2007 14:48:18 +0100 writes:
>>> ##-- 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)
Adam> Here is the results:
>> xB <- c(2000,1e6,1e50,Inf)
>> for(df in c(0.1, 1, 10))
Adam> + for(ncp in c(0, 1, 10, 100))
Adam> + print(pchisq(xB, df=df, ncp=ncp), digits == 15)
Adam> Error in print.default(pchisq(xB, df = df, ncp = ncp), digits == 15) :
Adam> object "digits" not found
well, that's a typo - I think - you should have been able to fix
(I said "something like" ...).
Just do replace the '==' by '='
Martin
Adam> Thanks again...
Adam> 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.
