spatstat

    0熱度

    1回答

    我有一段距離序列,我想用spatstat繪製一條線。例如: library(spatstat) x <- c(0.3, 5) y <- c(3, 1.2) range.x <- c(0, max(x)+0.2) range.y <- c(0, max(y)+0.2) owin <- owin(range.x, range.y) the.line <- psp(x0 = x[1],

    6熱度

    2回答

    在R包spatstat(我使用的是當前版本,1.31-0)中,有一個選項use.gam。如果將此設置爲true,則可以在線性預測器中包含平滑項,這與使用R包mgcv的方法相同。例如, g <- ppm(nztrees, ~1+s(x,y), use.gam=TRUE) 現在,如果我想要攔截的置信區間,通常可以使用summary或vcov,當你不使用gam其工作,但失敗時,你使用GAM vcov

    0熱度

    1回答

    我有一個土地覆蓋的柵格文件,我已經減少到僅包含樹木覆蓋單元。我在柵格包中使用了clump來聚集()一起連續的森林區域。這使得所有單元格相互接觸相同的ID,因爲它們是同一個補丁的一部分。 然後我想找出每個叢塊的PatchStat(),我通過將我的叢塊柵格轉換爲as.matrix來完成。我試圖讓PatchStat()對柵格執行此操作,但只有在矩陣中才有效。 我現在想用補丁stat輸出做一個光柵,即「p

    0熱度

    2回答

    我試圖將文本文件導入到R,特別是spatstat包中。我已經加載了一個shp文件作爲窗口,並且工作得很好(顯示plot())。我不能讓ppp()命令運行。我一直運行pp <- ppp(X, Y, window=W) Error in ppp(X, Y, window = W) : 1 out of 904 coordinate values are NA or NaN 後收到此錯誤我

    0熱度

    1回答

    我仍然是R的初學者,並且存在以下問題。任何提示將高度讚賞 - 這將是一個很大的幫助開始。 我在ArcGIS中有一個國家的道路shapefile,我想用它作爲R中的網絡變量來應用網絡K函數(使用spatstat包查找網絡上的點隨機性)。我在R之前從來沒有用過這種類型的問題,我無法弄清楚如何將矢量線數據集轉換爲R中的網絡變量。我嘗試使用shapefile包,但在閱讀.shp,.shx和.dbf文件,我

    0熱度

    1回答

    我試圖按照教程"Handling shape files in the spatstat package" by Adrian Baddeley將SpatialPolygonsDataFrame的形狀轉換爲owin的集合。 形狀的關鍵屬性總結如下,我的目標是測試點的集合是否在多邊形集合內。 摘要(形狀) Object of class SpatialPolygonsDataFrame Coord

    0熱度

    1回答

    我希望測試均勻Thomas集羣進程的適合度, 理想情況下使用殘差圖,如diagnose.ppm和qqplot.ppm。 由於使用kppm的最簡單的方法適合托馬斯羣集過程是 ,我不能使用上述功能。任何人都可以建議 做到這一點? 例如轉化kppm到ppm或使用PPM,以適應托馬斯集羣過程中, 非常感謝

    0熱度

    1回答

    我正在使用包spatstat。我有一個座標數據框,我分爲兩組,健康(Mark =否)和患病(Mark =是)。我能找到患病點和所有其他點之間的距離: >D<-crossdist(diseased,healthy) 但現在我需要根據他們從患病點,有點像的距離此言我的觀點: >i<-length(D) >for n=1:i, > if D[n]<1, > mark(n)<-yes, >

    3熱度

    3回答

    所以我遇到的問題是非常基本的,但我似乎無法解決它。我一直在使用spatstat的F,G,J,K和L函數,我想在同一個繪圖上爲一個點過程繪製F函數,作爲不同點過程的F函數。這很容易,除了我的一個點進程評估函數r = 20,另一個r = 15('r'是x軸上的獨立變量)。我希望他們都被評估爲相同的r值,最好20.我該怎麼做?本質上,這是問題所在。 首先,我創建一個包含函數的值的數據幀: mydata.

    4熱度

    1回答

    這是我遇到的一個非常奇怪的情況。基本上,我試圖將累積分佈函數擬合到我的數據的G函數中。完成後,我想繪製模型和原始數據,並將其輸出爲PDF。我會讓代碼解釋(簡單的複製和粘貼): library(spatstat) data(swedishpines) mydata <- swedishpines mydata.Gest <- Gest(mydata) Gvalues <- mydat