search for: listof

Dear R developers A question about the class 'listof', defined in package 'stats'. Other than its definition and use in the code for 'anova', we can't see that the class 'listof' is used for anything else (in recommended packages, or elsewhere). In the spatstat package we have been using a 'listof' to repre...

...+Error(Run),data=assayvar) > str(summary(assayaov)) List of 2 $ Error: Run :List of 1 ..$ :Classes ?anova? and 'data.frame': 1 obs. of 3 variables: .. ..$ Df : num 3 .. ..$ Sum Sq : num 46.5 .. ..$ Mean Sq: num 15.5 ..- attr(*, "class")= chr [1:2] "summary.aov" "listof" $ Error: Within:List of 1 ..$ :Classes ?anova? and 'data.frame': 1 obs. of 5 variables: .. ..$ Df : num 8 .. ..$ Sum Sq : num 36.4 .. ..$ Mean Sq: num 4.55 .. ..$ F value: num NA .. ..$ Pr(>F) : num NA ..- attr(*, "class")= chr [1:2] "summary.aov" "listof&q...

...d(year,as.matrix(x[,1])) >lm.sol<-lm(x~year) >xtrend<-coef(lm.sol)[2,] # get the matrix of regression coefficient >t.test<- ? # also want to get a similar matrix of t-test value for the regression coefficient >class(summary(lm.sol)) > "listof" >class(summary(lm.sol)) >"listof" # the t-test values are in the obj "summary(lm.sol)", but how to get them out as a matrix similar as the above"xtrend"? Anyone...

...:List of 1 ..$ :Classes ‘anova’ and 'data.frame': 1 obs. of 5 variables: .. ..$ Df : num 4 .. ..$ Sum Sq : num 12.4 .. ..$ Mean Sq: num 3.1 .. ..$ F value: num NA .. ..$ Pr(>F) : num NA ..- attr(*, "class")= chr [1:2] "summary.aov" "listof" $ Error: subject:myfactor:List of 1 ..$ :Classes ‘anova’ and 'data.frame': 2 obs. of 5 variables: .. ..$ Df : num [1:2] 2 8 .. ..$ Sum Sq : num [1:2] 14.9 4.4 .. ..$ Mean Sq: num [1:2] 7.47 0.55 .. ..$ F value: num [1:2] 13.6 NA .. ..$ Pr(>F) : num [1:2] 0....

Dear r-users, Suppose I have three datasets: Dataset-1: Date x y Jan-1,2005 120 230 Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21 Jan-7,2005 11 299 Mar-5,2005 200 311 Dataset-2: Date x y Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21

...r: PPNR :List of 1 ..$ :Classes anova and `data.frame': 1 obs. of 5 variables: .. ..$ Df : num 70 .. ..$ Sum Sq : num 70.7 .. ..$ Mean Sq: num 1.01 .. ..$ F value: num NA .. ..$ Pr(>F) : num NA ..- attr(*, "class")= chr [1:2] "summary.aov" "listof" $ Error: PPNR:CONSENS:List of 1 ..$ :Classes anova and `data.frame': 2 obs. of 5 variables: .. ..$ Df : num [1:2] 1 70 .. ..$ Sum Sq : num [1:2] 0.254 70.746 .. ..$ Mean Sq: num [1:2] 0.254 1.011 .. ..$ F value: num [1:2] 0.251 NA .. ..$ Pr(>F) : num [1:2] 0.618...

...frame'': 2 obs. of 5 variables: ..$ Df : num [1:2] 1 18 ..$ Sum Sq : num [1:2] 0.000155 0.00058 ..$ Mean Sq: num [1:2] 1.55e-04 3.22e-05 ..$ F value: num [1:2] 4.82 NA ..$ Pr(>F) : num [1:2] 0.0415 NA - attr(*, "class")= chr [1:2] "summary.aov" "listof" - attr(*, "na.action")=Class ''omit'' Named int 5 .. ..- attr(*, "names")= chr "5" agradezco su atención y ayuda [[alternative HTML version deleted]]

...ame': 2 obs. of 5 variables: ..$ Df : num [1:2] 3 14 ..$ Sum Sq : num [1:2] 2.79e+10 3.08e+10 ..$ Mean Sq: num [1:2] 9.29e+09 2.20e+09 ..$ F value: num [1:2] 4.22 NA ..$ Pr(>F) : num [1:2] 0.0253 NA - attr(*, "class")= chr [1:2] "summary.aov" "listof" > attr(summary(pb), "Pr(>F)") NULL

...4 4 8 27 ..$ Sum Sq : num [1:5] 1232 4889 6472 4181 310 ..$ Mean Sq: num [1:5] 616.1 1222.4 1618 522.6 11.5 ..$ F value: num [1:5] 53.6 106.4 140.8 45.5 NA ..$ Pr(>F) : num [1:5] 3.96e-10 4.07e-16 1.14e-17 1.09e-13 NA - attr(*, "class")= chr [1:2] "summary.aov" "listof" I try to use results$Df or similar, but nothing seem to work Thanks in advance -- View this message in context: http://r.789695.n4.nabble.com/Problems-to-extract-data-from-anova-table-tp3723125p3723125.html Sent from the R help mailing list archive at Nabble.com.

...ist",representation(genDist = "numeric",roundDec = "integer")) Most things work fine, but when subsetting a haploList object an object of class list is returned. I realize that I need to write a function for subsetting this new object, and tried to find the code for '[.listof' or something similar could not find it, probably due to a suboptimal understanding of how it is organized. My question is, how could I define a extraction function for my new class that uses all the existing functionality of the Extract function for list? Thanks in advance, Albart Coster

...I suppose represent the effects on obs1 and obs2 combined. The example for summary.manova shows that > summary.aov(fit) would print the effects separately for each observable. However, I get the message Error in 1:object$rank : NA/NaN argument and > class(fit) [1] "aovlist" "listof" I suppose that this means that I shouldn't want to do this sort of test, and rather look at each observable separately. The question is: why? Is there a simple explanation comprehensible for the non-statistician? RenE Bertin. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-....

...ame': 2 obs. of 5 variables: ..$ Df : num [1:2] 1 269 ..$ Sum Sq : num [1:2] 60.2 1792.0 ..$ Mean Sq: num [1:2] 60.24 6.66 ..$ F value: num [1:2] 9.04 NA ..$ Pr(>F) : num [1:2] 0.00289 NA - attr(*, "class")= chr [1:2] "summary.aov" "listof" but if I try to extract it: test.summary[[2]]$'Pr(>F)' NULL I assume that the incantation has changed. Using R-2.3.1 - Tried this on both Windows XP and FC4 Linux Jim

...#39; and 'data.frame': 2 obs. of 5 variables: ..$ Df : num [1:2] 2 147 ..$ Sum Sq : num [1:2] 63.2 39 ..$ Mean Sq: num [1:2] 31.606 0.265 ..$ F value: num [1:2] 119 NA ..$ Pr(>F) : num [1:2] 1.67e-31 NA - attr(*, "class")= chr [1:2] "summary.aov" "listof" > tmp$1 Error: unexpected numeric constant in "tmp$1" > tmp$"1" NULL > tmp$"" NULL Thank you Liviu -- Do you know how to read? http://www.alienetworks.com/srtest.cfm Do you know how to write? http://garbl.home.comcast.net/~garbl/stylemanual/e.htm#e...

Dear all, Trying to extract a few rows for each element of a list of data.frames, I'm puzzled by the following behaviour, > d <- lapply(1:4, function(i) data.frame(x=rnorm(5), y=rnorm(5))) > str(d) > > lapply(d, "[", i= c(1)) # fine, this extracts the first columns > lapply(d, "[", j= c(1, 3)) # doesn't do nothing ?! > > library(plyr)

Dear all, Trying to extract a few rows for each element of a list of data.frames, I'm puzzled by the following behaviour, > d <- lapply(1:4, function(i) data.frame(x=rnorm(5), y=rnorm(5))) > str(d) > > lapply(d, "[", i= c(1)) # fine, this extracts the first columns > lapply(d, "[", j= c(1, 3)) # doesn't do nothing ?! > > library(plyr)

...o I wonder if you can tell me what went wrong with my function. > TukeyHSD(result.subject, c('Factor1','Factor2'), ordered=TRUE) Error in UseMethod("TukeyHSD") : no applicable method for 'TukeyHSD' applied to an object of class "c('aovlist', 'listof')" Also, I I attempted *pairwise.t.test* by averaging each participant's reading times in each condition (see the table below) with *p* values adjusted with Bonferroni method (*e.g. *pairwise.t.test(y,group, "bonferroni", pair=TRUE)), . It gave me significant results. But I...

...n-subjects), a number of subjects (Subj), and the dataframe is called WMU3C. The model is > RT.aov <- aov(RT~Cond + Error(Subj/Cond), WMU3C) I understand that TukeyHSD only works with an aov object, but that RT.aov is an aovlist object. > class(RT.aov) [1] "aovlist" "listof" I've tried to work around it using the "maiz" example in the MMC documentation of the HH package (a solution previously recommended), but I couldn't get it to work: My best shot was to calculate another aov avoiding the error term (I don't see how this could be a fea...

...#39;t work with a nested design!! (it does work if I do the exact same thing without the nested term. I get the following: anova(RSslopes,RSancova) #Error in UseMethod("anova") : # no applicable method for 'anova' applied to an object of class "c('aovlist', 'listof')" This has something to do with the nested design putting out an aov list, I think. Now I really lose it....because I know how to perform a Tukey HSD test to compare each combination of treatments if I had a regular non-nested ANCOVA, but when I try it with these tests it can't run....

...811 0.4811 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Residuals 30 15.78 0.5259 Althougth the Anova Model object it isn't a object of class "aov" but a object of class"aovlist": #> attr(AnovaModel.2,"class") [1] "aovlist" "listof" This class isn't suitable for the glht function neither TukeyHSD function. I can extract a object of class "aov" through: #> ls.str(pat="AnovaModel.2") AnovaModel.2 : List of 3 $ (Intercept):List of 9 $ patient :List of 9 $ Within :List of 6 #> An...

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