search for: szhan

Hello, R experts: I got data like this: group duplicate treatment A Y 5 A Y 3 A N 6 B Y 2 B N 4 B Y 1 How to sort the data and calculate the average treatment value for each group in two level of duplicate. Results like this: group duplicate treatment A Y 4 A N

Hello, R experts, I tried to retrieve all biological process GO terms at level 3 starting "biological process" as level 1 using the code as bellows: 1 library(GO) 2 library(GOstats) 3 level2<-getGOChildren("GO:0008150")$"GO:0008150"$Children 4 for ( i in 1:length(level2)) { 5 level3 <- getGOChildren(level2[i])$level2[i]$Children 6 for ( j in

Dear R experts: I tried two ways to install Package:XML on windows xp for R 1.9.1, all failed. Messages are given as bellows: 1> download from CRAN > install.packages("XML", CRAN = getOption("CRAN"), + contriburl = contrib.url("http://cran.r-project.org"), + available = NULL, destdir = NULL, + installWithVers = FALSE) trying URL

Dear R experts: Can you help me to overlap a histogram and theoretical density curve of poison distribution? for example data like this: The numbers of sanils found in each of 100 sampling quadrats in an area were as follows: number of snails, r 0 1 2 3 4 5 8 15 f=frequency of r 69 18 7 2 1 1 1 1 apparently this is not Poison but near to Poison, How to overlap the theoretical desity curve to

Dear R-users: I am runing R 1.6.2 with Windows XP. I try to calculate Pearson correlation and Spearman correlation of any pairwise columns of 8000 x 80 data matrix with missing values and randomize the matix 1000 times and calculate this two correlations again. The code bellow for Pearson is working fine but for Spearman got the error bellow for randomized data matrix and R console is stop

Hello R experts, I'd like to studentize a matrix (tmp1) by column using apply() function and skip some specific values such as zeros in the example below to tmp2 but not tmp3. I used the script below and only can get a matrix tmp3. Could you please help me to studentize the matrix (tmp1) without changing the zeros and generate a new matrix tmp2? Thanks, Joshua tmp1 [,1] [,2] [,3] [,4]

Hello R users, I have gene expression data of two groups of genes (large and small). Gene expression intensities of those genes are classified into 1 to 10 levels. What I want is to make a random set of genes that have the same levels as the small group from large group using sample(). I used smallvec to hold the number of genes in each levels (1 to 10) for small group, largevec for large group.

Dear R users: Is there a function or easier way to randomize the data between rows (considering a row as a whole unit)? For example, original data: A1 A2 A3 A4 B1 1 2 3 4 B2 5 6 7 8 B3 9 10 11 12 B4 13 14 15 16 randomized data: A1 A2 A3 A4 B4 13 14 15 16 B2 5 6 7 8 B1 1 2 3 4 B3 9 10 11 12 Thank you for your attention! Joshua

Hello, R users, I have a long list with over 20,000 elements (say foo), now I want to make a sub list of 30 elements from this list using the component names such as subname<-c( "a", "b", "d", "f", "g", "h", "o", "q",...). One way to do it is to make rest of elements in the list NULL, that is too silly. Is there

Hello, R experts, Sorry for asking this question again again since I really want a help! I have a two-factor experiment data and like to calculate estimates of interation contrasts say factor A has levels of a1, a2, and B has levels of b1, b2, b3, b4, and b5 with 3 replicates. I am not sure the constrast estimate I got is right using the script below:

Hello, R experts, I have a list called dataHP which has 30 elements (m1, m2, ..., m30). Each element is a 7x6 matrix holding yield data from two factors experimental design, with treatment in column, position in row. For instance, the element 20 is: dataHP[[20]] col1 col2 col3 trt1 trt2 trt3 [1,] 22.0 20.3 29.7 63.3 78.5 76.4 [2,]

...`.' 0.1 ` ' 1 Residual standard error: 0.2874 on 44 degrees of freedom Multiple R-Squared: 0.973, Adjusted R-squared: 0.9601 F-statistic: 75.55 on 21 and 44 DF, p-value: < 2.2e-16 Joshua Quoting "(Ted Harding)" <Ted.Harding at nessie.mcc.ac.uk>: > On 24-Aug-06 szhan at uoguelph.ca wrote: >> Hello, R users, >> I have two factors (treat, section) anova design experiment where >> there are 3 replicates. The objective of the experiment is to test if >> there is significant difference of yield between top (section 9 to 11) >> and botto...

Hello, R users, I have two factors (treat, section) anova design experiment where there are 3 replicates. The objective of the experiment is to test if there is significant difference of yield between top (section 9 to 11) and bottom (section 9 to 11) of the fruit tree under treatment. I found that there are interaction between two factors. I wonder if I can contrast means from levels of

...nandes en bpi.pt) > 46. connecting to running process possible? (Charles Cosse) > 47. Re: Overlaying lattice graphs (continued) (hadley wickham) > 48. RServe (java2R) question (Guanrao Chen) > 49. Re: Lattice: hiding only some strips (Deepayan Sarkar) > 50. interaction contrast (szhan en uoguelph.ca) > 51. Re: Barchart legend position (Deepayan Sarkar) > 52. How to run "mathematica" or "c" programs in R? (Zhang Jian) > 53. Re: Imputing missing values in time series (Horace Tso) > 54. Asteriscs in a plot to represent approximate size of p-value...

Hello, There: I read data from tab-delimted text file(1888.txt) in C:/temp, which has thousands rows and 80 colulms, using read.delim("C:/temp/1888.txt"). when I retrieved the numeric coulum which has real zero values and applies R buid-in function such as mean( ), sum() and got a result of NA. It seems that R (R 1.6.2 in windows) considers zero as missing value (NA). Thus my

Hello, Everyone, I am sorry that my message got truncated. I resend it again as below: Hello, R Users, I have a list (say listexp) of 10,000 elements, each of which consists of a matrix (5X6). It likes: $"a" trt1rep1 trt1rep2 trt2rep1 trt2rep2 ctlrep1 ctlrep2 [1,] 50 54 98 89 40 45 [2,] 60 65 76 79

Hello, R experts, Sorry for asking this question again since I really want a help! I have a two-factor experiment data and like to calculate estimates of interation contrasts say factor A has levels of a1, a2, and B has levels of b1, b2, b3, b4, and b5 with 3 replicates. I am not sure the constrast estimate I got is right using the script below:

Hello, R experts, I have a two-factor experiment data and like to calculate estimates of interation contrasts say factor A has levels of a1, a2, and B has levels of b1, b2, b3, b4, and b5 with 3 replicates. I am not sure the constrast estimate I got is right using the script below: score<-c(7.2,6.5,6.9,6.4,6.9,6.1,6.9,5.3,7.2,5.7,5.1,5.9,7.6,6.9,6.8,

Hello, R experts, Sorry for asking this question again since I really want a help! I have a two-factor experiment data and like to calculate estimates of interation contrasts say factor A has levels of a1, a2, and B has levels of b1, b2, b3, b4, and b5 with 3 replicates. I am not sure the constrast estimate I got is right using the script below: