similar to: "plinear"

Displaying 20 results from an estimated 1000 matches similar to: ""plinear""

2008 Jul 09
2
Port package
Hi When I type: > ?nls I come across this section: algorithm: character string specifying the algorithm to use. The default algorithm is a Gauss-Newton algorithm. Other possible values are '"plinear"' for the Golub-Pereyra algorithm for partially linear least-squares models and '"port"' for the 'nl2sol'
2005 Oct 31
2
nls() fit to Kahnemann/ Tversky function
Dear WizaRds, I would like to fit a curve to ten points with nls() for one unknown parameter gamma in the Kahnemann/ Tversky function, but somehow it won't work and I am unable to locate my mistake. p.kum <- seq(0.1,1, by=0.1) felt.prob.kum <- c(0.16, 0.23, 0.36, 0.49, 0.61, 0.71, 0.85, 0.89, 0.95, 1) ## how to find a function that fits these points nicely? plot(p.kum,
2018 May 05
0
Bug in profile.nls with algorithm = "plinear"
Dear sirs It seems like there is a bug in `profile.nls` with `algorithm = "plinear"` when a matrix is supplied on the right hand side. Here is the bug and a potential fix ##### # example where profile.nls does not work with `plinear` but does with # `default` require(graphics) set.seed(1) DNase1 <- subset(DNase, Run == 1) x <- rnorm(nrow(DNase1)) f1 <- nls(density ~ b1/(1 +
2008 Oct 02
1
nls with plinear and function on RHS
Dear R gurus, As part of finding initial values for a much more complicated fit I want to fit a function of the form y ~ a + bx + cx^d to fairly "noisy" data and have hit some problems. To demonstrate the specific R-related problem, here is an idealised data set, smaller and better fitting than reality: # idealised data set aDF <- data.frame( x= c(1.80, 9.27, 6.48, 2.61, 9.86,
2008 May 06
2
NLS plinear question
Hi All. I've run into a problem with the plinear algorithm in nls that is confusing me. Assume the following reaction time data over 15 trials for a single unit. Trials are coded from 0-14 so that the intercept represents reaction time in the first trial. trl RT 0 1132.0 1 630.5 2 1371.5 3 704.0 4 488.5 5 575.5 6 613.0 7 824.5 8 509.0 9
2005 Jun 21
2
nls(): Levenberg-Marquardt, Gauss-Newton, plinear - PI curve fitting
Hello, i have a problem with the function nls(). This are my data in "k": V1 V2 [1,] 0 0.367 [2,] 85 0.296 [3,] 122 0.260 [4,] 192 0.244 [5,] 275 0.175 [6,] 421 0.140 [7,] 603 0.093 [8,] 831 0.068 [9,] 1140 0.043 With the nls()-function i want to fit following formula whereas a,b, and c are variables: y~1/(a*x^2+b*x+c) With the standardalgorithm
2008 Jul 08
2
nls and "plinear" algorithm
hello all i havnt had a chance to read through the references provided for the "nls" function (since the libraries are closed now). can anyone shed some light on how the "plinear" algorithm works? also, how are the fitted values obtained? also, WHAT DOES THE ".lin" below REPRESENT? thanking you in advance ###################################### i have a quick
2008 Apr 29
0
nls plinear formula
I want to fit a nonlinear model of the form: Y=A+B*X1+C*X2+log(X3/(X3+D)) I think that the best way is to use the plinear algorithm, but I don't know how to specify the formula in the nls function. I've tried: Y~cbind(rep(1,times=length(Y)),X1,X2,log(X3/(X3+D))) But this fits the model: Y=A+B*X1+C*X2+D1*log(X3/(X3+D)) How can I specify the formula correctly? Thanks, M
2006 Jun 09
0
R CMD check and directory/package name
The NEWS for R 2.3.0 states that "R CMD check works for packages whose package name is different from the directory name in which it is located." However that hasn't been my experience. I ran R CMD check on package sources located in a directory with the same name as the package and it worked as expected. Then I renamed the directory and tried again. The first attempt got stuck
2015 Apr 29
0
Formula evaluation, environments and attached packages
Hi Milan, I expect I may be able to do something about the way the terms are evaluated, to ensure the evaluation is done in the gnm namespace (while still ensuring the variables can be found!). In the meantime, I think the following will work: Mult <- gnm::Mult f <- Freq ~ Eye + Hair + Mult(Eye, Hair) gnm::gnm(f, family=poisson, data=dat) Hope that helps, Heather On Wed, Apr 29, 2015,
2015 Apr 29
2
Formula evaluation, environments and attached packages
Hi! Some time ago, I replaced calls to library() with calls to requireNamespace() in my package logmult, in order to follow the new CRAN policies. But I just noticed it broke jackknife/bootstrap using several workers via package parallel. The reason is that I'm running model replicates on the workers, and the formula includes non-standard terms like Mult() which are provided by gnm. If gnm
2009 Apr 30
0
Categorical variable in a custom nonlin function with gnm
Hi all I want to construct a generalised nonlinear model (binomial family) using gnm, of the form: Response = a + b variable1 + c variable2 + d variable3 - d b variable4 - d c variable5, with the parameters b, c, and d appearing more than once. Hence, I think I need to use a custom nonlin function with gnm. One of my predictor variables is categorical, so I have created a dummy variable for
2013 Aug 22
1
Confusion about Depends:, Imports:, Enhances:, import(), inportFrom()
In checking my vcdExtra package, the following NOTE newly appeared (R-Forge, using R version 3.0.1 Patched (2013-08-20 r63635)) Package in Depends field not imported from: ?gnm? These packages needs to imported from for the case when this namespace is loaded but not attached. In the DESCRIPTION file, I have Depends: R (>= 2.10), vcd, gnm (>= 1.0.3) In NAMESPACE: # we are a vcd
2006 Nov 18
1
Why SAM has totally diffent results in R2.1.1 and R2.4.0
Hi, I am using SAM (from siggenes_1.2.11 package) method to select genes from a microarray data set. After installing the latest R2.4.0 on my computer, to my surprise the results are totally different from that calculated using R2.1.1. Even the example code doesn't work the same way under these two versions of R. Does anybody know what is going on? Thanks for any suggestions.
2007 Jan 16
1
nonlinear regression: nls, gnls, gnm, other?
Hi all, I'm trying to fit a nonlinear (logistic-like) regression, and I'd like to get some recommendations for which package to use. The expression I want to fit is something like: y ~ A * exp(X * Beta1) / (1 + exp(-(x + X * Beta2 - xmid)/scal)) Basically, it's a logistic function, but I want to be able to modify the saturation amplitude by a few parameters (Beta1) and shift the
2012 Aug 23
1
NLS bi exponential Fit
Hi everyone, I'm trying to perform a bi exponential Fit with the package NLS. the plinear algorithm seems to be a good choice see: p<-3000 q<-1000 a<--0.03 b<--0.02 t<-seq(0:144);t y<-p*exp(a*t) + q*exp(b*t)+rnorm(t,sd=0.3*(p* exp(a*t) + q*exp(b*t))) fittA <- nls(y~cbind(exp(a*t), exp(b*t)), algorithm="plinear",start=list(a=-.1, b=-0.2), data=list(y=y, t=t),
2006 Jan 08
1
confint/nls
I have found some "issues" (bugs?) with nls confidence intervals ... some with the relatively new "port" algorithm, others more general (but possibly in the "well, don't do that" category). I have corresponded some with Prof. Ripley about them, but I thought I would just report how far I've gotten in case anyone else has thoughts. (I'm finding the code
2010 Apr 15
0
[R-pkgs] vcdExtra 0.5-0 is released to CRAN
I'm pleased to announce the release of the vcdExtra package, v. 0.5-0 from R-Forge to CRAN, on its way to a CRAN server near you. vcdExtra was originally designed to serve as a sandbox for introducing extensions of mosaic plots and other visualizations for categorical data, particularly those that apply to (poisson surrogate) loglinear models fitted using glm() and related, generalized
2012 Sep 30
0
New package: logmult (log-multiplicative models)
This is a wrapper around gnm by Turner and Firth to make fitting log-multiplicative models as convenient as possible: it provides simple functions, good starting values, jackknife or bootstrap standard errors, and direct plotting of the results. In addition, it makes it possible to identify scores from RC(M) association models, which gnm does not allow without computing the SVD yourself. It will
2012 Sep 30
0
New package: logmult (log-multiplicative models)
This is a wrapper around gnm by Turner and Firth to make fitting log-multiplicative models as convenient as possible: it provides simple functions, good starting values, jackknife or bootstrap standard errors, and direct plotting of the results. In addition, it makes it possible to identify scores from RC(M) association models, which gnm does not allow without computing the SVD yourself. It will