similar to: Residual DF of NLS models (PR#14194)

Colleagues, I am using "nls" successfully (2.11.1, OS X) but I am having difficulties retrieving part of the output - residual sum of squares. I have assigned the output to FIT: > > FIT > Nonlinear regression model > model: NEWY ~ PMESOR + PAMPLITUDE * cos(2 * pi * (NEWX - POFFSET)/PERIOD) > data: parent.frame() > PMESOR PAMPLITUDE POFFSET >

move or copy the directories mass, nnet and class to the library directory - then execute link.html.help() Now execute library(MASS) and data(petrol) Should work. "Troels Ring"

Dear friends. I have trobpe with path. Below I have tried to paste a path and file name but it will not work > path <- "D:/rw0641/own/own/procedures/GFR_TIME" > path [1] "D:/rw0641/own/own/procedures/GFR_TIME" > path1 <- paste(path,"gfr.txt",sep="/") > path1 [1] "D:/rw0641/own/own/procedures/GFR_TIME/gfr.txt" -- this is good

Hi All, I was wondering whether there have been some listening experiments done to test how well spatial information is preserved in the celt signal, e.g. comparison of sound localization performance for the original uncompressed sound and the celt sound (most probably for different bit rates). Best, Pablo -- Pablo F. Hoffmann PostDoc Acoustics Dept. of Electronic Systems Aalborg

Hi, I'm using windows xp and R 2.8.0 - I wonder what is the command to put in a script that has the same effect as when in a plot you choose menu "History" and "Recording". Best wishes Troels -- Troels Ring - - Department of nephrology - - Aalborg Hospital 9100 Aalborg, Denmark - - +45 99326629 - - tring at gvdnet.dk

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi In the Samba HOWTO (ch. 24 section 24.3.2) is written "When winbind is used, the default condition is that the local user george will be treated as the account DOMAIN\george and the foreign (non-member of the domain) account will be treated as MACHINE\george because each has a different SID." I have enable 'winbind use default

I was attempting to use the residual resampling approach to generate 999 bootstrap samples of alpha and beta and find their confidence intervals. However, I keep getting the error message:Error in nls(resample.mp ~ cases/(alpha + (beta * cases)), start = init.values, : singular gradientafter R has only produced a few bootstraps.Could anyone suggest where I am going wrong? Would greatly

Dear friends - I'm on winXP, R 2.71 - I have with some help dveloped this multivariate normal model, which gives very plausible results in WinBUGS even without any initial values specified. However, when I then try to run the same model via the bugs function in R2WinBUGS with inits specified as inits=NULL the program stops in a dead end. So I have tried to make inits for the bugs function

Dear Friends, I'm using winXP and R 1.7.1 and plotting some data using dates on the x-axis, and wanted to use identify to show some points but was told by identify that the x and y vectors producing a fine graph with 84 points were not equal in length. Below are the Dato for date - and length(Dato) finds 9 but str finds 84 as known. Will identify not work in this context ? Best wishes

We would like to fit parameters using a simulation with stochastic processes as theoretical values. We generate a simple exemple with nls.lm to see the logic and the problem: First without stochasticity (it is a dummy example, the fited value is simple the mean of a set of 10 numbers): #Ten numbers x <- 1:10 #Generate 10 Gaussian random number with mean=3 sd=1 simy <- rnorm(length(x),

Dear friends, I was drawing a circle with centrum in (1,-1) and radius 5 to show my girl that the line y=3*x+1 goes through (1,4) and (-2,-5) of the circle, but on Windows XP, R 2.20 the drawing was not good at all, and the known solutions were not shown in the graph. I guess I got it wrong? Is this use not intended ? Best wishes Troels Ring, MD Aalborg, Denmark

Dear friends - I'm using Tinn-R in Windows and find it OK. My son is running MacOS and is just starting using R and I know nothing about MacOS and R but we seem to have difficulties in finding something similar to Tinn for MaC. What are your suggestions? Thanks a lot! Troels Ring, MD Aalborg, Denmark

Hi, I was trying to run the example of Indomethacin kinetics from the book: ## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus, ## Springer, Second Printing 2001, Section 6.2 library(nlme) plot(Indometh) fm1Indom.nls <- nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh) summary(fm1Indom.nls) plot(fm1Indom.nls,Subject~resid(.),abline=0) ## .... the last plot command gives me the

Dear friends. I just downloaded tseries dated 27.10.99 (from Aalborg, Denmark) and it seems it lacked something, library ts was not available ? What to do ? > library(tseries) Loading required package: ts Warning: There is no package called `ts' > > Best wishes Troels Ring, M.D Department of Nephrology Aalborg, Denmark tring at mail1.stofanet.dk

Hi list, I tried to use nls to do some nonlinear least square fitting on my data with 340 observations and 10 variables, but as I called nls() function, I got this error message: Error in qr.qty(QR, resid) : 'qr' and 'y' must have the same number of rows Then I traced back a little bit into nls() function, the error seemed to happen when calling nlsiter() internal function, but

I hope you will forgive me this simple question on titration. I'm trying to find very small values from the algorithm below, which I believe is correctly formatted, and the constants are also correct. When SID goes over ATOT, fitted vales are much too low compared to the literature. I guess I must be using optimize in a wrong way but cannot find out how to improve it. I'm on windows, R

Dear R helpers I am trying to convert a list into a data frame but when I try, I get a stack overflow error (Error: protect(): stack overflow). My list contains about 17000 rows and looks like shown at the bottom. The reason that I want to convert it in to a data frame is that I want to export it to a mysql database with the dbWriteTable function. The function that I use is

Dear list I would like to be able to group terms in a formula using a function that I will call tvar(), eg. the formula Y ~ 1 + tvar(x:A) + tvar(z) + u + tvar(B) + tvar(poly(v,3)) where x,u and v are numeric and A and B are factors - binary, say. As output, I want the model.matrix as if tvar had not been there at all. In addition, I would like to have information on the grouping, as a vector

Dear R-experts! I have y(x) data, dim(y) == dim(x) == c(2000, 2) I'd like to fit them with nls: fit.result <- nls ( y ~ f(x, p1, p2, p3), start = list(p1 = ... , p2 = .. , p3 = ..) ) Actually I want to fit y[,1] ~ x[,1] and y[,2] ~ x[,2] *simulaneously*, with the same parameters set {p1, p2, p3}. I tried to feed R tha above formula, R errors with:

Dear friends - I hope you will accept a naive question on lm: R version 3.4.1, Windows 10 I have 204 "baskets" of three types corresponding to factor F, each of size from 2 to 33 containing measurements, and need to know if the standard deviation on the measurements in each basket,sdd, is different across types, F. Plotting the observed sdd versus the sizes from 2 to 33, called