similar to: Where is the bioDist package?

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

Hi, I am observing a deadlock with llvm-lit on windows 7. When I attached a debugger, the communicate() call is blocked. In file utils/lit/lit/TestRunner.py > # FIXME: There is probably still deadlock potential here. Yawn. > procData = [None] * len(procs) > procData[-1] = procs[-1].communicate() I am invoking python directly on windows to run the unit tests.

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

I am confused by the logic of renaming: # Rename individual list items? Empl<-list(employee="Anna",spouse="Fred") names(Empl)<-c("empl","spo") names(Empl) #[1] "empl" "spo" # worked like a charm... but names(Empl[1])<-"newempl" # no error message, yet .... names(Empl) #[1] "empl" "spo" #

I can't get the PERL subexpression translated to R. Following, for example, B. Ripley's http://finzi.psych.upenn.edu/R/Rhelp02a/archive/58984.html I am using sub, but it looks like an ugly substitute. Assume I want to extract the first alpha part and the first numeric part, but only if they are in sequence. Do I really have to use the sub twice, first extracting the first variable, then

Hi there, I hope I can get advice regarding the calculation of leverage values or studentized residual values of a non-linear regression model. It seems like rstudent() does not work on a nls object. Many thanks in advance! Best regards, Xingli

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)

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 =

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

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

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

Dear List, creating factors in a given non-default orders is notoriously difficult to explain in a course. Students love the ifelse construct given below most, but I remember some comment from Martin M?chler (?) that ifelse should be banned from courses. Any better idea? Not necessarily short, easy to remember is important. Dieter data = c(1,7,10,50,70) levs =

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 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 R and lme4 users- I am trying to fit a mixed-effects model, with the glmer function in lme4, to right-skewed, zero-inflated, non-normal data representing understory grass and forb biomass (continuous) as a function of tree density (indicated by leaf-area). Thus, I have tried to specify a Gamma distribution with a log-link function but consistently receive an error as follows: >

I read somewhere that vector graphics such as eps or dpf are more favorable than alternatives (jpeg, bmp or png) for publication because vector graphics scale properly when enlarged. However, my problem is that the file generated from a graph of fixed size is too large (in the order of 10MB) because of many data points in multiple scatterplots. Any suggestions? Thanks, Gang [[alternative HTML

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

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

Dear List, After having installed R 2.6.0/Windows 2000, my trusted command script to build and install my package fails with an error message I do not understand. I have installed a fresh copy of Rtools, and made sure that these are the first on the path. For testing, I am using a package version from an older backup where I am sure that everything else was ok.