[R] Increase Decimal Accuracy?

2006-03-02 Thread Ann Hess
Is there any way to increase the decimal accuracy for probability distributions. For example (in R): 1-pchisq(90,5) [1] 0 But in Maple, I find that the value is 0.67193x10-17. I need to compare some really small p-values...is there a way to increase the decimal place accuracy beyond 1x10-16

[R] contrasts for lm

2005-12-07 Thread Ann Hess
I would like estimate a number of contrasts from a one-way ANOVA model. I see that the lm command has a contrasts option, but I can't figure out how to use it! Any help that can be offered would be greatly apreciated. Here is my model statement: Model-lm(log2PM~P+T+P*T) where P has 16

[R] .Random.seed

2003-10-16 Thread Ann Hess
I am writing a function for the purposes of a simulation. Due to memory problems, the function sometimes crashes. In order to get around this problem, I would like to include to be able to save the last seed, so I can pick up with the next run of the simulation after a crash. I am having