2017-06-14 108 views
0

幫助density`功能「時間迭代沒有超過iterMax;增加DTT或增加iterMax」中的R爲density功能得到錯誤`spatstat ::有關「iterMax」錯誤

如何增加iterMax?

我正在使用庫spatstat。

這是我的腳本:

D660 < - 密度(取消標記(control23_network.ppp),660)

線性網絡上我的觀點模式有:

1969年點 線性網絡20126頂點和21363線 圍窗口:矩形= [-87.63141,-87.55547]×[41.75817,41.83858]單元

enter image description here

它的工作的一個例子是在與樣品數據集中的spatstat庫:

數據(芝加哥)

芝加哥

D60 < - 密度(取消標記(芝加哥),60)

plot(d60)

d60

該樣品數據集具有:線性網絡 116點

點圖案 多類型,與可能的類型:攻擊,盜竊,cartheft,損壞,搶劫,盜竊,侵入 線性網絡338點的頂點和503線 圍窗口:長方形= [0.3894,1281.5523] x [153.1035,1276.5602]英尺

+0

有兩件事要嘗試:首先是發佈[MCVE]。其次,爲了增加SO的知識而不太理想,可以嘗試在調用「spatstat :: density.lpp」時添加iterMax。 (在R中大寫函數名稱不正確會導致誤導新手的不良副作用。) –

回答

1

這是一個關於spatstat::density.lpp的問題,density針對線性網絡上的點模式(類別lpp)的方法。

選擇的帶寬sigma太大了。包圍矩形大約爲0。1個單位寬,而帶寬是sigma=660單位。這將需要大量的算法迭代,所以它會拒絕。

在這個例子中,帶寬sigma的合理值將在0.001和0.05單位之間。

1

平滑帶寬sigma與網絡規模(包圍矩形的近似邊長爲0.1個單位)相比似乎很龐大(660個單位)。你有沒有試過類似sigma = 0.001或類似的東西?值得注意的是sigma註定會出現問題。幫助說:

計算時間短,但與西格瑪二次增加。

也許這不是根本問題,但至少在做很多其他事情之前值得嘗試。

最後,你的R和spatstat版本是什麼? density.lpp的代碼最近改變了很多。