[R] need help for t.test using "by"

I am trying to do ttest for each plate which has equal number of disease and controls. by searching this forum I found one posting suggesting OP to use "by(eo,PlateID, function(.sub) t.test(IL1Ra~Group,data=.sub))". when i modified this for my use I used to get the pvalues for each plate, recently

Re: [R] need help for t.test using "by"

any one 1Rnwb wrote: > > I am trying to do ttest for each plate which has equal number of disease > and controls. by searching this forum I found one posting suggesting OP to > use "by(eo,PlateID, function(.sub) t.test(IL1Ra~Group,data=.sub))". when i > modified this

Re: [R] Fitting a logistic regression

with my limited understanding, I am not surprised to see this data fitting nicely at the end just by eyeballing at it. the reaction at the early time point is not completed as the time passes which is close to 20 units the reaction generates more metabolite to be measured reliably your t=0 and t=1

Re: [R] need help for t.test using "by"

y"by(eo,eo$sex,function(.sub) t.test(prt1~Disease,data=.sub))" is also not working when I select the data for 40 plates. I used these two formulas earlier sucessfully. I am uploading a part of data as a delim file for sample http://www.nabble.com/file/p25051220/r..txt r..txt . Thanks 1Rnwb

[R] extracting pvalues from ttest

Hello list, I have a similar issue as this post http://tolstoy.newcastle.edu.au/R/e6/help/09/04/11438.html#options2 and I used the suggestion provided by Jorge with modifications to my data do.call(c,lapply(your_list_with_the_t_tests,function(x) x$p.value)) but I am getting the following error

Re: [R] extracting pvalues from ttest

okay fixed it by putting "c" in quote marks. 1Rnwb wrote: > > Hello, I am using "B" as a vector to store all the t.tests. since i am a > newbie to both R and statistics I am not sure if "B" is the list. Also I > see "c" used in the do.call

Re: [R] extracting pvalues from ttest

Is "B" a list? If not, I am not sure if lapply can take > it. Try aggregate(). > > On Fri, Aug 28, 2009 at 10:53 AM, 1Rnwb wrote: > >> >> Hello list, >> >> I have a similar issue as this post >> http://tolstoy.newcastle.edu.au/R/e6/help/09/0

[R] How to do adjust for sex, age, genotype for a data

Hello R gurus, I am biologist doing biomarker research and I have a data set where I have 6 proteins and close to 3000 samples, i have to look for differences between disease(Y) and controls(N) along with genetic risk, genotypes, sex and other demographic info available. however i do not know any

Re: [R] How to do adjust for sex, age, genotype for a data

then what will be the other factors needed to be adjusted and whether I should adjust or use them as covariates. Finally how these analysis will be done in R Harrell, Frank E wrote: > > 1Rnwb wrote: >> Hello R gurus, >> >> I am biologist doing biomarker research and I

Re: [R] dot plot with several points for 2 categories

I use the long way, you might have to play around with my script to get is correct. once you get it to work, you can add as many points to reflect, median, percentile etc. ## generatiing vectors for Group1 a<-LMMP8[Self_T1D=="N"] w <- a[!is.na(a)] length(w) c<-mean(a, na.rm=TRUE) c

Re: [R] heatmap plot

since heatmap is a graphic image it needs the "graphics" library. you can see these two in the examples for heatmap. require(graphics); require(grDevices) Markus Mühlbacher wrote: > > Dear R community! > > I am trying to create a heatmap based on the following data. As you can > see the diago

[R] Re gression using age and Duration of disease as a continous factors

Please explain me as what it means and how this analysis can be done using R and which library(ies) are needed. Thanks -- View this message in context: http://www.nabble.com/Regression-using-age-and-Duration-of-disease-as-a-continous-factors-tp24574133p24574133.html Sent from the R help mailing

Re: [R] Re gression using age and Duration of disease as a continous factors

01.408.8111 > > >> -Original Message- >> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- >> project.org] On Behalf Of 1Rnwb >> Sent: Monday, July 20, 2009 11:30 AM >> To: r-help@r-project.org >> Subject: [R] Re gression using age and Du

Re: [R] Re gression using age and Duration of disease as a continous factors

I thought this forum is for help. now i know what the statistician in my dept does all day long Rolf Turner-3 wrote: > > > On 21/07/2009, at 5:30 AM, 1Rnwb wrote: > >> >> Please explain me as what it means and how this analysis can be >> done using R >&g

Re: [R] Re gression using age and Duration of disease as a continous factors

the analysis. after all I do have to explain all these things to my boss as well. Thanks Steve Lianoglou-6 wrote: > > > On Jul 20, 2009, at 5:30 PM, 1Rnwb wrote: > >> I have read that multiple times without understanding anything. > > If that's the case, then per

Re: [R] R Codes for fitting Logistic Regression for Multivaraite Data

usually i go thru the manuals of the libraries and the data attached with the libraries to get the code working for my data. not a good way but usually it works out. Nimal Fernando wrote: > > Hi > > It would be a great help if someone can direct me to access some existing > R > codes for fittin

Re: [R] Antw: How to extract the upper xlim and ylim of my plot?

use hilox<-range.x=range(x variable) and hiloy<-range.y=range(y variable) or store these values in vectors to be later on user as xlim=hilox, ylim=hiloy Frank Bloos wrote: > > You may want to use the function corner.label from the plotrix-package. > > Frank > Mark Na 21.07.2009 23:03 >

Re: [R] Newbie in R: Reading .txt files and storing the 'numbers' in a vector

I am also a newbie, but here are the few thing i do eo=read.table(file="path where you have saved your text file/xxx.txt",header=T, ) #this will only read excelfiles saved as *.csv file attach(eo)

[R] characters, mathematical expressions and computed values

Hello, I have to add "Age (bar(x)=14.3) as a title on a chart. I am unable to get this to working. I have tried bquote, substitute and expression, but they are only doing a part of the job. new<- c(14.3, 18.5, 18.1, 17.7, 18, 15.9, 19.6, 17.3, 17.8, 17.5, 15.4, 16.3, 15, 17.1, 17.1, 16.4, 15.2,

Re: [R] characters, mathematical expressions and computed values

gt; From: [hidden > > email]<http://user/SendEmail.jtp?type=node&node=4645934&i=0>[mailto:[hidden > email] <http://user/SendEmail.jtp?type=node&node=4645934&i=1>] On Behalf > > Of 1Rnwb > > Sent: Thursday, October 11, 2012 2:32 PM > > To: [hidden e

[R] subtraction based on two groups in a dataframe

Hello I have a data set like below: plate.id well.id Group HYB rlt1 1 P1 A1 Control SKOV3hyb 0.190 2 P1 A2 Control SKOV3hyb 0.210 3 P1 A3 Control SKOV3hyb 0.205 4 P1 A4 Control SKOV3hyb 0.206 5 P1 A5 Control SKOV3hyb 0.

Re: [R] subtraction based on two groups in a dataframe

Thanks for the help. Sharad On Mon, Sep 27, 2010 at 9:12 PM, Remko Duursma [via R] < ml-node+2716469-935075351-6...@n4.nabble.com > wrote: > Try something like this: > > > > dfr <- read.table(textConnection("plate.id well.id Group HYB > rlt1 > 1 P1 A1 Control SKOV3hyb 0.

[R] loess with missing data points

Hello List, I have a longitudinal samples on multiple individuals, for initial analysis, i am doing some plots to see how the variable changes with the age. In some of the individuals one or two data points are missing, IL1Ra Age.at.Sample.Collection Subject.ID 6.575

Re: [R] subtraction based on two groups in a dataframe

Hi The reshape suggestions works great on my previous data, but I am unable to make is work on the new dataset. It actually works but only gives me the output of single row, instead of 96 rows. The dataset has two control groups control 1 and control 2, two disease groups viz disease 1 and disease

Re: [R] subtraction based on two groups in a dataframe

P1 B2 SKOV3hybNANANA > 7 P1 B3 SKOV3hybNANA 0.371 > 8 P1 B4 SKOV3hybNANANA > rlt1.control2 > 1NA > 2NA > 3 NA >

[R] Plot creates a straigth line

Hello all, I am using 'plot' to create standard curves for elisa data. when I use 'plot' with type='b' i get the points connected with lines and one straigth line from the lowest datapoint to the highest data point. how can i avoid/remove it from the figure. i am using R2.9.1, below is the example

[R] grep

Hi I have a script which is designed to gather data from individual columns from a file, which is an output from an instrument. the file has multiple sections and each a section has data under each column (vars), I am using the name of the column as a variable to gather the column ID using vidx<

Re: [R] Plot creates a straigth line

Thanks Jim for the explanation. I will modify the code accordingly -- View this message in context: http://r.789695.n4.nabble.com/Plot-creates-a-straigth-line-tp3004090p3005569.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-

Re: [R] exact pattern match in grep for column headers

my question is exactly as described in this post http://tolstoy.newcastle.edu.au/R/e2/help/07/02/9812.html i am doing the pattern search using a vector 'vars[vi]' where i cannot use '^vars[vi]$' for a exact pattern match -- View this message in context: http://r.789695.n4.nabble.com/grep-tp30044

Re: [R] exact pattern match in grep for column headers

Thanks I agree i did not asked the question correctly, but I solved the issue using the suggestions from this post. http://www.mail-archive.com/r-h...@stat.math.ethz.ch/msg23646.html -- View this message in context: http://r.789695.n4.nabble.com/grep-tp3004422p3006236.html Sent from the R help

Re: [R] Adding rows to column

sorry i got clik happy df2<-df1[, c(3,5,7,9,11,13,15)] df2<-df2[grep('ID', df2$Group), ] -- View this message in context: http://r.789695.n4.nabble.com/Adding-rows-to-column-tp3005607p3006302.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] exact pattern match in grep for column headers

the current code is like this vidx<-grep(paste('^',vars[vi],'$',sep=''),gsub("[[:punct:]]","",strrl1[[datbeg-1]]),ignore.case=T) where as the old one was vidx<-grep(vars[vi],gsub("[[:punct:]]","",strrl1[[datbeg-1]]),ignore.case=T) -- View this message in context: http://r.789695.n4.nabble.com/

Re: [R] Adding rows to column

If I understand correctly you want to create a new dataframe with selected columns which can be achieved this was as well, it will right away create a new dataframe with column headers df2<-df1[ ,c(3,7,9,11,13,15)] -- View this message in context: http://r.789695.n4.nabble.com/Adding-rows-to-c

[R] understanding the 4 parameter logisitc regression

I have questions regarding test=data.frame(cbind(conc=c(25000, 12500, 6250, 3125, 1513, 781, 391, 195, 97.7, 48.4, 24, 12, 6, 3, 1.5, 0.001), il10=c(330269, 216875, 104613, 51372, 26842, 13256, 7255, 3049, 1849, 743, 480, 255, 241, 128, 103, 50))) nls(log(il10)~A+(B-A)/(1+(conc/xmid )^scal),da

Re: [R] understanding the 4 parameter logisitc regression

Thanks Dieter for the help. This is how I want plot(log(test$conc),fn(test$conc,15,3.5,600,1/2.5),type="l") # looks good points(log(test$conc),log(test$il10)) regards and happy holidays sharad -- View this message in context: http://r.789695.n4.nabble.com/understanding-the-4-parameter-logis

Re: [R] Title for y-axis on right side

mtext(4, ylab="your title", font) -- View this message in context: http://r.789695.n4.nabble.com/Title-for-y-axis-on-right-side-tp3092531p3092631.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https:

[R] merging multiple columns from two dataframes

Hello, I have data in a dataframe with 139104 rows which is multiple of 96x1449. i have a phenotype file which contains the phenotype information for the 96 samples. the snp name is repeated 1449X96 samples. I haveto merge the two dataframes based on sid and sen. this is how my two dataframes loo

[R] log transformation and mean question

I have question about log2 transformation and performing mean on log2 data. I am doing analysis for ELISA data. the OD values and the concentration values for the standards were log2 transformed before performing the lm. the OD values for samples were log2 transformed and coefficients of lm were a

[R] A better way to do this

Hello gurus, I have a dataframe containing two groups viz., 'control' and 'case', each of these groups contains longitudinal data for 100 subjects. I have to plot all these subjects on a single chart and then put a regression line for each of the group for all the subjects. I have written a functi

Re: [R] A better way to do this

here is the data set, yes i am doing spagghetti plots for each Subject for MCP1 with respect to the Age.at.Sample.Collection, with a final of all the controls and all the cases Subject.ID sample Group Age.at.Sample.CollectionMCP1 19 00173-0 3455 control11.767282 2

[R] matching by gender and age

Hello R gurus, I have a data set from which i have to extract the gender and age matched rows from controls and disease group disease<-paste(rep(c('y','n'),11)) gender<-paste(rep(c('m','f'),11)) mcp<-rnorm(700,1400) age<-rnorm(32,34) dat<-data.frame(disease=disease,sex=gender,Dr_age=age[1:22],MCP=

Re: [R] matching by gender and age

I hope this modified example will work set.seed(100) disease<-paste(rep(c('y','n'),50)) gender<-c(paste(rep(c('m','f'),25)), paste(rep(c('f','m'),25))) mcp<-rnorm(100, mean=1000,sd=600) age<-rnorm(100,mean=32,sd=20) dat<-data.frame(disease=disease,sex=gender,Dr_age=age,MCP=mcp

Re: [R] matching by gender and age

Thanks, for pointing out the package e1071, the example for matchControls is exactly what I am looking for, however how can I add additional factors to match for. library(e1071) Age.case <- 40 + 5 * rnorm(50) Age.cont <- 45 + 10 * rnorm(150) Age <- c(Age.case, Age.cont) Sex.case <- sample(c("M", "

Re: [R] Clearing Console; of weeks of codes!

click in console window and type rm(list=ls(all=TRUE)), everything will be gone -- View this message in context: http://r.789695.n4.nabble.com/Clearing-Console-of-weeks-of-codes-tp3447506p3448041.html Sent from the R help mailing list archive at Nabble.com. _

[R] error for ttest

Hello all, I have arranged my data as per Dennis's suggestion in this post http://www.mail-archive.com/r-help@r-project.org/msg107156.html. the posted code works fine but when I try to apply it to my data, i get "> u2 <- ddply(xxm, .(plateid, cytokine), as.data.frame.function(f)) Error in t.test.f

Re: [R] error for ttest: solved

I knew that the "NA's" in my data were the root cause of the trouble, but did not find out how to get rid of them. Untill I found your another post mentioning to use 'na.omit' to remove the lines containing 'NA's" and the problem got fixed after that. Thanks for the help and all the trouble you hav

Re: [R] Clearing Console; of weeks of codes!

Thanks, I thought that removing the list would take care of it. the question is I do not see a .Rhistory file in my current working directory, so where it is stored. it is not visible in C:\Program files\R either. Serarching the C;\ and D:\ drives shows some old .Rhistory files but not the recent o

[R] msProcess, creating msSet file(s)

I have 120 MS chromatograms containing two columns'Time' and 'Peak height'. how do in create a msSet file which is required by all the functions of 'msProcess' to process these files. I would appreciate help very much. Thanks Sharad -- View this message in context: http://r.789695.n4.nabble.com/m

[R] splitting sample names

I have a data frame like this xx<-data.frame(cbind(Sample=c('Ctrl_6h','1+0_6h','1+200_6h','1+5k_6h','Ctrl_5K_6h','ConA_6h'), IFIT1=c(24,25,24.7,24.5,24.2,24.8))) grep('[[:digit:]]h',xx$Sample) yy<-xx$Sample strsplit(yy,"_") I have to extract the time information separated by '_

Re: [R] splitting sample names

Thanks to all of you for the suggestions and corrections. Sharad -- View this message in context: http://r.789695.n4.nabble.com/splitting-sample-names-tp3753712p3755297.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.o

Re: [R] matching by gender and age

the examples work great but it gives me this error on this file http://r.789695.n4.nabble.com/file/n3603737/il1ra_L.csv il1ra_L.csv > Age.co<-c(dat2b$Dr_Age) > Sex.co<-as.factor(dat2b$Sex) > casecont.co<-as.factor(dat2b$Self_T1D) > stat.co<-as.factor(dat2b$status) > > m <- matchControls(Sex.co~A

Re: [R] matching by gender and age

getting errors when running this file http://r.789695.n4.nabble.com/file/n3605311/il1ra_status.txt il1ra_status.txt > Age.co<-c(dat2b$Dr_Age) > Sex.co<-as.factor(dat2b$Sex) > casecont.co<-as.factor(dat2b$Self_T1D) > stat.co<-as.factor(dat2b$status) > > m <- matchControls(stat.co~Sex.co+Age.co,

[R] Help with code

hello gurus, i have a data frame like this HTN HTN_FDR Dyslipidemia CAD t1d_ptype[1:25] 1Y YY T1D 2 T1D 3 Ctrl_FDR 4 T1D 5Y

Re: [R] Help with code

I do not know how to use dput, i am attaching the txt file for the data http://r.789695.n4.nabble.com/file/n4222616/foo.txt foo.txt c1<-read.dlim('foo.txt') c2<-c1 any_comp<-NULL for( i in 1:dim(c1)[1]) { num_comp<-0 for (j in 1:dim(c1)[2]) if (c1[i,j]==2) num_comp=num_comp+1 #"Y"

Re: [R] Help with code

here is the dump and the code once again, sorry for creating so much noise. c1<-structure(list(HTN = structure(c(2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("", "Y"), class = "factor"), HTN_FDR = structure(c(2L, 1L, 1L, 1L, 2L,

Re: [R] Help with code

Thanks Bill, the output of dput was very similar to your example, I was not sure so did not put it on the post. however i uploaded the foo.txt file which contains the part of the data. sharad -- View this message in context: http://r.789695.n4.nabble.com/Help-with-code-tp4218989p4222947.html Sen

Re: [R] Help with code

this is how the ouput from the code should be structure(list(HTN = 1:10, HTN_FDR = structure(c(4L, 2L, 1L, 2L, 3L, 1L, 1L, 2L, 3L, 2L), .Label = c("Ctrl_noc", "T1D_noc", "T1D_oc", "T1d_w"), class = "factor"), Dyslipidemia = structure(c(3L, 2L, 1L, 2L, 4L, 1L, 1L, 2L, 4L, 2L), .Label = c("Ctrl_

[R] correlation matrix

Hello Gurus I have two correlation matrices 'xa' and 'xb' set.seed(100) d=cbind(x=rnorm(20)+1, x1=rnorm(20)+1, x2=rnorm(20)+1) d1=cbind(x=rnorm(20)+2, x1=rnorm(20)+2, x2=rnorm(20)+2) xa=cor(d,use='complete') xb=cor(d1,use='complete') I want to combine these two to get a third

Re: [R] correlation matrix

okay so fixed what i need to do this way finit=0 for(ri in 1:dim(xa)[1]) { finit=finit+1 xc[ri,1:finit]<-xa[ri,1:finit] xc[1:finit,ri]<-xb[1:finit,ri] } but getting error in heatmap.2 > mycol <- colorpanel(n=40,low="red",mid="white",high="blue") > heatmap.2(xc, breaks=pairs.breaks, col=mycol, Ro

Re: [R] correlation matrix

Thank you all for your suggestions. Sharad -- View this message in context: http://r.789695.n4.nabble.com/correlation-matrix-tp3891085p3894329.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://st

Re: [R] calculating ratios from all combinations

so here is the code I need help to fix the line for 'st' when it reaches maximum, that's where i am getting stuck. thanks Sharad set.seed(100) d=data.frame(x=rnorm(20)+5, x1=rnorm(20)+5, x2=rnorm(20)+5, x3=rnorm(20)+5, x4=rnorm(20)+5, x5=rnorm(20)+5, x6=rnorm(20)+5, x7=rnorm(20)+5, x8=rno

Re: [R] calculating ratios from all combinations

Thanks Micheal, that's what I wanted, I did not quite catch which variable is unused. Another part is I have my variable values in log scale so it generates '0's' instead of '1', how do i get rid of those cols. Thanks for you patience Sharad -- View this message in context: http://r.789695.n4.n

Re: [R] calculating ratios from all combinations

thanks for the help and pointer. I am modifying it like this x=which(n[,1]==n[,2]) n=n[-x,] to get rid of combinations which will generate '0' or ratio of 1. Thanks once again. sharad -- View this message in context: http://r.789695.n4.nabble.com/calculating-ratios-from-all-combinations-tp39125

[R] help with paste

Hello gurus, I have some variables, and i am creating combinations for analysis in the end i need these variables to be displayed like "LEPTIN+SAA+PTH". currently i am using loop to perform this. I would appreciate any pointers to do it without the loop. > mols=c("LEPTIN","SAA","PTH","sEGFR") > s

Re: [R] dotPlot with diagonal

is this what you are looking for library(lattice) (Seq <- matrix(c(1, 1, 6, 1, 2, 2, 5, 4, 3, 3, 4, 3, 4, 4, 3, 2, 5, 5, 2, 5, 6, 6, 1, 6), ncol = 6)) plot(Seq[1,], Seq[2,], main = "Sequenz 1 und Sequenz 2", asp = 1) abline(a=0,b=1,col='red',lwd=2) -- View this message in context: http://r.789

Re: [R] help with paste

Hi Joshua, Thanks for the suggestion. Sharad -- View this message in context: http://r.789695.n4.nabble.com/help-with-paste-tp3942818p3943753.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://s

[R] help with parallel processing code

Hello R gurus, I have the code below for which i need help and pointers to make it run in parallel on a dual core win7 computer with R 2.13.x, using foreach, iterators,doMC. library(scatterplot3d) # Loads 3D library. library(fields) library(MASS) library(ROCR) library(verification) library(caret)

Re: [R] help with parallel processing code

sorry for noise the simulated data should be like this d=data.frame(replicate(9, rnorm(40)+10),rep(c('y','n'),20)) colnames(d)<-c("LEPTIN","SAA","PTH","sEGFR","IGFBP6","MMP2","OPG","IGFBP3","PDGFAABB","group") -- View this message in context: http://r.789695.n4.nabble.com/help-with-parallel-pr

Re: [R] package for time series

if you are a beginer then try this http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm there are lots of books and their exercises available on the authors website. also search in google books http://books.google.com/books?id=0fRcUmyhUIQC&printsec=frontcover&source=gbs_ge_summary_r&c

Re: [R] help with parallel processing code

my modification gives me error > rows<- c(1:nrow(mat)) > scores <- c() > labels <-c() > itr<-1000 > chnksz<-ceiling(itr/getDoParWorkers()) > smpopt=list(chunkSize=chnksz) > foreach(icount(itr),.combine=cbind,.options.smp=smpopts)%dopar% + { + train <- sample(rows, length(rows)-1) + label =

Re: [R] help with parallel processing code

yes the library(MASS) was loaded. -- View this message in context: http://r.789695.n4.nabble.com/help-with-parallel-processing-code-tp3944303p3945972.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list ht

Re: [R] help with parallel processing code

If i understand correctly you mean to write the line as below: foreach(icount(itr),.combine=combine,.options.smp=smpopts,.packages='MASS')%dopar% -- View this message in context: http://r.789695.n4.nabble.com/help-with-parallel-processing-code-tp3944303p3945954.html Sent from the R help mailin

Re: [R] help with parallel processing code

the part of the question dawned on me now is, should I try to do the parallel processing of the full code or only the iteration part? if it is full code then I am at the complete mercy of the R help community or I giveup on this and let the computation run the serial way, which is continuing from p

[R] help with formula for clogit

I would like to know if clogit function can be used as below clogit(group~., data=dataframe) When I try to use in above format it takes a long time, I would appreciate some pointers to get multiple combinations tested. set.seed(100) d=data.frame(x=rnorm(20)+5, x1=rnorm(20)+5, x2=rnorm(20)+5,

Re: [R] row names of a heatmap on left hand side

bump -- View this message in context: http://r.789695.n4.nabble.com/row-names-of-a-heatmap-on-left-hand-side-tp4426353p4447786.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailm

[R] labeling rows in heatmap.2

how can i put the row labels on the left hand of the heatmap in heatmap.2? abnr<-structure(c(1, 0.678622097406395, 0.670294749652918, -0.0016314464654279, -0.000519068106572792, 0.199581999119988, -0.0106623494189115, 0.0840111691399559, -0.0461494399639137, 0.249279171677728, NA, 1, 0.75711406

[R] height of plots

Hello R gurus, I have to create 12 plots, I have been using the following script, which leaves a large white space between two plot. I would appreciate if someone can suggest an alternative to reduce the white space. par(mar=c(3,3,.5,.5)) split.screen(c(6,2))# split display into two scree

[R] row names of a heatmap on left hand side

Hello R Gurus, I have a matrix for which I am doing a heatmap using heatmap.2. I want to put the rownames on the lefthand side instead of the right side of the heatmap. how can i put the rownames on left hand side: I have already tried axis but could not make it work ccc<-structure(c(1, 0.283300

Re: [R] creating standard curves for ELISA analysis

9 4.5 2 8 2097693.5 16.528 13 274.5 C7 35 257 4.5 25.517 16.59 12 10 52 10 5 3 7 1456750 24 23 11 70.5 Thanks for the help 1Rnwb wrote: > > H

[R] ttest in R

Dear list, i am a biologist who needs to do some ttest between disease and non disease, sex, genotype and the serum levels of proteins on a large number of individuals. i have been using excel for a long time but it is very tedious and time consuming. i am posting the data below and ask your hel

Re: [R] ttest in R

genotype ,data = data1) > summary(aov1) > > > If you insist on t.test, here is the way: > t.test(data ~ gender ,data = data1) > t.test(data ~ disease ,data = data1) > t.test(data ~ genotype ,data = data1) > > > Cheers, > Tal > > >

Re: [R] ttest in R

since the estimation is not done pairwise so i cannot use pairwise.t.test, how do i apply tukeyHSD Thomas Petzoldt-4 wrote: > > 1Rnwb schrieb: >> Thanks for the help, but ANOVA will give me a single pvalue, then how i >> can >> make sure which group is showing the

[R] creating standard curves for ELISA analysis

Hello R guru's I am a newbie to R, In my research work I usually generate a lot of ELISA data in form of absorbance values. I ususally use Excel to calculate the concentrations of unknown, but it is too tedious and manual especially when I have 100's of files to process. I would appreciate some h