[R] please help me

James W. MacDonald jmacdon at med.umich.edu
Thu Nov 8 23:16:29 CET 2007


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 at 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



More information about the R-help mailing list