similar to: extraction of p-value

Dear R community, I'm currently using R to analyze functional Magnetic Resonance Imaging data. Each analysis involves running ~120,000 repeated-measures ANOVAs. I would like to know if there is any automatic way to access the F- and P-value data that are associated with each of these 120,000 ANOVAs. For example, if the summary output (for the 1st ANOVA of 120,000) shows the following value

Hello dear R members. I have been learning the Anova syntax in order to perform an SS type III Anova with repeated measures designs (thank you Prof. John Fox!) And another question came up: where/what are the (between/within) residuals for my model? ############ Play code: phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),

Hi, I've done a repeated ANOVA looking as follows: model<-aov(y~a*b*c+Error(d)) Now I want to get rid of some or all of the interactions. But neither AIC nor anova work. What is the appropriate function to test my model?

Dear Madame, Dear Sir, I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from a 'summary' of 'aov'. The following example shows my steps. ## Initialize rm(list = ls()) # remove (almost) everything in the working environment utils::data(npk, package="MASS") # get data model <- yield ~ block + N*P*K ## Using lm npk.lm <-

Hi all, Perhaps this in an inappropriate post, but I've found a bug in proj I'd like to track down a bit further before making a formal bug report. The example below shows the problem, if you change the rownames proj fails. The problem seems to be that there is a mismatch in the rownames in the qr objects constructed by aov and the rownames that proj is expecting them to have.

Dear listmembers, I would like to define a class with a slot that takes either an object of class aov or NULL. I have been reading "S4 Classes in 15 pages more or less" and "Lecture: S4 classes and methods" #First I tried with list and NULL setClass(listOrNULL") setIs("list", "listOrNULL") setIs("NULL", "listOrNULL") #doesn't

Hi, This might be a minor bug, or it could be that I'm just mis-using the tools (again). If you modify the example for model.tables to introduce an NA, N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0) P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0) K <- c(1,0,0,1,0,1,1,0,0,1,0,1,0,1,1,0,0,0,1,1,1,0,1,0) yield <-

A non-text attachment was scrubbed... Name: not available Type: text Size: 3126 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19991130/5cb00c0f/attachment.pl

R-help What are the strengths and weakness of 'aov' in 'car' package? My model looks something like this : library(car) aov(lm(fish.length~zone*area,data=my.data)) Thank you Luis Ridao Cruz Fiskiranns??knarstovan N??at??n 1 P.O. Box 3051 FR-110 T??rshavn Faroe Islands Phone: +298 353900 Phone(direct): +298 353912 Mobile: +298 580800 Fax:

Hello!!! I'm trying to do multiple ANOVA tests with R (testing the affect off different factors on the same response). As a result I get many ANOVA tables, and I want to extract a list of the "Pr(>F)" from all the tables. Maybe someone have an idea how to do this? Thanks Imri -- View this message in context:

Doea anyone know how to compute components of variance analysis? For example, I have the score of pupils on a test. Each pupil belongs to a school and within each school I may have several classes? How can I estimate the variance of the pupils, schools, classes and the errror variance? Any examples? Any help will be appreciated. Jason Dr. Iasonas Lamprianou Department of Education The

I have DHCP and TFTP on a Solaris machine, pxelinux.0 loads, but it can't load the configuration file. The tftp server on Solaris tries to send the file, but data transfers do not receive ACKs. Previous posts to this list have indicated that the suspect was the PXE BIOS (ref "PXELinux can't load the config file?" dated March 20, 2003). Does that still appear to be the case

Dear statistitians / R users, I was told to analyse the effects of the four factors using lm or update for aov. The following is the result from lm. As you can see that in the last few lines enclosed coefficients. 1)Does this mean I can write a formula, like yield =0.26 + 0.03*H + 1.48*T + 0.04*L + 0.004*C 2) in the two levels design, is there any difference between I use (-1, 1) to

Hello, I noticed that model.tables fails when applied to an aov() fit if called inside a function. The problem seems to occur when as.formula is used inside a function on a string containing "<formula> + Error( x / y )" The reason I tried to use as.formula is to generate dynamic calls to aov(). Here is a minimal example illustrating the problem: ## Example test <-

Is there a way to prevent the re-ordering of factors by aov? I do have a three-way interaction that I do want to fit before a two-way interaction (different factors, so they are not nested), but R moves the two-way interaction to the front. I know it generally makes sense to fit the two-way interactions first, but in this case I think I know what I'm doing. + hv

I'm baffled. When I run make check after installing from source, I get a Error 2. From my understanding of how these things work, it would appear to be coming from this (as at the end of base-Ex.Rout.fail: > has.VR <- require(MASS, quietly = TRUE) Attaching package 'MASS': The following object(s) are masked from package:base : confint confint.lm nclass.FD nclass.scott

Is there an easy way of converting an aov.summary into a matrix in which the rows are the factor names and the columns are Df, Sum Sq, Mean Sq, F value and Pr. For example, convert Df Sum Sq Mean Sq F value Pr(>F) block 5 343.29 68.66 4.4467 0.015939 * N 1 189.28 189.28 12.2587 0.004372 ** P 1 8.40 8.40 0.5441 0.474904 K

hello all, Quick question, how do I get the p value out of the anova? Thanks, Paul > pb<-aov(as.numeric(diff[5,16:33]) ~ grF) > summary(pb) Df Sum Sq Mean Sq F value Pr(>F) grF 3 2.7860e+10 9.2867e+09 4.2236 0.02534 * Residuals 14 3.0783e+10 2.1988e+09 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 '

Hi all I have the next data frame year STOD SLAGNR TAL TALT TALVEKT 1 2002 2120006 57 1 NA 1 2 1997 97030032 57 NA NA NA 3 1997 97030071 27 1 NA NA 4 1997 97030005 57 1 NA NA 5 1997 97020127 60 NA 1 NA 6 2001 1160025 27 1 NA 1 7 1998 98020069 60 1 NA NA 8

R-help I have a martix with missing values( in which I want the sample size by column) When I : apply(matrix,2,length) I get the length of the vector regardless of missing values. I can't pass an argument to length in apply. Alternatively I could ifelse ( is.na ( matrix [, "columns in matrix " ] ) , 0 , 1) Is there any easier way? Thank you