similar to: sorting the data~

here is the data: y<-c(5,2,3,7,9,0,1,4,5) id<-c(1,1,6,6,7,8,15,15,19) t<-c(50,56,50,56,50,50,50,60,50) table1<-data.frame(y,id,t)//longitudinal data what I want to do is to use the linear model for each id ,then get the estimate value,like: fit1<-lm(y~t,data=table1,subset=(id==1)) but ,you can see the variable "id" is quite irregular,they are not arranaged in order

Hi: I need to create many variables at one time,how to do this in R? for eg ,X1,X2.......X100? Thanks~ [[alternative HTML version deleted]]

Dear all, I have a date frame like this: X Y Z 22 24 4.3 2.3 3.4 5.3 ..... .... 57.2 23.4 34 What my purpose is: to sort the data frame by either X, Y or Z. sample output is (sorted by X) : X Y Z 2.3 3.4 5.3 ..... ...... 22 24 4.3 ... 57.2 23.4 34 I have no idea how to use sort, order or rank functions. Please help me out. Thanks! Leaf

I need to create 100 variable ,whose name is id.1,id.2~~~~id.100 then I need to let a vector say id<-c(id.1,id.2....id.100) any easy way to do this? thanks a lot~ [[alternative HTML version deleted]]

Hi everyone: I have a density function f(x|theta)=theta*x^(theta-1),where 0<x<1,0<theta<infinite I want to pratice on R to find the MLE of this function,here is my code: x <- (0:10)/10 f<-function(theta) prod(theta*x^(theta-1)) mle(f) and r gave me :Error in eval(expr, envir, enclos) : argument is missing, with no default what mistake I just made?and how to add a

hi: just a very simple quesion,how to do multiple multiplication in R x<-c(1,2,3,4,5,6) how to get 1*2*3*4*5*6? I checked the Arithmetic Operators in R,but did not found the operators for this function,anyway can do this except the loop? [[alternative HTML version deleted]]

Dear all, I would like to sort the values of "member" for each group, and obtain a variable to indicate its rank within the group. For example, we have original dataset as follows: df <- data.frame(group = c(rep(1, 3), rep(5, 3)), member = c(30, 10, 22, 21, 44, 15)) group member 1 1 30 2 1 10 3 1 22 4 5 21 5 5 44 6 5 15 After

a table with three varialbe:treatment (low,high),size(large,small),response y I want all response y with treatment=low and size=large, I try to write as : treatment<-c("low","low","high","high") size<-c("small","large","small","large") y<-c(1,2,4,5) table1<-data.frame(treatment,size,y)

here is the data: y<-c(5,2,3,7,9,0,1,4,5) id<-c(1,1,6,6,7,8,15,15,19) t<-c(50,56,50,56,50,50,50,60,50) table1<-data.frame(y,id,t)//longitudinal data the above is only part of data. what I want to do is to use the linear model for each id ,then get the estimate value,like: fit1<-lm(y~t,data=table1,subset=(id==1)) but ,you can see the variable "id" is quite

Hi everyone: I try to use r to do the Cholesky Decomposition,which is A=LDL',so far I only found how to decomposite A in to LL' by using chol(A),the function Cholesky(A) doesnt work,any one know other command to decomposte A in to LDL' My r code is: library(Matrix) A=matrix(c(1,1,1,1,5,5,1,5,14),nrow=3) > chol(A) [,1] [,2] [,3] [1,] 1 1 1 [2,] 0 2 2

hello, I'd like know how to sort a data frame in R for example how I should do to sort by Catholic with swiss data frame like below thanks Fertility Agriculture Examination Education Catholic Infant.Mortality Courtelary 80.2 17.0 15 12 9.96 22.2 Delemont 83.1 45.1 6 9 84.84 22.2

Hello, I am not only interested in finding out which genes are the most highly up- or down-regulated (which I have done using the linear models and Bayesian statistics in Limma), but I also want to know which genes are consistently highly transcribed (ie. they have a high intensity in the channel of interest eg. Cy5 or Cy3 across the set of experiments). I might have missed a straight forward

Hi all, I am attempting to run a script in which I permute my data and run a Wilcoxon rank sum test on the data 1000 times and compare my original test statistic to the permuted test statistics to more accurately estimate a significance level for the trends I am observing. In the process of doing this, I need to sort a data frame that measures 3 x 1001 based on the value of the test statistic

I've noticed that a frequently asked question on R-help is how to sort a data frame by multiple columns. Since this question is asked so often, making this task easier for users seems a worthwhile goal. At a minimum, the following changes to the documentation would surely be helpful and reduce email postings to R-help: 1. In the "See Also" section of the help page for

Hi everyone: I have been struggling with this repeated data type for whole afternoon,I sent two emails to server for help,many people kindly responded , hereby thank you so much,but since I dont want to write to much in email,so I divide the problem in parts,so far this seem did not work out very well,so this is my whole problem~ first I have example of data here:

I need to sort users by their rank how do I do this? -- Posted via http://www.ruby-forum.com/.

Derar R-people I have troubles understanding what the function sort.list() is doing. On the homepage it says that it returns a permutation which rearranges a vector into ascending or descending order (like order() but on a vector instead of a sequence). > sort.list(c(0, 2, 10, 11, 4)) [1] 1 2 5 3 4 which does not make sense to me. In fact I am getting the same (non-sensical) result using

Hi, I recently ran into the problem that I needed a Siegel-Tukey test for equal variability based on ranks. Maybe there is a package that has it implemented, but I could not find it. So I programmed an R function to do it. The Siegel-Tukey test requires to recode the ranks so that they express variability rather than ascending order. This is essentially what the code further below does. After the

Hello, First time posting to this mail list. I'd like to use R in the most efficient way. I'm accomplishing what I want, but feel there is a more R'ish way to do it. Just learning R. *My goal: get ranks of value across rows with row names and column names intact.* I'm guessing one of the [?]apply functions will do what I need, but I couldn't sort out which one (after a lot

Hi, Is there a function like rank but that solves the ties by randomly assigning a value (doesn't average ranks of ties). This is what I actually need: I want to make NA all elements of each column in an array that are ranked in a position larger that rankmax for each column. # Say I've got an array b: b<-cbind(c(1:5,5:1),c(1,12,14,2,5,4:8)) #> b # [,1] [,2] #[1,] 1 1