similar to: aov summary to matrix

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

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 <-

I would appreciate it if someone could explain the results of the example from the aov() help file. The output given below is different from book Venables and Ripley - MASS The results In R1.5.1 (Under windows) is as follows: > 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 <-

Hi, I am working on table suing R and Latex. I am writing .Rnw file first in which i m reading input file and storing into dataframe. After filtering certain values from this dataframe. I am planning to display it. I don't want to use xtable since i need to change column names. *Sample .Rnw file* @ <<echo=FALSE>>= cat("\begin{table}[ht]")

Does anyone know how to obtain a plot of residuals by predicted values for a main-effects aov? I want to check that the residuals are distributed equally across treatment means. Thanks, Jeff [[alternative HTML version deleted]]

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr <- ts(t(tkr), deltat=1/12) or tstkr <- ts(c(tkr), deltat=1/12) and use the stl function I get the following error: Error in

Hi, I am trying to find the best way to read 85 tick data files of format: > head(nbbo) 1 bid CON 09:30:00.722 09:30:00.722 32.71 98 2 ask CON 09:30:00.782 09:30:00.810 33.14 300 3 ask CON 09:30:00.809 09:30:00.810 33.14 414 4 bid CON 09:30:00.783 09:30:00.810 33.06 200 Each file has between 100,000 to 300,300 rows. Currently doing nbbo.list<-

The build system has rolled up R-2.15.2.tar.gz (codename "Trick or Treat") at 9:00 this morning. This is a maintenance release; see the list below for details. You can get it from http://cran.r-project.org/src/base/R-2/R-2.15.2.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

The build system has rolled up R-2.15.2.tar.gz (codename "Trick or Treat") at 9:00 this morning. This is a maintenance release; see the list below for details. You can get it from http://cran.r-project.org/src/base/R-2/R-2.15.2.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

The build system has rolled up R-2.15.2.tar.gz (codename "Trick or Treat") at 9:00 this morning. This is a maintenance release; see the list below for details. You can get it from http://cran.r-project.org/src/base/R-2/R-2.15.2.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

Hello, We are planning to moving most of our servers to ESX but before buying our SAN, we want to do some I/O stats to see if iSCSI is enough or if we have to go with FC. So I found a plugin for Nagios that can log I/O stats with iostat. So far it's fine with single disk/one partition servers, but on our Oracle Database 10g server, we have two drives in RAID 1 (/dev/sda) and 4 other

Hi, I think that it would be a good thing to know more about those 2 projects (and also the future patent free format). I think that many people as me know about Lame, but not about Vorbis, and vice-versa. It would be fine that someone (perhaps the maintainer) of every project would introduce to both group of people those projects. 2 things would be interesting (to my mind): - to know about the

I continue to struggle with mixed models. The square zero version of the problem that I am trying to deal with is as follows: A number (240) of students are measured (tested; for reading comprehension) on 6 separate occasions. Initially (square zero) I want to treat the test time as a factor (with 6 levels). The students are of course ``random effects''. Later I want to look at

Hi folks: I've recently migrated one of my production servers (bad Jason!) from WBEL to CentOS 3.3 using Donavan's instructions in the FAQ (http://www.centos.org/modules/smartfaq/faq.php?faqid=19). While it was as painless as the document suggests, I've run into problems with postgrey (a postfix greylisting daemon) on my server. Whenever I try to start postgrey, I get the

Greetings, I downloaded R for my mac from software, R binaries, rm123, base and then downloaded "rm123.sit.hqx" Now, when I open R and try to type anything in the blank, I get a message that I need CarbonLib 1.20 or newer to run R. Please tell me what this means and where to go to get it. By the way, when I type anything in and recieve the error message, R will then refuse to close

Hello, I have a small problem with developing design matrix X, which I use in estimation the log-likelihood of a multinomial logit model. I have the data: number of observation - 289 number of choice alternative- 3 number of choice specific variables in matrix X -4 matrix X =289x4 I tried to use the function createX, I know that I have to get design matrix 289x12 (am I right?) but

CentOS Errata and Enhancement Advisory 2019:1653 Upstream details at : https://access.redhat.com/errata/RHEA-2019:1653 The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 6110d3f449df1361f3ec449ef3763e808bee85252a7d1aa40f453ad00115dbd9 microcode_ctl-1.17-33.14.el6_10.i686.rpm x86_64:

Dear All, I am working on a large matrix of dimension 20000x700 say 'mat'. I have calculated pearson correlation for the rows of the matrix and their p-values using rcorr function in library Hmisc. Now I wish to filter out those pairs who's PCC value is above 0.8 cut off and p-value is less than 0.05. >library(Hmisc) >mat_cor=rcorr(t(mat),type="pearson")

Send CentOS-announce mailing list submissions to centos-announce at centos.org To subscribe or unsubscribe via the World Wide Web, visit https://lists.centos.org/mailman/listinfo/centos-announce or, via email, send a message with subject or body 'help' to centos-announce-request at centos.org You can reach the person managing the list at centos-announce-owner at centos.org When

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =