search for: jeaneid

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

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

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

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

...m1, Terms = Beta1-Beta2=0) but I get the error: Object " Beta1" must be assigned locally before replacement. I also tried >anova (m1, Terms = 1-2 = 0) and I get: Invalid assginment: No object name : 1-2 = 0 What am I doing wrong? Thanks for your help!! >From: Jean Eid <jeaneid at chass.utoronto.ca> >To: r-help at stat.math.ethz.ch >Subject: Re: [R] ksmooth in SPLUS vs R Date: Mon, 22 Sep 2003 20:23:32 >-0400 > >Please do forget my questions as they are really trivial and I do not kno >what I was thinking of. > >Thank you thomas for clarifying...

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

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

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 >

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),

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,

...ympotic standard errors for your variance components estimates using the delta method. J.R. Lockwood 412-683-2300 x4941 lockwood at rand.org http://www.rand.org/methodology/stat/members/lockwood/ ------------------------------ Message: 20 Date: Thu, 26 Jun 2003 10:58:01 -0400 From: Jean Eid <jeaneid at chass.utoronto.ca> Subject: [R] lm diagnostics and qr (fwd) To: r-help at stat.math.ethz.ch Message-ID: <Pine.SGI.4.40.0306261057160.1716165-100000 at origin.chass.utoronto.ca> Content-Type: TEXT/PLAIN; charset=US-ASCII I have been struggling to find some informaation on what lm exa...

...R -- Dr Peter Dunn (USQ CRICOS No. 00244B) Web: http://www.sci.usq.edu.au/staff/dunn Email: dunn @ usq.edu.au Opinions expressed are mine, not those of USQ. Obviously... ------------------------------ Message: 27 Date: Wed, 10 Nov 2004 19:46:21 -0500 From: Jean Eid <jeaneid at chass.utoronto.ca> Subject: [R] polr probit versus stata oprobit To: r-help at stat.math.ethz.ch Message-ID: <Pine.SGI.4.40.0411101930330.52337307-100000 at origin.chass.utoronto.ca> Content-Type: TEXT/PLAIN; charset=US-ASCII Dear All, I have been struggling to understand why for th...

...such data.frame/header pairs. Is it as simple as using a connection and printing to it everything (including tabs and endlines, etc.), or is there some other way to go about this? Thanks, Sean ------------------------------ Message: 45 Date: Tue, 07 Mar 2006 17:15:27 -0500 From: Jean Eid <jeaneid at chass.utoronto.ca> Subject: Re: [R] reading in only one column from text file To: mark salsburg <mark.salsburg at gmail.com> Cc: R-help at stat.math.ethz.ch Message-ID: <440E05FF.8080403 at chass.utoronto.ca> Content-Type: text/plain; charset=ISO-8859-1; format=flowed You might w...

Dear All, I have been experiencing the following problem when I use a lot of system memory when using R heavily. Everytime I shutdown the process (Xemacs and ESS), my system remains slow and sometime does not respond. I use a debian testing machine. Even when I shutdown (halt) the system, it will crash and on reboot it will check the disk for errors. I have also used gc() after I removed

I have this somewhat large data set that is given to me in a fixed width format. The file itself is already 100MB (Maybe R can actually handle this but I am trying to gain some experience in postgres and RODBC). I am using postgres to preprocess the file and connect to the database through the RODBC package. My question is much before the processing of the database in R (Iknow this is somewhat

This is probably a weird question but I need to know if there is a way... I run an R batch job without saving the variables at each step to the disk. Is there a way to invoke another session of R and link it to the same environment for read only. The problem is that I am running optim with every step getting the parameters into the global env using <<- However, I forgot to issue a

I have a problem either understanding what loess is doing or that loess has a problem itself. As the x-axis variables become more concentrated on a particular point,the estimated loess tends to zero????. the examples below show what i am talking about, why is that? my intution tells me that it should tend to the mean of the variable which is been smoothed. Here's a worked up example x

there does not seem to be a package "tcltk" on CRAN for 2.0.0. I have successfully installed the same package for 1.9.1. In essence I require the package for a GUI interface to setwd. All work fine with 1.9.1. This is on a Linux Debian unstable kernel 2.4.20 The build version of R 2.0.0 (issuing version whithin R) has the following platform i686-pc-linux-gnu arch i686 os