similar to: nlsrob fails with puzzling error message on input accepted by nls

I have a question regarding robust nonlinear regression with nlrob. I would like to place lower bounds on the parameters, but when I call nlrob with limits it returns the following error: "Error in psi(resid/Scale, ...) : unused argument(s) (lower = list(Asym = 1, mid = 1, scal = 1))" After consulting the documentation I noticed that upper and lower are not listed as parameter in

Hi all, I'm trying to fit a nonlinear regression by "nlrob": model3=nlrob(y~a1*x^a2,data=transient,psi=psi.bisquare, start=list(a1=0.02,a2=0.7),maxit=1000) However an error message keeps popping up saying that the function psi.bisquare doesn't exist. I also tried psi.huber, which is supposed to be the default for nlrob: model3=nlrob(y~a1*x^a2,data=transient,psi=psi.huber,

Dear all, I am not sure if this mail is for R-help or should be sent to R-devel instead, and therefore post to both. While using nlrob from package 'robustbase', I ran into the following problem: For psi-functions that can become zero (e.g. psi.bisquare), weights in the internal call to nls can become zero. Example: d <- data.frame(x=1:5,y=c(2,3,5,10,9)) d.nlrob <-

Dear all, ? I implemented a non-linear model using nls(...) and it works just fine. I now tried to run the same model using nlrob(...) which basically does the same but uses a more robust estimation procedure. My problem: I cannot seem to get nlrob(...) running. Irrespective of how I try to call the function, I always get the error message "Error in is.null(x) : 'x' is missing".

What am I failing to understand here? The script below works fine if the dataset being used is DNase1 <- DNase[ DNase$Run == 1, ] per the example given in help(nlrob). Obviously, I am trying to understand how to use nls and nlrob to fit curves to data using R. #package=DAAG attach(codling) plot(pobs~dose) #next command returns 'step factor reduced below min factor

R-help, Sorry for posting the same question (almost) once again. I'm using a lm model for some data. The code is below: fitData <- glm(pm ~ age + mw , family = binomial(link = logit), maxit = 30, data = inputModel1) "age" is a factor and "mw" is a continous variable. The summary function does not provide the parameter estimates for the first level of the

Hello, Are there any R functions available for performing a serial correlation test for short time series (e.g, series having between 10-14 observations)? Many thanks! Isabella R. Ghement, Ph.D. Ghement Statistical Consulting Company 301-7031 Blundell Road, Richmond, B.C., Canada, V6Y 1J5 Tel: 604-767-1250 Fax: 604-270-3922 E-mail: isabella at ghement.ca Web: www.ghement.ca -----Original

I'm working with 3 different data sets and applying this non-linear regression formula to each of them. nls(Y ~ (upper)/(1+10^(X-LOGEC50)), data=std_no_outliers, start=list(upper=max(std_no_outliers$Y),LOGEC50=-8.5)) Previously, all of the regressions were calculated in Prism, but I'd like to be able to automate the calculation process in a script, which is why I'm trying to move to

Hi all, I made sure that it's "env$sRes1$nPositionsOptimizedM" that's correct... not the "env$sRes1$nPositionsOptimized"... But it seems both point to the same memory area... This is very dangerous because I have used naming conventions such as: MyLongVariableNameForA MyLongVariableNameForB MyLongVariableNameForC ... ... Then if internally they are actually the

I'm testing a new release of survival, executing the following piece of code: .... for (testpkg in survdep) { z <- testInstalledPackage(testpkg, outDir="tests") cat(testpkg, c("Ok", "Failed")[z+1], "\n", file="progress", append=T) } The vector survdep contains the names of all 156 packages listed as reverse depends on the CRAN

hi everyone a good basher around here? I try in a script: _rdiffBack="rdiff-backup -v5 --tempdir /tmp/ --no-eas --exclude-other-filesystems --exclude-symbolic-links" _rdffiExclude="--exclude '**/~*' --exclude '**.tmp' --exclude-regexp '(.glusterfs|.trashcan|temp)'" _execCom=${!2} _sourceDir=${1} _backupTo=${3} __backMeUp() { for _sourceDir in

Hi, I cannot make a reproducible example easily for my problem, so I'll describe it as best as I can. I merged 2 dataframes but was surprised when one line on the x dataframe did not get a match in the y dataframe, because I knew such a match existed. There was only one "by" variable in the merge, in Date format: in x: $ p?riode : Date, format: "2009-06-09"

I wrote an R program that does heavy computations with hundreds of lines of code. It's running fine both interactively and in batch mode on a Mac OS X computer. The program also has no problem running on a Linux system (Fedora 14) interactively. However, when I try it on the terminal in batch mode on the Linux system, it chokes in the middle of the computation with the "Execution

Hi, I just completed a really big rsync described earlier. Ie about 13,945 directories transfered about 600GB of data. Of 13,945 directories, 13,9441 directories transfer with matching du -b sizes of the preimage to the size of the destination machine image. of the 4 remaining directories i found source vs destination in bytes --------- a) 20480 vs 34922496 b) 28672 vs

Hi everyone, I am running a monte carlo and am getting an error that I haven't the slightest clue where to begin figuring it out. The error is as follows: Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases In addition: Warning message: In ltm.fit(X, betas, constraint, formula, con) : Hessian matrix at convergence is not positive definite;

Hi all, As you can see from below, the result is strange... I would imagined that the bb result should be much higher and close to 1, any way to improve the fit? Any other classification methods? Thank you! data=data.frame(y=rep(c(0, 1), times=100), x=1:200) aa=glm(y~x, data=data, family=binomial(link="logit")) newdata=data.frame(x=6, y=100) bb=predict(aa, newdata=newdata,

Hi, It is not a Samba issue, but a RH9 one. Try to replace in /etc/sysconfig/i18n the default charset utf8 to something easier (and faster) to handle for the system (look at the same file on your RH73). HTH, J?r?me -------- Message d'origine-------- De: Stefan G. Weichinger [mailto:monitor@oops.co.at] Date: mar. 9/30/2003 20:23 ?: samba@lists.samba.org Cc: Objet: Re:

Hi RedHat 7.2, ext3 on /, kernel 2.4.18p1. whilst updatedb was running, i had these messages appear... Jan 6 22:18:42 jaguar kernel: EXT3-fs error (device ide0(3,3)): ext3_readdir: bad entry in directory #147553: rec_len %% 4 != 0 - offset=0, inode=1651076143, rec_len=19527, name_len=85 Jan 6 22:18:42 jaguar kernel: EXT3-fs error (device ide0(3,3)): ext3_readdir: bad entry in directory

Greetings folks: Stepped away from a win/lin dual boot system and spent the better part of the past week resetting a windows only arrangement. Before the switch I was running 2.3.1 without any strange behaviours. When I reached the point of putting it back on the machine I would get this puzzling warning the first time I would try to select a cran mirror: > chooseCRANmirror() Warning

Dear all, I'm trying to connect to an MSAccess database (ArcGIS personal geodatabase). I keep getting an error about the channel when using sqlQuery(). However, sqlTables() does not complain about the channel and lists all tables in the database. If I try sqlFetch(), then R crashes. I'm happy to hear suggestions on how to solve this. Best regards, Thierry > MDB <-