similar to: Individual rename of list items

We have extracted parameters from physiological measurements by fitting SSlogis-like curves with nlsList and nlme. We presented residuals plot in a paper, but a referee argues that these cannot be included (too technical), and r-square values should be given instead to compare the goodness of fit with those of other authors. I remember that 30 years ago in my stat 101, I learned that r-square is

Is there an easy way to remove data frame columns by name instead of by index? The following gives the idea remove<-c("subj","drug") data[-remove] I found a solution with a few evals and substitutes, similar to that used in reshapeLong, but there must be an easier way out. Dieter --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel

Using a formula converted with as.formula with lme leads to an error message. Same works ok with lm, and with lme and a fixed formula. # demonstrates problems with lme and as.formula demo<-data.frame(x=1:20,y=(1:20)+rnorm(20),subj=as.factor(rep(1:2,10))) demo.lm1<-lme(y~x,data=demo,random=~1|subj) print(summary(demo.lm1)) newframe<-data.frame(x=1:5,subj=rep(1,5))

Dear R-Listers, to measure the psychometric curve of pitch discrimination, one sequentially presents two tones of slightly different pitch to an observer (animal will do), and asks "which is higher". The pschometric curve is the fraction of correct responses plotted against the pitch difference. It passes through 50% (pure guessing) at zero and normally approaches 100% at large

Dear friends of nlme, Running quinModel (Pinheiro/Bates page 380) on R (current release, windows) gives: Nonlinear mixed-effects model fit by maximum likelihood Model: conc ~ quinModel(Subject, time, conc, dose, interval, lV, lKa, lCl) Data: Quinidine Log-likelihood: -497 Fixed: lV + lKa + lCl ~ 1 lV lKa lCl 5.382 -0.273 2.470 Random effects: Formula: list(lV ~ 1, lCl ~ 1)

deriv seems to have problems with a minus-sign before a bracket. Below are four examples of the same function, the top one is wrong, all others are correct (hopefully). Rest of expression not shown, it is the same for all versions. _ platform i386-pc-mingw32 arch x86 os Win32 system x86, Win32 status major 1 minor 3.0 year 2001 month 06 day 22 language R

I have a two-dimensional loess fit, and need to calculate the gradient field from it. Even after looking at loess.c and loess.f, I don't understand the meaning of the returned polynomial coefficients. Or is the brute force method of using a tangential approx to the fitted values the way to go? Dieter Menne --------------------------------------- Dr. Dieter Menne Biomed Software 72074

I have problem getting the concept of a nested fixed variable into the nlme scheme. I fear the question is very stupid. In the past I had asked this before, and never got a reply (in other cases, the response was within hours). I also checked the S-list, where several similar enquiries of other people are orphaned. We have a cross-over design, where patient are treated two weeks with placebo,

Dear Mixed Effect Friends, Somehow, R(1021, Windows) seem to run differently from S Plus: The soybean example from Pinheiro/Bates on page 290 fails in R. (Soybean1 is Soybean with the NA and "critical" case removed. Same procedure with full Soybean). > fm1Soy.lis<-nlsList(weight~SSlogis(Time,Asym,xmid,scal),data=Soybean1) > fm1Soy.nlme<-nlme(fm1Soy.lis) Error: Singularity

Dear R-supporters, I have results of lme similar to those shown in Fig. 1.21, p.51 of Pinheiro/Bates. However, In my data set, one of the panels is missing, leading to an ugly frame shift of the following panels. How can I replace one of the panels by an empty one to restore the raster? Dieter Menne --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel

I try to display the grouping variable in coplot. It work, but it's special solution and rather ugly. Any better idea? # Simulate my data frame data(state) x77<-data.frame(state.x77) x77$region<-state.region coplot(Life.Exp ~ Income | region, data=x77, show.given=F, subscripts=T, panel = function(x, y,subscripts, ...) { panel.smooth(x, y, span = 1., ...)

Dear R's, I have a survey where customers rank a set of 5 packages for a product, so the response variable looks like a d b c a c d b d b a c Predictors variables are 4 socio-economic parameters. I have modelled the FIRST choice of each subject as a multinomial model, similar to the housing example in MASS ch7.3, , but I would prefer to use the whole rank-set instead. Can someone give me a

Hi Dieter, Sorry my mistake. I wanted to convert them into Decimal (not Hexadecimal). Given this string, the desired answer follows: > ascii_str <- "ORQ>IK" 79 82 81 62 73 75 > ascii_str2 <- "FDC" 70 68 67 - Gundala Viswanath Jakarta - Indonesia On Mon, Dec 22, 2008 at 5:49 PM, Dieter Menne <dieter.menne at menne-biomed.de> wrote: > Gundala

Thanks to Dieter Menne and Spencer Graves I started to get my way through lsoda() Now I need to use it in with nls() to assess parameters I have a go with a basic example dy/dt = K1*conc I try to assess the value of K1 from a simulated data set with a K1 close to 2. Here is (I think) the best code that I've done so far even though it crashes when I call nls()

Hi Dieter and R community: I tried both of these three versions with ylim as suggested, none work: I am getting only single (pch = 16) not overlayed (pch =3) everytime. *vs 1* require(lattice) xyplot(Sepal.Length ~ Sepal.Width | Species , data= iris, panel= function(x, y, subscripts) { panel.xyplot(x, y, pch=16, col = "green4", ylim = c(0, 10)) panel.lmline(x, y, lty=4, col =

I have a dialplan that looks like this: [dorecord] exten => _*99XX,1,Verbose(2,Doing custom record) exten => _*99XX,n,Answer() exten => _*99XX,n,Verbose(2,Doing custom record - before wait) exten => _*99XX,n,Wait(0.5) exten => _*99XX,n,Verbose(2,Doing custom record - before record) exten => _*99XX,n,Record(/tmp/prompt${EXTEN:3}.gsm) exten => _*99XX,n,Verbose(2,Doing custom

I created a report with Sweave today, that displayed perfectly on screen, but crashed both Adobe Pro 9.3.1 and Adobe Reader 9.0 on Windows 7. Output with Foxit Reader was flawless. I was able to reproduce a "minimal" example, which is not really minimal but the smallest I could get after 2 hours of wasting paper. Required: latticeExtra calling both panel.xyplot(...) # both

Hi folks, I've noticed that there's an issue with the recursive creation of directories that reside on network shares. For example: > dir.create('\\\\SERVERNAME\\Empl\\Home1\\active\\e\\ecortens\\thisisatest', recursive = TRUE) Warning message: In dir.create("\\\\SERVERNAME\\Empl\\Home1\\active\\e\\ecortens\\thisisatest", : cannot create dir

Dear UseR, I do not know if this a problem with me, my data or cph/survest in package design. The example below works with a standard data set, but not with my data, but I cannot locate the problem. Note that I am using an older package of survival to avoid a problem with the newly renamed function in survival meeting Design. Dieter # First, check standard example to make sure library(Design)

Hello everyone, I need some help with lists inside lists (a good way to emulate a struct)\ Assume that there is a small list called fred: fred <- list(happy = 1:10, name = "squash") and a big list called bigfred that includes fred list 5 times bigfred <- rep(fred,5) Is it possible somehow to index all these sublists(fred) inside bigfred with a more direct way like this: