search for: purkayastha

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

...cters: 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]]

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

...ion 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]]

...//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]]

...i 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]]

...e 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, 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]]

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

...ve 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]]

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

...cks 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]]

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

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

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

Hi, I need to plot $time on the x-axis and $f on the y-axis for the following data: time f 0h 0.00 0.5h 0.54 1h 1.15 2h 2.33 4h 1.57 6h 2.19 18h 1.45 24h 1.79 I want the 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

...-7000-182 100 943 24 MJ-8000-184 100 318 25 MJ-8000-185 80 92 26 MJ-8000-186 80 161 27 MJ-8000-188 60 158 28 MJ-8000-191 60 103 29 MJ-8000-192 40 468 30 MJ-8000-193 40 432 -- ============================= 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]]

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