🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

Q&A

1回答

941閲覧

「ダミー変数の1」を違う数値に置き換えたい

退会済みユーザー

退会済みユーザー

総合スコア0

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

0グッド

0クリップ

投稿2020/12/10 08:18

編集2020/12/10 08:28

前提・実現したいこと

「R」で経度や緯度が入っているデータを用いて距離計算をしています。
データの内容はある国の世帯データです。
距離計算を行い、他の世帯が200m以内であれば1、それ以上であれば0が入るダミー変数を作成しました。

最終的には、ダミー変数に1を各家計の所得に置き換えたいです。

例)A,B,C,D,Eの5世帯があったとして、Aさんの家から200m以内の家はB,Cだとします。
B,Cには1が入り、D,Eには0が入ります。
現状、ここまでは以下のコードで完成しています

これから、200m以内の家計(ダミーに1が入っている家計)の所得が知りたいので
1を所得に置き換えたいです。

現状、このような形でダミー変数を作成しています。
この1を所得に置き換えたいです。
||所得|A|B|C|D|E|
|:--|:--:|--:|
|A|1000|1|1|1|0|0
|B|2000|1|1|1|0|0
|C|2000|1|1|1|0|0
|D|1500|0|0|0|1|1
|E|1000|0|0|0|1|1

完成イメージ
||所得|A|B|C|D|E|
|:--|:--:|--:|
|A|1000|1000|2000|2000|0|0
|B|2000|1000|2000|2000|0|0
|C|2000|1000|2000|2000|0|0
|D|1500|0|0|0|1500|1000
|E|1000|0|0|0|1500|1000

発生している問題・エラーメッセージ

ダミー変数を作成したところで止まっています。

該当のソースコード

#200m以内にある世帯のダミーを作成 for (i in 1:28){ for (j in 1:28){ colname <- gps_v1_v2[i, "hhid"] gps_v1_v2[j,colname ] <- distGeo(gps_v1_v2[i, c("longitude", "latitude")], gps_v1_v2[j, c("longitude", "latitude")]) gps_v1_v2[j,colname ] <- ifelse(gps_v1_v2[j,colname ] <= 200,1, 0) } }

試したこと

<= 200,1, 0)
の1の部分に所得データの変数名である”income”
と入力しましたが、1の所に「1行目の世帯の所得」がすべての1に代入されてしまいました。

補足情報(FW/ツールのバージョンなど)

ここにより詳細な情報を記載してください。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

KojiDoi

2020/12/10 11:57

質問に書いてあるデータと該当のソースコードで参照しているデータの内容が食い違っているため、何をしたいのかがわかりません。矛盾を解消し、質問と直接関係ない部分は省いてください。
退会済みユーザー

退会済みユーザー

2020/12/10 12:50

ご指摘ありがとうございます。 まず、表はあくまでも例で掲載したものになります(実際のデータとは異なります)。 実際に使用してるデータ(gps_v1_v2)は世帯(hhid)が28世帯入っているデータになります。 そこに経度、緯度があるので距離計算コード(distGeo)で世帯間の距離を計算し 距離が200m以内であれば1が入るダミー変数を作成しています。 経度、緯度以外に28世帯の所得が入っている列もあるため、ダミー変数で1が入っているところを各世帯の所得に置き換えることは可能なのか。 という質問でございます。 日本語が拙く分かりにくい質問となってしまい申し訳ありません。 よろしくお願いいたします。
KojiDoi

2020/12/11 05:01

指摘したことを全く理解していただいてないようで、残念ですね。 言い方を変えます。 colname <- gps_v1_v2[i, "hhid"] このcolnameにはどのようなデータが入っているのですか? その次の行で、gps_v1_v2[j,colname ]という記述があるので、普通に考えれば列番号ですが、そんなものが元データに格納されているのは不自然なので、たぶん違うのでしょう。しかし、このようにコードに書いてあるものと質問の内容が矛盾しているので、結局、対象のデータがどのような形をしているのか、読んでいる方には全くわかりません。 >ダミー変数で1が入っているところを各世帯の所得に置き換えることは可能なのか。 可能か? といわれたら、たぶん可能なのでしょう。しかし、上記のような問題から、どうすれば要求を満たすコードとなるのか判断できません。日が変わってもだれからも回答がつかないのはそれが原因です。 私は本物のデータを丸ごと示せと言っているのではありません。しかし、ダミーデータであっても、書いてみたコードや質問の内容を過不足なく反映しているものを見せてくれないと、どうしようもありません。我々はエスパーではないので。 というわけで、質問が質問の体を為していないのです。修正してください。
guest

回答1

0

質問内容の例のデータフレームのうち,一つ目が完成しているとして
ifelse()を使って,次のように実装しました.

質問では,ifelse()に入れるvectorの長さが異なっているのが原因では?
と感じました.

R

1#データ作成 2alphabet<-c("A","B","C","D","E") 3Data<-data.frame(matrix(c(1000,2000,2000,1500,1000, 4 1,1,1,0,0, 5 1,1,1,0,0, 6 1,1,1,0,0, 7 0,0,0,1,1, 8 0,0,0,1,1),5,6,dimnames = list(alphabet,c("所得",alphabet))))

R

1#行単位でみていく 2for(i in 1:nrow(Data)){ 3 Data[i,2:ncol(Data)]<-ifelse(Data[i,2:ncol(Data)]==1,Data$所得,0) 4}

投稿2020/12/11 11:16

yyki

総合スコア38

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問