search for: asymp

Hi folks, I am tearing my hair out on this one. I am using an example from Pinheiro and Bates. ### this works data(Orange) mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange ) ### This works mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange, fixed = Asymp + xmid + scal ~ 1, start = fixef(mod.lis) ) ### I try a slightly different model specification for the fixed...

...to write a small package that based on nlme however my code does not work. R always appears this message: Error in eval(expr, envir, enclos) : object "y" not found where y is the response variable. Please help me out! This is my code: require(nlme) AMPmixed<-function(y, x, S1=c("asymptotic","logistic"), random,data, S2=c("asymptotic","logistic"), start) { logist.logist<-function(x,alpha,delta,psi.l,tau.l,gamma,h){ delta+(alpha-delta+gamma*(x-(h-1))/exp(x))/(1+exp(-(x-tau.l)/psi.l))} logist.asymp<-...

...n, I have encountered the similar problem Dr. Stevens and Dr. Graves observed (please see the posts: https://stat.ethz.ch/pipermail/r-help/2006-May/105832.html ). I have tried Dr. Stevens's original example, the problem is still there, > mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), + data=Orange ) > > mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal), + data=Orange, + fixed = Asymp + xmid + scal ~ 1, + start = fixef(mod.lis) ) > > mod <- nlme(circumference ~ SSlogis(age, Asymp,...

Hi, I have a program written in R which is good on the version 1.2, but for the fallowing versions of R, an error always is at the same place. That is at the level of the fallowing line: Sur<- getInitial(res2[m:M,2]~SSasymp(res2[m:M,1],Asymp,resp0,lrc),data=res2) Error in eval(expr,envir,enclos):numeric envir arg not of length one I don't know at all this langage for the instant. Thanks, Christelle ------------------------------------------------------------- NetCourrier, votre bureau virtuel sur Internet...

Hi list, I tried to load the package SAGx and failed because it complains it's looking for the Biobass which is not there. Then I looked up the package list and Biobass is not found. I'm trying to run the Jonckheere-Terpstra test and from what I see in the R archive, SAGx is the only place it's been implemented. > library(SAGx) Loading required package: multtest Loading required

...the 2 variables. The overall chisquare value is much lower than the value obtained using R, though still significant. pre02 & post02 pre02 post02 0 1 0 311 20 1 119 9 Test Statistics(b) pre02 & post02 N 459 Chi-Square(a) 69.094 Asymp. Sig. .000 a Continuity Corrected b McNemar Test Any assistance is much appreciated, Bob Green

...especially the part that does the exact distribution. I'm looking forward to comments on improving the programming efficiency for this problem. (The "return(...)" seems to be necessary in R only.) Knut #-------------------------------------------------------------------------- # asymp.SMN.pvalue(pP,qP,pX,qX,pQ,qQ) # exact.SMN.pvalue(pP,qP,pX,qX,pQ,qQ) #-------------------------------------------------------------------------- # pP,qP = number of PP,PQ children of PP~PQ parents # pX,qX = number of PP,QQ children of PQ~PQ parents # pQ,qQ = number of PQ,QQ children of PQ~QQ parents...

An earlier post had posed the question: "Does anybody know what is relation between 'T' value calculated by 'wilcox_test' function (coin package) and more common 'W' value?" I found the question interesting and ran the commands in R and SPSS. The W reported by R did not seem to correspond to either Mann-Whitney U, Wilcoxon W or the Z which I have more