search for: anjan

...data in column a on the y-axis and the corresponding data in column id on the x-axis, so plot(a~id). However I would like to colour these points according to the data in column b. column b data may be colour coded into the following bins: 0-9; 10-19; 20-29. Any idea on how to accomplish this? TIA, Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma 02138 phone-703.740.6939 =================================== [[alternative HTML version deleted]]

...me column members are made of alphanumeric characters: T1, T2, T3 etc. I would like to acess only that subset of the data-frame with Name == T44. X <- dataframe[dataframe$Name=='T44'] does not work. Any ideas on how to do this? I'm sure I'm missing a simple concept here. Thanks, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...nse as I have only 8 data points. The other option is to use (95*8)/100 = 7.6th data point (which can be found by interpolation between the 7th and 8th data points). Reportedly, the second formula is not too accurate. However in my case the first formula does not make much sense. Any advice? TIA, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...o Mark Fisheh of Oracle Corp. & posted this issue at OTN under Linux thread this morning. I hope that someone among you might have experienced this and can help. On that basis, I am sending this to you too. I am stuck & will really appreciate if you can shed some light on this. Thanks. Anjan *********************************************************************************************************** I have a 2 node CentOS 4.5 86_64 system (kernel 2.6.9-55.EL). On this I installed Oracle OCFS2 1.2.7-1 (with exact kernel matching). After this I installed Oracle CRS 10.2.0.1 and that in...

...gene" stores an Id; the rest of the columns store intensity data. I would like to extract the rows for gene Ids i1, i2, i3 ( I know a priori that those rows exist). So I do this: subset(d, gene %in% c(i1, i2, i3)). This does not give me the required data. Any ideas where I am going wrong? TIA, Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma 02138 phone-703.740.6939 =================================== [[alternative HTML version deleted]]

Hi, I have two vectors, each of length 45000. How do I compare the vectors to ascertain if they are identical. Secondly if they are NOT identical, how do I determine the indices of positions at which the vectors differ? Thanks, Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma 02138 phone-703.740.6939 =================================== [[alternative HTML version deleted]]

Most books on R I come across describe running statistical procedures in R. Any suggestions on a good book that teaches *programming* in R? Thanks, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...and: boxplot(dataset, col= col1) expecting to see 47 boxplots, each of a different colour. Here is the problem: the boxplots are drawn correctly but it seems that only the first few colours in col1 are being used in a repeated pattern. Anybody has any ideas on how to tackle this? Thanks in advance, Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma 02138 phone-703.740.6939 =================================== [[alternative HTML version deleted]]

Hi all, Here is my problem: I have 3 matrices , A, B, C. Each is an nXn matrix. I need to create matrix D such that : D[i,j]= (A[i,j]+B[i,j]+C[i,j])/3. Being a newbie this is proving to be a challenge. Any ideas on how best to accomplish this? Thanks! Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...y bridge seems to be Statistics-R-0.03<http://search.cpan.org/%7Ectbrown/Statistics-R/lib/Statistics/R.pm>. Would anyone like to share their experience with this Perl-R bridge? I'd like to install it in a Mac OS X. Suggestions on alternate solutions will be appreciated. Thanks in advance, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

Hi, Given that the R forum is mostly a site to discuss issues related to the language, are there any forums that discuss statistical issues? I am specifically interested in the statistical issues in experimental design. TIA, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...For example I would like to extract all rows that have time == 0h and strain == ROC. So, t <- subset(r, (r$time == "0h" && r$strain == "ROC"), select= c(time, strain, gene, deltact)) returns an empty data-frame. Is it not possible to subset based on two criteria? TIA Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma phone-703.740.6939 =================================== [[alternative HTML version deleted]]

...ave a data matrix: repeat1 repeat2 1 30 34 2 12 23 3 50 13 4 56 98 . . . I would like to plot for the mean and standard error of mean of repeat 1 and repeat2 for each row. Any pointers on how to do this in R? TIA, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...els= sampleNames(data)) produces an X-axis with 50 ticks and labels drawn from the sampleNames namespace. Problem is the labels at each tick are printed parallel to the x-axis and thus run into each other. How do I specify for the labels to be printed perpendicular to the x-axis? Thanks in advance. Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma 02138 phone-703.740.6939 =================================== [[alternative HTML version deleted]]

Hi, my data table has 38939 rows. R prints the first 11111 columns and then prints an error message:[ reached getOption("max.print") -- omitted 27821 rows ]]. is it possible to set the maxprint parameter so that R prints all the rows? tia, anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...order of the data-points to be retained and the x-coordinates of each of the 8 data-points to be equally spaced. however, using plot($time, $f) completely scrambles the ordering of the points. Any pointers to this newbie on how to plot these data points accurately will be much appreciated. Thanks, Anjan [[alternative HTML version deleted]]

hi, am new to R. what command do i use to check of certain packages have been loaded? Also, it seems that most packages come with "vignette"- is this a document that describes the usage of the package? thanks for your help, anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

I need to uninstall R 2.7.1 from my Mac. What is the best way to uninstall it? Simply delete the R icon in the Applications folder? Or is it more involved? TIA, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

...an index. I would like to calculate the mean of the numbers stored in each of the rest of the columns. So, here is my data matrix: col1 col2 col3 1 23 34 2 45 56 3 23 56 4 34 68 For each row I would like to calculate the means of the numbers stored in col2 and col3. How can this be done in R? TIA, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142 purkayas [at] wi [dot] mit [dot] edu 703.740.6939 [[alternative HTML version deleted]]

..., 0.5h, .....24h. $deltact holds real numbers. I want to subset this dataframe by the levels of $time. For all the rows which have r$time == "0h" I would like to convert the present value of r$deltact to 0. Any pointers to a simple method to accomplish would be greatly appreciated. TIA, Anjan -- =================================== anjan purkayastha, phd. research associate fas center for systems biology, harvard university 52 oxford street cambridge ma phone-703.740.6939 =================================== [[alternative HTML version deleted]]