私はあなたがこの問題を解決するのに役立つことを願っています。申し訳ありませんが、私はこの記事を書いている間にいくつかの間違いを犯した場合、私の英語は少し錆びています。LinearKのR関数を使ったタコの捕獲の分析R
ここに質問があります。私はRで分析したい.shpデータを持っています。.shpは、タコを捕まえるために設定したトラップ行を表す行か、それらの行の直上にあるポイントで、どこで成熟したかを表します。
私は答えることをしようとしている質問です:は、統計的にグループ化されたかどうかタコか?
調査の結果、RaptorとそのlinearK関数を使用して、Maptools、SpatStat、Spというライブラリを使用してその質問に答える必要があると思われます。
t1<- as.linnet(readShapeSpatial("./20170518/t1.shp"))
は、私は次の警告を取得し、それをトラックとlinnetがオブジェクトを作成するライブラリ
library(spatstat)
library(maptools)
library(sp)
をロード
:ここ
は私がRStudioで使用しているコードです。動作するように思われるWarning messages:
1: use rgdal::readOGR or sf::st_read
2: use rgdal::readOGR or sf::st_read
すべてが、私はここに同じ警告を受けるが、ときに私の実際の問題は、開始ポイント
p1<- as.ppp(readShapeSpatial("./20170518/p1.shp"))
とのPPPオブジェクトの作成
plot(t1)
OKであることを確認するためにそれをプロットそれをプロットしようとしてください:
> plot(p1)
Error in if (!is.vector(xrange) || length(xrange) != 2 || xrange[2L] < :
missing value where TRUE/FALSE needed
In addition: Warning messages:
1: Interpretation of arguments maxsize and markscale has changed (in spatstat version 1.37-0 and later). Size of a circle is now measured by its diameter.
2: In plot.ppp(x, ..., multiplot = FALSE, do.plot = FALSE) :
All mark values are NA; plotting locations only.
3: In plot.ppp(x, ..., multiplot = FALSE, do.plot = FALSE) :
All mark values are NA; plotting locations only.
4: In plot.ppp(x, ..., multiplot = FALSE, do.plot = FALSE) :
All mark values are NA; plotting locations only.
5: In plot.ppp(x, ..., multiplot = FALSE, do.plot = FALSE) :
All mark values are NA; plotting locations only.
6: In plot.ppp(x, ..., multiplot = FALSE, do.plot = FALSE) :
All mark values are NA; plotting locations only.
7: In plot.ppp(x, ..., multiplot = FALSE, do.plot = FALSE) :
All mark values are NA; plotting locations only.
今すぐ何残っていることは、LPPオブジェクト内のオブジェクトに参加するとlinearK機能
> pt1 <- lpp(p1,t1)
> linearK(pt1)
Function value object (class ‘fv’)
for the function r -> K[L](r)
......................................
Math.label Description
r r distance argument r
est {hat(K)[L]}(r) estimated K[L](r)
......................................
Default plot formula: .~r
where “.” stands for ‘est’
Recommended range of argument r: [0, 815.64]
Available range of argument r: [0, 815.64]
これが今の私の状況であると、それを分析することです。私が知らないのはなぜプロット関数がpppオブジェクトで動作しないのか、そしてlinearK関数のリターンを控える方法です。ヘルプ(linearK)はヒントを提供しませんでした。私は多くのトラックを持っているので、それぞれのポイントがあり、私の希望する結果は、分析されたxトラックのような要約であり、グループ化され、分散され、未知のものです。
ありがとうございました。私がこの問題を解決するのを手伝っていただければ幸いです。
編集:ここには、1日分のshpファイル、トラックとポイント、およびコード付きのtxtファイルを含むzipファイルへのリンクがあります。 https://drive.google.com/open?id=0B0uvwT-2l4A5ODJpOTdCekIxWUU
あなたははるかに例が完全に再現可能にするあなたの良い答えを得るために可能性が高いです。 – Axeman
あなたは、大丈夫です。私は私の質問の一番下にリンクを追加しました。 ありがとうございます。 –