Well Sadeghian, I have to give you credit for perseverance. This must be the fifth email you have sent to R-help for this one question.
Unfortunately, with each additional message, your chances of getting any help are becoming asymptotically closer to zero. So some advice: > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. I think you definitely need to read the posting guide, and then maybe you could try to come up with some code using a small portion of your data set that would illustrate your problem? As it stands, nobody can attempt to reproduce your error to see what the problem is, so how can anybody possibly help? Best, Jim azadeh sadeghian wrote: > Dear list. > I am student M.S. statistics in department statistics .I have thise problem > in work with the function variofit and nls and dont know how to solve it. > var1<-variog(data,option="bin") > var2<-variog(data,option="cloud") > v1<-var1$v > u1<-var1^u > v2<-var2$v > u2<-var2$u > variofit(var1,ini.cov.pars=c(0.005,1.5),cov.model="power",fix.nugget=F,weight="equal") > variofit: weights used: equal > variofit: minimisation function used: optim > Error in if (loss > (.Machine$double.xmax^0.5) | loss == Inf | loss == : > missing value where TRUE/FALSE needed > In addition: Warning message: > unreasonable initial value for sigmasq + nugget (too low) in: variofit(var1, > ini.cov.pars = c(0.005, 1.5), cov.model = "power", > ....................................................................................................................... > variofit(var2,ini.cov.pars=c(0.005,1.5),cov.model="power",fix.nugget=F,weight="equal") > variofit: weights used: equal > warning: minimisation function nls can not be used with given cov.model. > changing for "optim". > variofit: minimisation function used: optim > Error in if (loss > (.Machine$double.xmax^0.5) | loss == Inf | loss == : > missing value where TRUE/FALSE needed > In addition: Warning message: > In variofit(var2, ini.cov.pars = c(0.005, 1.5), cov.model = "power", : > unreasonable initial value for sigmasq + nugget (too low) > ............................................................................................................................... > nls(v2~c0+ce*(1-exp(-(u2^2)/(ae^2))),start=list(c0=0,ae=30000,ce=20000)) > Error in nls(v2 ~ c0 + ce * (1 - exp(-(u2^2)/(ae^2))), start = list(c0 = 0, > : > number of iterations exceeded maximum of 50 > best regards. > Sadeghian. > > __________________________________________________ > > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > 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, reproducible code. -- James W. MacDonald, M.S. Biostatistician Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 ______________________________________________ R-help@r-project.org mailing list 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, reproducible code.
