Specifying the initial values of likfit function of geodata in R

[Cynderella]

I have a geostatistical data. For analysis, I am using `geoR`. I want to estimate the parameters. But `likfit` or `variofit` function requires initial values for the covariance parameters: \(\sigma^2\) (partial sill) and \(\phi\) (range parameter).

But I can't guess the initial values for my data set.

How can I specify the initial values?

 

[Dason]

Are you saying you don't know any value to plug in or you've tried some values and it doesn't end up converging?

 

[Cynderella]

Are you saying you don't know any value to plug in or you've tried some values and it doesn't end up converging?

My data is in thousand (say, 3503.198, 1821.964, etc.) and coordinates are latitude and longitude. If I plug any values, it converges. But the estimates change. Say I put initial value \(\sigma^2=1\) and \(\phi=0.5\). Then the estimated value of \(\sigma^2=0\). But if I plug \(\sigma^2=22000\) and \(\phi=0.5\), then the estimated value of \(\sigma^2\) is a very big number.