[R] glmmPQL in 2.3.1

2006-09-25 Thread Justin Rhodes
Dear R-help, I recently tried implementing glmmPQL in 2.3.1, and I discovered a few differences as compared to 2.2.1. I am fitting a regression with fixed and random effects with Gamma error structure. First, 2.3.1 gives different estimates than 2.2.1, and 2.3.1, takes more iterations to con

[R] power analysis for 2-way anova interaction

2006-04-13 Thread Justin Rhodes
Dear R-help, Does anyone have a function that I could use to determine power for testing significance of an interaction term in a 2-way ANOVA. The model I am considering has 2 levels in each factor. I would like to be able to input sample size (equal in all groups), MSE, means for the 4 gro

[R] power analysis for 2-way anova interaction

2006-04-13 Thread Justin Rhodes
Dear R-help list, Does anyone have a function that I could use to determine power for testing significance of an interaction term in a 2-way ANOVA. The model I am considering has 2 levels in each factor. I would like to be able to input sample size (equal in all groups), MSE, means for the 4

Re: [R] R-square n p-value

2005-09-04 Thread Justin Rhodes
ce now recently and it is incredibly helpful and fast. Thanks again. It is much appreciated. cc Dirk Best, Justin At 03:43 AM 9/4/2005, Peter Dalgaard wrote: >Dirk Eddelbuettel <[EMAIL PROTECTED]> writes: > > > On 3 September 2005 at 17:59, Justin Rhodes wrote: > >

[R] R-square n p-value

2005-09-03 Thread Justin Rhodes
Dear R-help, Can someone please help me discover what function or code will give me a p-value from the input: 1) R-square statistic from a simple linear regression, and 2) sample size, n This would be greatly appreciated. I need this because I am using a database that gives me R-square and sa

[R] Indexing within scan

2005-07-19 Thread Justin Rhodes
something like: results<-numeric(24) for (i in 1:24) { p<-scan("ai.txt") # where the i is an index for each of the 24 files results[i]<-mean(p) } Thanks for any help, Justin Rhodes Behavioral Neuroscience Oregon Health & Science University VA Medical Center (R &

[R] Question about setting tick mark limits in Plot

2004-12-29 Thread Justin Rhodes
p. I really appreciate it. I have wasted too much time trying to figure this out on my own. par(xaxp=c(0,240,8)) par(mar=c(13,1,13,1)) plot(p,q, pch="|", yaxt="n", main=subject,xlab="Minutes", ylab="") Take care, Justin Rhodes Behavioral Neuroscience