Re: [R] install problem from Github or downloaded file

2017-03-01 Thread Jianling Fan
ty", not having been through CRAN, might simply be broken? >> >> Yes the directories 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

[R] install problem from Github or downloaded file

2017-03-01 Thread Jianling Fan
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

Re: [R] about data problem

2016-09-20 Thread Jianling Fan
inimal, self-contained, reproducible code. >> >> >> >> -- >> Gregory (Greg) L. Snow Ph.D. >> 538...@gmail.com >> > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org m

Re: [R] means by year, month and day

2016-07-17 Thread Jianling Fan
> > __ > 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,

Re: [R] mixed-effects models with (g)lmer in R and model selection

2016-02-19 Thread Jianling Fan
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.

[R] LSD value

2016-01-18 Thread Jianling Fan
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

Re: [R] LSD value

2016-01-18 Thread Jianling Fan
LSD > [1] 13.4704 >> out$statistics[4] > LSD > 13.4704 >> out[[1]][4] > LSD > 13.4704 > > ----- > David L Carlson > Department of Anthropology > Texas A University > College Station, TX 77840-4352 > > > -Orig

[R] find max. value for a equation by optimize

2015-10-29 Thread Jianling Fan
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) {

Re: [R] find max. value for a equation by optimize

2015-10-29 Thread Jianling Fan
> > On Thu, Oct 29, 2015 at 4:19 PM, Jianling Fan <fanjianl...@gmail.com> wrote: >> >> 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 >> opt

Re: [R] How to coerce a parameter in nls?

2015-09-22 Thread Jianling Fan
t;- c(Rm1=1.01, Rm2=1.01, Rm3=1.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 <fanjianl...@gmail.com> > wrote: >> >

Re: [R] How to coerce a parameter in nls?

2015-09-22 Thread Jianling Fan
will probably mean 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 >> packag

Re: [R] How to coerce a parameter in nls?

2015-09-22 Thread Jianling Fan
start = c(Rm1=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 > <ggrothendi...@gmail.com> wrote: >> >> Just write out the 20 terms. >> >> On Mon, Sep 21, 2015 at 10:

Re: [R] How to coerce a parameter in nls?

2015-09-21 Thread Jianling Fan
60 0.8480684 6 > 2980 0.9268951 6 > 30 100 0.9723207 6 > 31 120 0.9939966 6 > 32 140 0.9992400 6" > > dproot <- read.table(text = Lines, header = TRUE) > > > > On Mon, Sep 21, 2015 at 12:22 PM, Jianling Fan <fanjianl...@gmai

Re: [R] How to coerce a parameter in nls?

2015-09-21 Thread Jianling Fan
gt; 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

Re: [R] How to coerce a parameter in nls?

2015-09-21 Thread Jianling Fan
; 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 > 2860 0.8480684 6 > 2980 0.9268951 6 > 30 100 0.9723207 6 > 31 120 0.9939966 6 > 32 140 0.9

Re: [R] How to coerce a parameter in nls?

2015-09-20 Thread Jianling Fan
.OO#. .OO#. rocks...1k > --- > Sent from my phone. Please excuse my brevity. > > On September 18, 2015 9:54:54 PM PDT, Jianling Fan <fanjianl...@gmail.com> > wrote: >>Hello, everyone, >> >>I am using a nls regression with 6 groups data. I am trying

[R] How to coerce a parameter in nls?

2015-09-18 Thread Jianling Fan
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

[R] How to do global curve fitting in R

2015-09-09 Thread Jianling Fan
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

Re: [R] How to do global curve fitting in R

2015-09-09 Thread Jianling Fan
or on a statistics site, 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

Re: [R] How to do global curve fitting in R

2015-09-09 Thread Jianling Fan
ergence 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 <fanjianl...@gmail.com> wrote: >>

[R] nls model parameter compare?

2015-08-25 Thread Jianling Fan
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

[R] About nls.

2015-07-30 Thread Jianling Fan
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