search for: menne

...but I cannot find a referee-proof citation for this. Probably SAS does it (I know what Bill Venables will say about this argument...). Can anybody provide me a citation on r-square for/against non-linear fits, or tell me that I am wrong? Dieter --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send &quo...

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

...stead 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 (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "i...

...es(Empl[1])<-"newempl" # no error message, yet .... names(Empl) #[1] "empl" "spo" # where's my newempl??????? If this line names(Empl[1])<-"newempl" was stupid, why is there no error message? --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "i...

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

...vide the counts for a stimulus difference equally between positive and negative branch, with correct/wrong inverted for the negatives. The solution is not clean if we have uneven counts. Is this a valid approach? Alternatively, could I force a glm/binomial/logistic fit to go through 0/0.5? Dieter Menne --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read...

...assume I want to make a minor local change in a package that is supplied as .rda. For example, I want to get rid of the non-verbose-protected "Iteration" message in GLMM/lme4. Probably I have to load / change / save the package, but could someone help me to get the syntax right? Dieter Menne

...level=0) # ok demo.lm2<-lme(fma,data=demo,random=~1|subj) # looks ok, but isn't print(summary(demo.lm2)) #Fixed effects: Error in as.vector(x, "list") : cannot coerce to vector #predict(demo.lm2,newframe,level=0) # does not work --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "i...

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 =

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

...on: is the source of nlme_one_comp_open (or similar) available? I would like to learn how to do these models in c, but could not find the code. (To avoid RTFM: I know about r-ext.pdf, but a special modelxXX example would make life easier). Dieter --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-..-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "...

...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 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read h...

...with print(...,split,more=T) or (...position). Sometimes I wish I had something like the old "par(mfcol=2, mfrow=3)" mechanism, where the next free viewport is automatically chosen. I tried fiddling with grid-viewports, but could not find an easy solution. Did I miss something? Dieter Menne

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

...<.0001 #TreatPla 5.106889 1.096139 18 4.658980 0.0002 #Seq2 0.466712 2.265797 18 0.205981 0.8391 #TreatPla:Seq2 3.525737 1.550174 18 2.274413 0.0354 # Sequence as random ??????? glu.lme<-lme(Glu~Treat,data=glu,random=~1|Seq/Patient) print(summary(glu.lme)) Dieter Menne --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read...

...(~offv+Asym/(1+exp((xmid-x)/scal)),"x") expression({ .expr3 <- exp((xmid - x)/scal) ## Coorrect > deriv(~offv+Asym/(1+exp(-1.0*(x-xmid)/scal)),"x") expression({ .expr5 <- exp(-1 * (x - xmid)/scal) ## Correct --------------------------------------- Dr. Dieter Menne Biomed Software 72074 Tübingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne@menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "in...

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 s...

...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 Viswanath <gundalav <at> gmail.com> writes: > >> For example I have the following ASCII strings: >> >> > ascii_str <- "ORRRIROOQRQQOLORRRRRRQRROKK>IKFAA" >> > ascii_str2 <- &...

...green4") panel.xyplot (iris$Petal.Length [subscripts], iris$Petal.Width[subscripts], pch= 3, col = "red", ylim= c(0,10) ) panel.lmline(iris$Petal.Length [subscripts], iris$Petal.Width [subscripts], col = "red") }, as.table=T, subscripts=T) Help please: From: Dieter Menne <dieter.menne@menne-biomed.de> > Date: Wed, Jul 27, 2011 at 8:44 AM > Subject: Re: [R] lattice overlay > To: r-help@r-project.org > > > > > Ram H. Sharma wrote: > > > > I want to overlay lattice scatter plot: I do not know why the following > > code &...

...lution below let's me wait (there are about 10* cases in the real data set). Any suggestions to speed this up? Yes, I could wait for the result in real life, but just curious if I did something wrong. In real life, data set is ordered by groups, but how can I use this with a data frame? Dieter Menne # Generate test data ncol = 20 nrow = 20000 ngroup=nrow %/% 20 colrow=ncol*nrow group = factor(floor(runif(nrow)*ngroup)) sc = data.frame(group,matrix(ifelse(runif(colrow) > 0.1,runif(colrow)>0.3,NA), nrow=nrow)) # aggregate system.time ({ s = aggregate(sc[2:(ncol+1)],list(group = gr...