1

再捕獲データの Minimum Number Alive (MNA) の計算を手伝ってくれる人がいるかどうか疑問に思っています。

MNA は次のように計算されます。

時点 (t) で捕獲された個体の実際の数 + 存在し、時点 (t) で捕獲されなかったが、その後捕獲された個体 (Krebs, 1966)

101 個体の 44 の二次捕獲期間で構成される 9 つの一次捕獲期間のバイナリ データ (捕獲の場合は 1、非捕獲の場合は 0) があります。セカンダリ期間のキャプチャ データを、プライマリ キャプチャ期間ごとに 0 または 1 に要約する必要があります。このようにして、各プライマリ キャプチャ期間の MNA を計算できます。

とても感謝しております、

ターニャ

4

1 に答える 1

0

これはあなたの二次期間のデータを凝縮しているようです。ただし、動物 6 と 10 は、検出されなかったにもかかわらず、第 2 原始期に生きていることが知られていることに注意してください。また、動物 5 は、最初の初生期の後に検出されたという理由だけで、最初の初生期に生存しているとは知られていません...

この例では、3 つのプライマリ ピリオドと 3 つのセカンダリ ピリオドがあります。

地理的および人口学的閉鎖に関するあなたの仮定が不明であるため、最小生存数を推定しようとはしていません. ただし、以下のデータの概要が役立つことを願っています。

my.data <- as.data.frame(matrix(rbinom(10*9,1,0.2), nrow = 10, byrow=FALSE))
my.data

#   V1 V2 V3 V4 V5 V6 V7 V8 V9
#1   0  0  0  0  0  0  0  1  1
#2   0  0  0  0  0  0  0  0  0
#3   0  0  0  0  0  0  0  0  0
#4   0  0  1  0  0  0  0  0  0
#5   0  0  0  0  0  1  0  1  1
#6   0  1  0  0  0  0  0  0  1
#7   1  1  1  0  0  0  0  0  0
#8   1  0  0  0  1  0  0  0  0
#9   0  0  0  0  1  0  1  1  1
#10  1  0  0  0  0  0  0  0  1

year.sums <- sapply(1:3, function(i)  apply(my.data[,(3*(i-1)+1):(3*(i-0))], 1, function(x) sum(x)) )
year.sums

#      [,1] [,2] [,3]
# [1,]    0    0    2
# [2,]    0    0    0
# [3,]    0    0    0
# [4,]    1    0    0
# [5,]    0    1    2
# [6,]    1    0    1
# [7,]    3    0    0
# [8,]    1    1    0
# [9,]    0    1    3
#[10,]    1    0    1

year.sums[year.sums>0] <- 1
year.sums

#      [,1] [,2] [,3]
# [1,]    0    0    1
# [2,]    0    0    0
# [3,]    0    0    0
# [4,]    1    0    0
# [5,]    0    1    1
# [6,]    1    0    1
# [7,]    1    0    0
# [8,]    1    1    0
# [9,]    0    1    1
#[10,]    1    0    1
于 2013-09-20T05:47:04.717 に答える