_
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reprod
.
>>
>>
>>
>> --
>> Gregory (Greg) L. Snow Ph.D.
>> 538...@gmail.com
>>
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> h
Hello,
I am trying to do a nls regression with R. but I always get a error
as "Error in numericDeriv(form[[3L]], names(ind), env) : Missing
value or an infinity produced when evaluating the model". I googled
it and found someone said it is because of the improper start value. I
tried many times
Hello everyone,
I am doing nonlinear regression using a same sigmoidal model for
different treatments. for each treatment, I got a set of estimated
parameters (a1, b1, c1 for treatment 1; a2, b2, c2 for treatment 2;
a3, b3, c3 for treatment 3). And I want to compare these parameters
for different
Hello all,
I am trying to fit my data to a nls model. I have many sets of data
and each can fit well for the curve. but I want to fit them at once by
sharing 2 of 3 parameters of the model. I know it is a typical global
curve fitting problem, but I don't know how to do it by R?
Does anyone know a
te, like
> stats.stackexchange.com.
>
> Cheers,
> Bert
> Bert Gunter
>
> "Data is not information. Information is not knowledge. And knowledge
> is certainly not wisdom."
>-- Clifford Stoll
>
>
> On Wed, Sep 9, 2015 at 9:04 AM, Jianling Fan wrote:
>
convergence: 5
> Achieved convergence tolerance: 4.536e-06
>
> (See the references listed in help(nls) to learn more about the statistical
> issues.)
>
>
> Bill Dunlap
> TIBCO Software
> wdunlap tibco.com
>
> On Wed, Sep 9, 2015 at 9:28 AM, Jianling Fan wrote:
>>
>
Hello, everyone,
I am using a nls regression with 6 groups data. I am trying to coerce
a parameter to 1 by using a upper and lower statement. but I always
get an error like below:
Error in ifelse(internalPars < upper, 1, -1) :
(list) object cannot be coerced to type 'double'
does anyone know h
phone. Please excuse my brevity.
>
> On September 18, 2015 9:54:54 PM PDT, Jianling Fan
> wrote:
>>Hello, everyone,
>>
>>I am using a nls regression with 6 groups data. I am trying to coerce
>>a parameter to 1 by using a upper and lower statement. but I always
>
; would consider giving this a try.
>
> JN
>
>
> On 15-09-20 01:19 PM, Jianling Fan wrote:
>>
>> no, I am doing a regression with 6 group data with 2 shared parameters
>> and 1 different parameter for each group data. the parameter I want to
>> coerce is for o
360 0.950 5
> 2480 0.995 5
> 25 100 1.000 5
> 2620 0.4345774 6
> 2740 0.6654726 6
> 2860 0.8480684 6
> 2980 0.9268951 6
> 30 100 0.9723207 6
> 31 120 0.9939966 6
> 32 140 0.9992400 6"
>
> dproot &l
348 4
> 19 100 6.5780660 4
> 20 120 6.6032064 4
> 2120 0.614 5
> 2240 0.827 5
> 2360 0.950 5
> 2480 0.995 5
> 25 100 1.000 5
> 2620 0.4345774 6
> 2740 0.6654726 6
> 28 60 0.8480684 6
&
=1.01, Rm2=1.01, Rm3=1.01, Rm4=6.65, Rm5=1.01, Rm6=1,
> d50=20, c=-1))
>
>
> On Tue, Sep 22, 2015 at 7:04 AM, Gabor Grothendieck
> wrote:
>>
>> Just write out the 20 terms.
>>
>> On Mon, Sep 21, 2015 at 10:26 PM, Jianling Fan
>> wrote:
>>>
>>&
having a "translation" routine to glue things
> together.
>
> JN
>
>
> On 15-09-21 12:22 PM, Jianling Fan wrote:
>>
>> Thanks Prof. Nash,
>>
>> Sorry for late reply. I am learning and trying to use your nlmrt
>> package since I got your emai
.01, Rm4=6.65, Rm5=1.01, d50=20, c=-1)
> Rm6 <- 1
>
> fm.nls <- nls(fo, dproot2, start = st)
>
> AIC(fm.nls)
> summary(fm.nls)
>
>
> On Tue, Sep 22, 2015 at 12:46 PM, Jianling Fan
> wrote:
>>
>> Hello Prof. Nash,
>>
>> My regression works goo
Hello, everyone,
I have a specific equation, When I draw it out, I found there is a
max. value for it at about x=5. but I couldn't point it out by
optimize(). Does anyone know why? Thanks!
my code is:
d50<-18.04
c<- -1.276
dm<-147
sm <- function (x) {
1/(1+(x/d50)^c)+(1-1/(1+(dm/d50)^c))*x/
Oops!
Thanks a lot! Jean
On 29 October 2015 at 15:45, Adams, Jean wrote:
> If you read the help file for optimize, you will see that you need to set
> the argument maximum to TRUE.
>
> ?optimize
>
> optimize(sm, c(0, 50), maximum=TRUE)
>
> Jean
>
> On Thu, Oct 29,
Hello, everyone,
I am using LSD.test() from package "agricolae" to do my anova
analysis. I know I can calculate LSD value by its equation
t*sqrt(MSE*2/n), but I am wondering if there is a code or something
that can give the LSD value more directly in R?
Thanks!
Julian
__
out$statistics[4]
> LSD
> 13.4704
>> out[[1]][4]
> LSD
> 13.4704
>
> -----
> David L Carlson
> Department of Anthropology
> Texas A&M University
> College Station, TX 77840-4352
>
>
> -Original
Hello, Wilbert,
You did give a good procedure for lme model selection! thanks! I learn some.
I am also working on similar problem recently, maybe you can take a
look at "glmmLasso" package, which allows model selection in
generalized linear mixed effects models using the LASSO shrinkage
method.
Hello everyone,
I can install packages from CRAN, but I got some problem for
installing packages from Github or downloaded package files. I got
this problem in recent weeks.
example:
> devtools::install_github('bart6114/artyfarty')
Downloading GitHub repo bart6114/artyfarty@master
from URL https
tories you ask about are only present temporarily while the
>> package is checked and moved into your package library, so they are
>> "temporary" even though they are not in your TMP or TEMP directories. The
>> missing file is almost certainly the R file, not the
22 matches
Mail list logo