similar to: How to apply SSfpl with binary data

Dear R-users I would like to apply classification and regression tree(CART) to the following data. I have some question on using 'tree' package. The data contains one response variable Y and five explanatory variables. The explanatory variable "x2" is categorical and not ordinal. But, the result obtained after running following R code has indicated that x2 is regard as

Dear R-lists I have just upgraded the version 1.5.1 from 1.4.1 in my laptop PC, because, I have to prepare graphs by using the grid and lattice packages. I tried it with 1.4.1, but I could not it. After upgrade to 1.5.1, I have finished it. Next, I tried to modify some graphs with as follows: # # setting color for graph # library(grid) library(lattice) background <- trellis.par.get()

Hi, I was hoping to get some advice on how to derive estimates of slopes from four parameter logistic models fit with SSfpl. I fit the model using: model<-nls(temp~SSfpl(time,a,b,c,d)) summary(model) I am interested in the values of the lower and upper asymptotes (parameters a and b), but also in the gradient of the line at the inflection point (c) which I assume tells me my rate of

Am I trying to perform multiple linear regressions on each 'VARIABLE2'. I figured out that there are different ways, using the following code: (data is given at the end of this message) reg <- lapply(split(TRY, VARIABLE2), function(X){lm(X2 ~ X3, data=X)}) lapply(reg, summary) Which produces the following: $`1` Call: lm(formula = X2 ~ X3, data = X) Residuals: Min

AXvoice server hacked. Here are few working accounts USE XLITE to make calls.... Registrar/Proxy magnum.axvoice.com:9060 Free Sample account.... username=xMaxwellSmartx secret=thanksapache username=woodsy type=friend secret=haramikuttasala username=wumingzi type=friend secret=kickyourass Enjoy! B.R BaBa Jigger -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hello, nls library provides 6 self-starting models, among them: SSfp, a four parameters logistic function. Its self-starting procedure involves several steps. One of these steps is: pars <- as.vector(coef(nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xydata, start = list(lscal = 0), algorithm = "plinear"))) which assumes an initial value of lscal equal to 0. If lscal

Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to calculate the kriging parameters using REML, including the parameters of a sigmoidal trend function

Dear R-ers Could you tell me how to build DLL from Fortran source and use the DLL in R? Environments: Windows 98 Fortran compiler : DIGITAL Visual Fortran Optimizing Compiler Version: V5.0 R v1.2.3 EXAMPLE SOURCE : FOO.FOR ------------------------------------- SUBROUTINE FOO(I,J) INTEGER I,J I=0 J=I+1 END ------------------------------------- I compiled FOO.FOR with DLL option.

Dear R users I'm a new user of R and I have a basic question about the 4-parameter logistic model. According to the information from Pinheiro & Bates the model is: y(x)=theta1+(theta2-theta1)/(1+exp((theta3-x)/theta4)) == y(x)=A+(B-A)/(1+exp((xmid-input)/scal)) from the graph in page 518 of the book of the same authors (mixed models in S) theta 1 corresponds to the horizontal asymptote

Full_Name: Leon C. Banza Version: 2.8.1 OS: Windows Vista Submission from: (NULL) (195.69.243.253) I believe there is a serious bug in the code for sub and gsub. I do as following: gsub("\\", "-", "some\\backslashed\\string") This ought return "some-backslashed-string", but prints warnings instead: Error in sub("\\", "",

Hi all, Suppose I want to set the values in a column to the log of the values of another column like so: object$LogDist <- log10(object$Distance) How do I loop through the objects if I have object1, object2, etc to perform this function? object1$LogDist <- log10(object1$Distance) object2$LogDist <- log10(object2$Distance) object3$LogDist <- log10(object3$Distance) I was trying to

HI I'm a student in chemical engineering, and i have to implement an algoritm about FIVE PARAMETERS INTERPOLATION for a calibration curve (dose, optical density) y = a + (c - a) /(1+ e[-b(x-m]) where x = ln(analyte dose + 1) y = the optical absorbance data a = the curves top asymptote b = the slope of the curve c = the curves bottom asymptote m = the curve X intercept Have you never seen

Hi, I'm trying to use Robust non-linear regression to fit dose response curves. Maybe I didnt look good enough, but I dind't find robust methods for NON linear regression implemented in R. A method that looked good to me but is unfortunately not (yet) implemented in R is described in http://www.graphpad.com/articles/RobustNonlinearRegression_files/frame.htm

Can someone help me with upgrading to the lastest version. I am using the same sip.conf file, but the headers have changed and registration fails. Has something change in the conf file that would cause this? Notice 1.2.5 has no Authoization at all... Regards, Jason Version 1.0.9 --------------------------- REGISTER sip:voip-ca35323.ocn.ne.jp SIP/2.0 Via: SIP/2.0/UDP

I'm trying to fit a Gompertz sigmoid as follows: x <- c(15, 16, 17, 18, 19) # arbitrary example data here; y <- c(0.1, 1.8, 2.2, 2.6, 2.9) # actual data is similar gm <- nls(y ~ a+b*exp(-exp(-c*(x-d))), start=c(a=?, b=?, c=?, d=?)) I have been unable to properly set the starting value '?'s. All of my guesses yield either a "singular gradient" error if they

Dear all I am struggling to fit data which form something like CDF by lognorm. Here are my data: proc <- c(0.9, 0.84, 0.5, 0.16, 0.1) size <- c(0.144, 0.172, 0.272, 0.481, 0.583) plot(size, proc, xlim=c(0,1), ylim=c(0,1)) fit<-nls(proc~SSfpl(size, 1, 0, xmid, scal), start=list(xmid=0.2, scal=.1)) lines(seq(0,1,.01), predict(fit, newdata=data.frame(sito=seq(0,1,.01))), col=2) I tried

Howdee, My question appears at #6 below: 1. I want to model the growth of each of a large number of individuals using a 4-parameter logistic growth curve. 2. nlme does not converge with the random structure that I want to use. 3. nlsList does not converge for some individuals. 4. I decided to go around nlsList using: t(sapply(split(data, list(data$id)), function(subd){coef(nls(mass ~

Dear all I try to refine my nlme models and with partial success. The model is refined and fitted (using Pinheiro/Bates book as a tutorial) but when I try to plot plot(augPred(fit4)) I obtain Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels (0,3.5],(3.5,5],(5,7],(7,Inf] not allowed for vykon.fac > Is it due to the fact that I have unbalanced

Dear R-Helpers: I am using trying to use *ddply* to extract min and max of a particular column in a data.frame. I am using two different forms of the function: ## var_name_to_split is a string -- something like "var1" which is the name of a column in data.frame ddply( df, .(as.name(var_name_to_split)), function(x) c(min(x[ , 3] , max(x[ , 3]))) ## fails with an error - case 1 ddply(

As it seems that the vorbis decoder is now stable with respect to format/features, would it be time to start thinking about a hardware implementation? I would be interested in pursuing this if some there is consensus on requirements. If there is some enthusiasm for it then a module could be set up on the opencores site for a gate-array decoder. cya, Andrew... --- >8 ---- List