search for: 3df

...kages(all.available =3D TRUE))) + if(nchar(pkg)) library(pkg, character.only=3DTRUE)}) > library(GOAT) > initGOAT() [1] "Showing the Tk GUI..." [1] "Analyzing enrichment..." Done. > > write.table(1h.txt.proc, file=3D"1h.proc.xls", sep=3D"\t",quote=3DF, col.= names=3DF, row.names=3DF) Error: syntax error I checked the file names, but still I got that message. Additionally here the instructions I followed step by step: turn on R 1.9.1. load package tcltk type in: library(GOAT) initGOAT() #a new window will open select your inputfile and your...

Hello. I have come across a curious result that I cannot explain. Hopefully, someone can explain this. I am doing a 1-way ANOVA with 6 groups (example: summary(aov(y~A)) with A having 6 levels). I get an F of 0.899 with 5 and 15 df (p=0.51). I then do the same analysis but using data only corresponding to groups 5 and 6. This is, of course, equivalent to a t-test. I now get an F of 142.3

Dear Rologists, In order to combine plots I need to get access to the some "par"s specific to my plot prior to replot it with modified parameters. I have not found any option like "plot=F" associated with truehist and would like to know whether someone can point out how to overcome this problem. Thanks, Joh

...flush.console() # (5) cumulate <- 0; c <- 0; qra <-0 for (i in 1:runs) { a <- rnorm(700*700); dim(a) <- c(700,700) b <- 1:700 timing <- system.time({ qra <- qr(a, tol =3D 1e-10); c <- qr.coef(qra, b) #Rem: a little faster than c <- lsfit(a, b, inter=3DF)$coefficients })[3] cumulate <- cumulate + timing } timing <- cumulate/runs times[5, 1] <- timing cat(c("Linear regression over a 700x700 matrix (c =3D a \\ b') (sec): ", = timing, "\n")) remove("a", "b", "c", "qra") flush...

I do not understand, from a PCA point of view, the option center=F of prcomp() According to the help page, the calculation in prcomp() "is done by a singular value decomposition of the (centered and possibly scaled) data matrix, not by using eigen on the covariance matrix" (as it's done by princomp()) . "This is generally the preferred method for numerical accuracy"