similar to: R environment

I want to be able to exctract a matrix from a data frame that contains repeated measurements for individuals. i.e I want to exctract the second observation for each individual in the study. Is there a way of doing this? I guess what I am trying to ask is if there as a multidimensional version of aggregate() as in aggregate(data, list(data[,"ID"]), function(x) x[2,]) but aggregate will

Hi A while back I used gsub to do the following temp<-"000US00231" gsub("something here", "", temp) "00231" I think it involved the `meta characters' somehow. I do not know how to do this anymore. I know strsplit will also work but I remember gsub was much faster. In essence the question is how to delete all characters before a particular

Dear all, I am interested in finding out how to change the names of coefficients in the lm function. I have a design matrix which I called "design" where each variate has its own name. However when I issue the command: lm.1<-lm(response~design-1, weights=some.weights) and follow it with: summary(lm.1) it seems to paste as a character the names of the variates with design i.e I

Dear list, I try to update the prabclus package. R CMD check works nicely, no warnings, good results in all tests. However, building the package fails: ginkgo:/disk5/home/chrish/RAusw/libsrc R CMD build prabclus * checking for file 'prabclus/DESCRIPTION' ... OK * preparing 'prabclus': * checking whether 'INDEX' is up-to-date ... OK * removing junk files * building

Hi, Is there a function that does percentile-percentile plot. I do not mean the qqplot. I need to plot the percentiles rather than points themselves. I am hoping for a plot that tells me that the x percentile of one data set corresponds to the y percentile of the other. for example a point on the plot of (.5, .2) will tell me that the 50th percentile of the first data and the 20th percentile of

Hi AlL, I ahve this problem that my objective function is discontinous in the paramaters and I need to use methods such as nelder-mead to get around this. My question is: How do i compute standard errors to a problem that does not have a gradient? Any literature on this is greatly appreciated. Jean,

Hi all, hope you having a nice day, I ahve this weird results with identical (probably I am not understanding correctly what it does ...) I have these two data frames and I issue : > identical(temp, temp1) [1] FALSE However, these data frames are Nx2 and when I issue: > identical(temp[,2], temp1[,2]) [1] TRUE > identical(temp[,1], temp1[,1]) [1] TRUE and the results from str >

Dear R users, I am trying to install R 1.7.1 on a sparc-sun-solaris2.8 system. although gnu readline is installed and works fine on the parent directory, R is not recognizing it. I get the following in the config.log file. configure:11627: checking for rl_callback_read_char in -lreadline configure:11658: gcc -o conftest -g -O2 -L/usr/local/lib conftest.c -lreadline -ldl -ltermcap -lm

I have this awkward problem with trellis (lattice). I am trying to generate some plots through loops but the .eps file is empty. When I generate them in a list and print them outside the loop all is fine. this is an example below:( nothing shows up in foo.eps, but all show up in foo1.eps) R vesion 2.0.1, lattice version 0.10-16, on a debian 2.6.8-1 kernel. X <- data.frame(x=rnorm(10000),

Is there a dataload utility for linux. The link in genstat is down but I managed to find the utility at: http://gurukul.ucc.american.edu/econ/gaussres/UTILITYS/DATALOAD.HTM but this is a dos/windows version. Thank you Jean

Dear All, I have the following data.frame `data.frame': 1563 obs. of 4 variables: $ Model :Class 'AsIs' chr [1:40] "Astro" "Astro" "Astro" "Astro" $ Make :Class 'AsIs' chr [1:40] "Chevrolet" "Chevrolet" $ Production : num 11219 12384 1082 5409 5458 ... $ date :`POSIXlt', format: chr

Hi There I have a lm object with 4 parameters and I want to test wether 2 parameters are equal using a Wald test (basically b1=b2 or b1-b2 =0). In the help file from R it says that under ANOVA the optional arguments " Terms" or "L" test whether a linear combination is equal to 0. I tried; >anova(m1, Terms = Beta1-Beta2=0) but I get the error: Object " Beta1"

Dear All, I am trying to solve a Generalized Method of Moments problem which necessitate the gradient of moments computation to get the standard errors of estimates. I know optim does not output the gradient, but I can use numericDeriv to get that. My question is: is this the best function to do this? Thank you Jean,

Dear list, I have been struggling for some time now with this code... I have this vector of unique ID "EID" of length 821 extracted from one of my dataframe (skate). It looks like this: > head(skate$EID) [1] "896-19" "895-8" "899-1" "899-5" "899-8" "895-7" I would like to remove the complete rows in another dataframe

Dear All, I may look ridiculous, but I am puzzled at the behavior of the nls with a fitting I am currently dealing with. My data are: x N 1 346.4102 145.428256 2 447.2136 169.530634 3 570.0877 144.081627 4 721.1103 106.363316 5 894.4272 130.390552 6 1264.9111 36.727069 7 1788.8544 52.848587 8 2449.4897 25.128742 9 3464.1016 7.531766 10 4472.1360 8.827367 11

Dear list, I have this dataframe: > names(events) [1] "EID" "X" "Y" "trip" "tow" "catch" "effort" "depth" [9] "season" Where some of my unique ID "EID" appears more than once in 162 cases. > length(events$EID)-length(unique(events$EID)) [1] 162 I would like to combined

Hi, I'm trying to get a toy program making use of nmmin to run successfully. I've gotten to the point of compiling. However, when I attempt to run my executable, I guess a bus error. I see that someone else has asked about using nmmin before <http://tolstoy.newcastle.edu.au/R/help/06/03/23944.html>, but I haven't come across any replies. Is there some documentation on how to

Hi all, Since roxygen is a great help to document R packages, I am wondering if there exists an approach to go back from the raw Rd files to roxygen-documentation? E.g. turn "\author{Somebody}" into "@author Somebody". This sounds ridiculous, but I believe it helps in the long term for me to maintain R packages. Thanks! Regards, Yihui -- Yihui Xie <xieyihui at

The current paravirt startup_ipi hook for vmware commit: ae5da273fe3352febd38658d8d34484cbcfb3423 is quite frankly ridiculous. In the middle of wake_up_secondary_cpu: We have: /* * Paravirt / VMI wants a startup IPI hook here to set up the * target processor state. */ startup_ipi_hook(phys_apicid, (unsigned long) start_secondary, (unsigned

The current paravirt startup_ipi hook for vmware commit: ae5da273fe3352febd38658d8d34484cbcfb3423 is quite frankly ridiculous. In the middle of wake_up_secondary_cpu: We have: /* * Paravirt / VMI wants a startup IPI hook here to set up the * target processor state. */ startup_ipi_hook(phys_apicid, (unsigned long) start_secondary, (unsigned