similar to: Printing multiple graphics devices

Hi I have a column of concatenated information stored in an RG object in the limma package and I need to split this information and then paste the first two pieces of data in each case back into two columns of the RG object. This is how I am currently doing this gene.info.split<-strsplit(RG$genes$Name,",",fixed=TRUE) for (h in 1: length(gene.info.split)){

I need to specify error terms properly in a mixed-effects anova model. I know you can add error terms in aov using Error and can specify random factors in lme but I am not sure how these get treated. When making the calculations for fixed and random factors, are the correct error terms used and how can you get aov or lme to use different error terms for fixed and random effects? I'm

I am using contour to draw contour lines onto a photo (jpeg) of a leaf on a white background. I have two problems: 1) The contour gets plotted at right angles to the jpeg image. I guess this is a problem of referencing the start and end points of the image matrix but I can't see how to over come this other than loading in a second image that has been rotated 90 degrees and plotting the

Dear R-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to plot data. I have tried emailing a MAANOVA discussion group, as well as the author of the package, and have not yet received a response so I am hoping that someone on this listserv can be of assistance. There are several functions in MAANOVA (riplot, resiplot) which call the

I'm trying to install packages on windows XP and I have trouble with command Rcmd build (R version 1.8.1) : In the Windows console for package maanova for example, answer is : C:\Documents and Settings\dillies\Mes documents\ghis\packages>Rcmd build maanova * checking for file 'maanova/DESCRIPTION' ... OK * preparing 'maanova': * cleaning src *

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

Hi, I run a maanova analysis and found this message error: Error in ma.svd(X, 0, 0) : 0 extent dimensions I did a google search and found this: \item ma.svd: function to compute the sigular-value decomposition of a rectangular matrix by using LAPACK routines DEGSVD AND ZGESVD. \item fdr: function to calculate the adjusted P values for FDR control. I did a search for LAPACK and

Hello, R community, I have been using the lmer and mcmcsamp functions in R with some difficulty. I do not believe this is my code or data, however, because my attempts to use the sample code and 'sleepstudy' data provided with the lme4 packaged (and used on several R-Wiki pages) do not return the same results as those indicated in the help pages. For instance: > sessionInfo() R

Hello, I am basically familiar with R and am trying to import a module that someone else has written. I know that it must go into the R library but even after I place the file there R doesn't recognize it. The module is maanova, available from the Churchill lab group for analysis of microarray data, if anyone is familiar with it. Any ideas/help? Thanks, Heather --

Hi, I am trying to run an analysis with the package maanova and I am not getting success. I suppose that I am wrong on set up the formula, so the issue may not be related to R, properly. I have two varieties of plants (V1 and V2). A group of each ones were treated and another was not treated. After treatment, in three different time RNA was collected from treated and from not treated plants for

Dear R-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to get interactions. I am a newbie in both, R and maanova, and I do not have good knowledge in statistical methods. I have four effects: Effects Levels Var 2 Ind 2 Trat 2 Time 4 Sample 3 <-- biological replicate Spot 4 <-- technical replicate I had

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

Generator Microsoft Word 11 (filtered medium) Hi, I'm the mannova package maintainer. We used La.svd(X, method="dgesvd") in maanova package before. After R-2.3.0, the old La.svd() method was deprecated for option method="dgesvd". I changed maanova code correspondingly, which will call method="dgesdd" instead. But after that, we keep getting below error message

Hi all, It seems that setting userMask in the doveadm http api's "user" command to anything involving wildcards ?s or *s causes the API to fail. When using the API with a userMask set to the fully qualified email address, I get the expected fields back. From a python shell, with an asterisk causes the following behaviour (note that doveadm is a python module that wrappers the

Hello, unfortunately I have I big problem I can't solve. I have to analyse if a gene is tissue specific. For example for the gene xyz I have following expression values: Heart Liver Brain 8.998497 10.013561 12.277407 9.743556 10.137574 11.033957 For every tissue I have two values from two different experiments. Now I want to test if Heart is significant higher

Hi, Am using maanova package for doing anova.But am getting error like this..plz, help me regarding this.. > TGR=read.madata("rmaexpr.dat",designfile="design.dat") Reading one color array. Otherwise change arrayType='twoColor' then read the data again Warning messages: 1: In read.madata("rmaexpr.dat", designfile = "design.dat") : Assume that

Hi, I have taken one microarray experiment and trying to implement same statistical measures what they have done.I have taken datasets from GEO platform with accession number GSE1557. In the experiment,about half of double transgenic rats (dTGR) over-expressing the human renin and angiotensinogen genes die by age 7 weeks of terminal heart failure (THF); the other (preterminal) half develops

Hi all, I am having memory allocation problem with my R 2.2.1 for Mac OS. The following is the error message that I get. I do not get this message if I break down the large dataset in to sub datasets. I think breaking up the dataset is not a sustainable solution in the long run. The data that I am analysing is essentially big, and it would be reasonable to do the analyis on the whole dataset

Hi, I'm trying to fit a model using variance components in R, but if very new on it, so I'm asking for your help. I have imported the SPSS database onto R, but I don't know how to convert the commands... the SPSS commands I'm trying to convert are: VARCOMP RATING BY CHAIN SECTOR RESP ASPECT ITEM /RANDOM = CHAIN SECTOR RESP ASPECT ITEM /METHOD = MINQUE (1) /DESIGN

Dear Help I'm new to R (from matlab)...using windows XP. I've tried to work out, to no avail, 4 things: 1) rotating the numbers on axes...something to do with par(str) ?? 2) how to close a window having opened one e.g. windows(7,7) 3) how to manipulate the key (e.g. dots, lines etc) on the legend. Using pch just gives me the same key for all functions on a plot. i.e. legend