search for: amitrhelp

--- On Fri, 15/5/09, Amit Patel <amitrhelp at yahoo.co.uk> wrote: > From: Amit Patel <amitrhelp at yahoo.co.uk> > Subject: Help with loops > To: r-help at r-project.org > Date: Friday, 15 May, 2009, 12:17 PM > Hi > I am trying to create a loop which averages replicates in > my data. > The original data has...

Does anyone know an easy way to convert all the zero values in a imported csv table into NA's [[alternative HTML version deleted]]

Hi I needed some help with ANOVA I have a problem with My ANOVA analysis. I have a dataset with a known ANOVA p-value, however I can not seem to re-create it in R. I have created a list (zzzanova) which contains 1)Intensity Values 2)Group Number (6 Different Groups) 3)Sample Number (54 different samples) this is created by the script in Appendix 1 I then conduct ANOVA with the command >

Sorry i had a misprint in the appendix code in the last email Hi I needed some help with ANOVA I have a problem with My ANOVA analysis. I have a dataset with a known ANOVA p-value, however I can not seem to re-create it in R. I have created a list (zzzanova) which contains 1)Intensity Values 2)Group Number (6 Different Groups) 3)Sample Number (54 different samples) this is created by the

Hi I have found a little problem with an R script. I am trying to merge some data and am finding something unusual going on. As shown below I am trying to assign (MatchedValues[Value2,Value]) to (ClusteredData[k,Value]) which are two separate dataframes. 1) By the following command you can see that the value im transferring is 481844.03 > MatchedValues[Value2,Value] [1] 481844.03 6618

Hi I have tried many attempts but cant get the loop right, as I am not a strong programmer. What I am basically trying to do is compare 2 spreadsheets. The problem is that one of them only contain a portion of the overall data (TESTSAMP), where the other has a full datasetFULLSAMP. From the complete set I would like to remove the rows of data which are not in the TESTSAMP. Column 1 contains the

Hi I am trying to create a column in a data frame which gives a sigificane score from 0-7. It should read values from 7 different colums and add 1 to the counter if the value is <=0.05. I get an error message saying Error in if (ALLRESULTS[i, 16] <= 0.05) significance_count = significance_count + : missing value where TRUE/FALSE needed The script is included below it works if i

Hi i am using a script which involves the following loop. It attempts to reduce a data frame(zz) of 95000 * 41 down to a data frame (averagedreplicates) of 95000 * 21 by averaging the replicate values as you can see in the script below. This script however is very slow (2days). Any suggestions to speed it up. NB I have also tried using rowMeans rather than adding the 2 values and dividing by 2.

When doing the t-test in the below manner will r compare each element of the array with the relevant one. I.e. if i was comparing x and y would (1 and 0) and (1 and 9) be treated as separate variables. Or does it just assume one variable. # test data x <- c(1,1.1,1.15,1.2,1.21,1.23) y <- c(0.9,1,1.16,1.18,1.19,1.2) z <- c(1.4,1.42,1.43,1.44,1.45,1.46) ###  Student's t-test #

Hi I have done ANOVA on a dataset (See Below) but am having problems retreiving the p-value. I am assuming that Pr(>F) is the p-value but cannot get this value or in fact any other value (e.g. DF) from the summary.Any suggestions?? I have tried sum<-summary(zzz.aov) > sum$Pr(>F) Error: unexpected '>' in "sum$Pr(>" > sum$"Pr(>F)" NULL >

Hi I am attempting to use plsr which is part of the pls package in r. I amconducting analysis on datasets to identify which proteins/peptides are responsible for the variance between sample groups (Biomarker Spoting) in a multivariate fashion. I have a dataset in R called "FullDataListTrans". as you can see below the structure of the data is 40 different rows representing a

Hi I am analysing a dataset of 40 samples each with 90,000 intensity measures for various peptides. I am trying to identify the Biomarkers (i.e. most significant peptides). I beleive that PLS with jack knifing, or alternativeley CMV(cross-model-validation) are multivariateThe 40 samples belong to four different groups. I have managed to conduct the plsr using the commands: BHPLS1 <-

Hi When I call for the loadings of my plsr using the command, x <- loadings(BHPLS1) my loadings contain variable names rather than numbers. >str(x) loadings [1:94727, 1:10] -0.00113 -0.03001 -0.00059 -0.00734 -0.02969 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:94727] "PCIList1" "PCIList2" "PCIList3" "PCIList4" ... ..$ : chr

Hi I am attempting to do a loadings plot from a plsr object. I have managed to do this using the gasoline data that comes with the pls package. However when I conduct this on my dataset i get the following error message. >plot(BHPLS1, "loadings", comps = 1:2, legendpos = "topleft", labels = "numbers", >xlab = "nm") Error in

Hi I have implented boxplots in my script to create box plots BoxplotsCheck <- readline(prompt = "Would you like to create boxplots for any Feature? (y/n):")   if (BoxplotsCheck  == "y"){     BoxplotsFeature <- readline(prompt = "Which Feature would you like to create a Boxplot for?:")     BoxplotsFeature <- as.numeric(BoxplotsFeature)     BoxplotsData

Hi I have got a long script which will not run for me as i keep getting errors : > source("clusterfixV1_4.r") Error in source("clusterfixV1_4.r") : clusterfixV1_4.r: unexpected symbol at 158: eck[k,2] <- as.numeric(1) 159: #ClusterInfo[k,2] <- "Clustered I have sorted all the ones i can but i am having a problem here Can anyone tell me the

Hi I am trying to create a loop which averages replicates in my data. The original data has many rows. and consists of 40 column zz[,2:41] plus row headings in zz[,1] I am trying to average each set of values (i.e. zz[1,2:3] averaged and placed in average_value[1,2] and so on. below is my script but it seems to be stuck in an endless loop Any suggestions?? for (i in 1:length(average_value[,1])) {

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Paul         __________________________________________________ [[alternative HTML

Hi I am trying to conduct a Bartlett test between two groups Samp 1 and Samp 2, both of which are vectors of equal length. I cant find any information on how to do this. Does the data need to be in a structured list. Thanks in advance

Hi I Have a very large dataset that I would like to conduct ANOVA tests on. Im not a very strong programmer so any help would be appreciated. the format is Identifier             A1       A2        B1      B2       C1   C2      Norm1         Norm2 1234                  1        1            NA     NA      4       3        NA               NA 4567                  2        2              4