0

R の mlogit パッケージで動作するようにデータをセットアップしようとしています

次のコードで作成されたデータフレームがあります。

id <- 1:10
id <- rep(id, each=5)
site <- c("site1", "site2", "site3", "site4", "site5")
choice <- c("site3", "site5", "site1", "site4", "site2",
            "site4", "site3", "site5", "site2", "site1")
df <- cbind(id, site)

id の値ごとにサイトの選択を示すバイナリ変数を作成したいと考えています。id 変数は繰り返されるシーケンスであるため、新しい指標変数は、「site」が「choice」の関連する値と同等である行を除くすべての行で 0 である必要があります。id == 1 の場合、これは「choice」ベクトルの最初の要素になります。id == 2 の場合、選択ベクトルの 2 番目の要素になります。

変数が含まれている最終的なデータフレームは次のようになります。

      id   site   indicator
 [1,] "1"  "site1" "0"
 [2,] "1"  "site2" "0"
 [3,] "1"  "site3" "1"
 [4,] "1"  "site4" "0"
 [5,] "1"  "site5" "0"
 [6,] "2"  "site1" "0"
 [7,] "2"  "site2" "0"
 [8,] "2"  "site3" "0"
 [9,] "2"  "site4" "0"
[10,] "2"  "site5" "1"
[11,] "3"  "site1" "1"
[12,] "3"  "site2" "0"
[13,] "3"  "site3" "0"
[14,] "3"  "site4" "0"
[15,] "3"  "site5" "0"
[16,] "4"  "site1" "0"
[17,] "4"  "site2" "0"
[18,] "4"  "site3" "0"
[19,] "4"  "site4" "1"
[20,] "4"  "site5" "0"
[21,] "5"  "site1" "0"
[22,] "5"  "site2" "1"
[23,] "5"  "site3" "0"
[24,] "5"  "site4" "0"
[25,] "5"  "site5" "0"
[26,] "6"  "site1" "0"
[27,] "6"  "site2" "0"
[28,] "6"  "site3" "0"
[29,] "6"  "site4" "1"
[30,] "6"  "site5" "0"
[31,] "7"  "site1" "0"
[32,] "7"  "site2" "0"
[33,] "7"  "site3" "1"
[34,] "7"  "site4" "0"
[35,] "7"  "site5" "0"
[36,] "8"  "site1" "0"
[37,] "8"  "site2" "0"
[38,] "8"  "site3" "0"
[39,] "8"  "site4" "0"
[40,] "8"  "site5" "1"
[41,] "9"  "site1" "0"
[42,] "9"  "site2" "1"
[43,] "9"  "site3" "0"
[44,] "9"  "site4" "0"
[45,] "9"  "site5" "0"
[46,] "10" "site1" "1"
[47,] "10" "site2" "0"
[48,] "10" "site3" "0"
[49,] "10" "site4" "0"
[50,] "10" "site5" "0"

私はこれを何度も試みましたが、それを理解することも、オンラインで関連する答えを見つけることもできません.

前もって感謝します ：）

4

3 に答える 3

0

指定された行列 (df) を使用して、指標は次のように計算できます。

indicator <- as.numeric(choice[as.numeric(df[,"id"])] == df[,"site"])
# Final matrix 
df <- cbind(df,indicator)

必要なマトリックス:

id   site    indicator
 [1,] "1"  "site1" "0"      
 [2,] "1"  "site2" "0"      
 [3,] "1"  "site3" "1"      
 [4,] "1"  "site4" "0"      
 [5,] "1"  "site5" "0"      
 [6,] "2"  "site1" "0"      
 [7,] "2"  "site2" "0"      
 [8,] "2"  "site3" "0"      
 [9,] "2"  "site4" "0"      
[10,] "2"  "site5" "1"      
[11,] "3"  "site1" "1"      
[12,] "3"  "site2" "0"      
[13,] "3"  "site3" "0"      
[14,] "3"  "site4" "0"      
[15,] "3"  "site5" "0"      
[16,] "4"  "site1" "0"      
[17,] "4"  "site2" "0"      
[18,] "4"  "site3" "0"      
[19,] "4"  "site4" "1"      
[20,] "4"  "site5" "0"      
[21,] "5"  "site1" "0"      
[22,] "5"  "site2" "1"      
[23,] "5"  "site3" "0"      
[24,] "5"  "site4" "0"      
[25,] "5"  "site5" "0"      
[26,] "6"  "site1" "0"      
[27,] "6"  "site2" "0"      
[28,] "6"  "site3" "0"      
[29,] "6"  "site4" "1"      
[30,] "6"  "site5" "0"      
[31,] "7"  "site1" "0"      
[32,] "7"  "site2" "0"      
[33,] "7"  "site3" "1"      
[34,] "7"  "site4" "0"      
[35,] "7"  "site5" "0"      
[36,] "8"  "site1" "0"      
[37,] "8"  "site2" "0"      
[38,] "8"  "site3" "0"      
[39,] "8"  "site4" "0"      
[40,] "8"  "site5" "1"      
[41,] "9"  "site1" "0"      
[42,] "9"  "site2" "1"      
[43,] "9"  "site3" "0"      
[44,] "9"  "site4" "0"      
[45,] "9"  "site5" "0"      
[46,] "10" "site1" "1"      
[47,] "10" "site2" "0"      
[48,] "10" "site3" "0"      
[49,] "10" "site4" "0"      
[50,] "10" "site5" "0"
于 2016-04-07T14:02:01.683 に答える
0

Akrun が提案したようdata.frameに、df を定義するために使用します。

df <- data.frame(id, site)

次に、次のようにします。

df$indicator <- (df$site == choice[df$id])*1

*1T/F の結果を 1 と 0 に変換します。

結果：

   id  site indicator
1   1 site1         0
2   1 site2         0
3   1 site3         1
4   1 site4         0
5   1 site5         0
6   2 site1         0
7   2 site2         0
8   2 site3         0
9   2 site4         0
10  2 site5         1
11  3 site1         1
12  3 site2         0
13  3 site3         0
14  3 site4         0
15  3 site5         0
16  4 site1         0
17  4 site2         0
18  4 site3         0
19  4 site4         1
20  4 site5         0
21  5 site1         0
22  5 site2         1
23  5 site3         0
24  5 site4         0
25  5 site5         0
26  6 site1         0
27  6 site2         0
28  6 site3         0
29  6 site4         1
30  6 site5         0
31  7 site1         0
32  7 site2         0
33  7 site3         1
34  7 site4         0
35  7 site5         0
36  8 site1         0
37  8 site2         0
38  8 site3         0
39  8 site4         0
40  8 site5         1
41  9 site1         0
42  9 site2         1
43  9 site3         0
44  9 site4         0
45  9 site5         0
46 10 site1         1
47 10 site2         0
48 10 site3         0
49 10 site4         0
50 10 site5         0

as.character変換したい列で数字や係数の代わりに文字列を使用したい場合

于 2016-04-07T13:09:53.030 に答える
0

split「id」で「サイト」を取得し、「choice」の対応する値と比較して論理インデックスを取得するために使用できMapます。

df$indicator <- +(unlist(Map(`==`, split(df$site, df$id), choice), use.names=FALSE))

または で「id」の頻度を取得しtabulate、「choice」を複製し、「site」と比較してバイナリに変換します

+(rep(choice,tabulate(df$id))==df$site)

データ

df <- data.frame(id, site)
于 2016-04-07T12:14:11.373 に答える