<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-AU link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Dear list,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am aware there are many post dealing with singular model variogram fits but I have not been able to fix my problem even after re-scaling the coordinates and variables.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am using fit.variogram(gstat) on the file ‘trans.csv’ (attached) and have tried several values for the range, sill and nugget and different models. The following gives me an ok variogram (from visual plot) but still a ‘singular model’ warning. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>library(sp)<o:p></o:p></p><p class=MsoNormal>library(gstat)<o:p></o:p></p><p class=MsoNormal>points <- read.csv("trans.csv")<o:p></o:p></p><p class=MsoNormal>coordinates(points)=~x+y<o:p></o:p></p><p class=MsoNormal>depth.ivgm <- vgm(model='Gau',range=1000, psill=var(points$ln_depth))<o:p></o:p></p><p class=MsoNormal>depth.vgm <- fit.variogram(variogram(ln_depth~x+y, points), model=depth.ivgm)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Any ideas on why I am getting a singular fit? Is it the small number of samples used?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Cheers,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Javier<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>