[R] geoR/variog4() not returning all directions

Paulo Justiniano Ribeiro Jr paulojus at est.ufpr.br
Thu Sep 16 22:44:06 CEST 2004


Hi Martin

The only pssible explanation I can find from your description is that the
code may not be fiuding pairs in one of the directions.

however...I never discard the possibility of a bug :)

Could you please send me the data and the piec of relevant code such that
i can rteproduce the error here?

Thanks
P.J,



On Thu, 16 Sep 2004, Martin Henry H. Stevens wrote:

> Mac OS 10.3.5, R 2.0.0
> latest version of geoR
>
> I have an incomplete 5 x 20 spatial array of samples (60 out of 100
> possible locations) for which I would like to calculate directional
> variograms using variog4(). Unfortunately, I can't get it to return all
> 4 directions. It returns variograms for 45, 90, and 135 degrees,
> omitting 0 degrees (pi/4, pi/2, 3pi/4, omitting 0). If I specify 0
> degrees, it returns a variogram with (binned) distances of 1, 2, 3, and
> 4.
>
> Any ideas what I am doing wrong, or what other information I could
> provide to better elucidate my problem?
> Many thanks, as usual.
> Hank Stevens
>
> Dr. Martin Henry H. Stevens, Assistant Professor
> 338 Pearson Hall
> Botany Department
> Miami University
> Oxford, OH 45056
>
> Office: (513) 529-4206
> Lab: (513) 529-4262
> FAX: (513) 529-4243
> http://www.cas.muohio.edu/botany/bot/henry.html
> http://www.muohio.edu/ecology/
> http://www.muohio.edu/botany/
> "E Pluribus Unum"
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
>

Paulo Justiniano Ribeiro Jr
Departamento de Estatística
Universidade Federal do Paraná
Caixa Postal 19.081
CEP 81.531-990
Curitiba, PR  -  Brasil
Tel: (+55) 41 361 3573
Fax: (+55) 41 361 3141
e-mail: paulojus at est.ufpr.br
http://www.est.ufpr.br/~paulojus

 /"\
 \ /  Campanha da fita ASCII - contra mail html
  X   ASCII ribbon campaign - against html mail
 / \




More information about the R-help mailing list