[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