初心者のための R および RjpWiki に関する質問コーナー

以下の 4 つの「●項目」のどれかをクリック!!


新規投稿欄: --- その前に,「投稿における注意事項」を読んでから!
このページの目次 --- 質問への回答・コメントの参照
最新のスレッド --- 最も最近に書き込まれた「親」記事(最新のコメントが付いた記事ではありません)
Q&A (初級者コース) 過去の目次 --- 今までの 16 個の書庫の参照


●● 目次 ●●

参照は,個々の項目をクリック(質問への回答・コメントの参照も,個々の項目をクリック)


まずは,「投稿における注意事項」を読んでください。
これに従わない場合には回答が得られにくい場合があります。

●● 新規投稿 ●●  




参考書籍

ムーミン? (2018-04-21 (土) 18:25:20)

よろしくお願いします
エクセルからのインポートはできかしたが、そのデーターからグラフを作成しょうとWEBをみて回りましたがさっぱりわかりません。(途中が省かれたり解説なかったり、例がインポートしたものではなかったりなど)”Rによるデーターサイエンス”の本図書館から借りてみましたが上級者向けで理解できませんでした。ホントの初心者です、エクセルからインポートしたデーターをグラフ化することが詳しく解説されている書籍を教えてください。
めげそうです。

  • インポートしたデータであろうがなかろうが,データフレームとして格納されたデータをグラフ化するという点では同じでしょう。
    まずは,どのようなグラフを作りたいのか? -- 河童の屁は,河童にあらず,屁である。? 2018-04-21 (土) 21:54:40
  • http://cse.naro.affrc.go.jp/takezawa/r-tips/r/48.html 参考にどうぞ -- 2018-04-24 (火) 11:15:50
  • 金先生の「Rによるデーターサイエンス」は初心者向けです。初心者向けなので、上級者向けの詳細をばっさり削られて逆に理解できないという状態でしょうか。まずは、書籍よりもRに付属するマニュアルでRグラフィックの基礎の基礎を理解されてはいかがでしょうか。Rをインストールしていればマニュアルを参照できますが、オンラインで閲覧したければ、ここを参照して下さい。具体的には、An Introduction to Rの12 Graphical proceduresになります。日本語訳は間瀬先生がこつこつ作成されています。このWikiシステムのどこかにリンクがあったと思います。なお、視覚化する対象データがもともとエクセルのデータであろうが、SPSSのデータであろうが、視覚化作業には全く関係がないです。 -- S? 2018-04-24 (火) 11:20:52
  • http://cse.naro.affrc.go.jp/takezawa/r-tips/r/48.html 参考にどうぞ -- 2018-04-24 (火) 14:56:56
  • 回答ありがとうございます -- ムーミン? 2018-04-24 (火) 18:29:02

R言語へのエクセルからのデーターのインポート

ムーミン? (2018-04-18 (水) 20:14:01)

どこに書き込んだらよいかわからなかったので新規投稿に書きました
よろしくお願いします。
手順は
1)data4=read.table("clipboard",header=TRUE)を記述、
2)データーをコーピー
3) 1)をエンターキーを押しました。
又は
1)、2)を逆にし実行

下記が結果です。何も起こりません


R version 3.4.4 (2018-03-15) -- "Someone to Lean On"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)

R は、自由なソフトウェアであり、「完全に無保証」です。
一定の条件に従えば、自由にこれを再配布することができます。
配布条件の詳細に関しては、'license()' あるいは 'licence()' と入力してください。

R は多くの貢献者による共同プロジェクトです。
詳しくは 'contributors()' と入力してください。
また、R や R のパッケージを出版物で引用する際の形式については
'citation()' と入力してください。

'demo()' と入力すればデモをみることができます。
'help()' とすればオンラインヘルプが出ます。
'help.start()' で HTML ブラウザによるヘルプがみられます。
'q()' と入力すれば R を終了します。

data4=read.table("clipboard",header=TRUE)

  • その操作の後,単に data4 と 5 文字入力した後に,リターンキーを押したらどうなりますか??? -- 河童の屁は,河童にあらず,屁である。? 2018-04-18 (水) 22:04:59
  • ありがとうございます。できました。 -- ムーミン? 2018-04-19 (木) 18:22:36

症例対照研究のロジスティック回帰について

AAA? (2018-04-17 (火) 22:42:38)

お世話になります。
glm関数で症例対照研究のロジスティック回帰を実行しようとしたのですが、症例対照研究は症例群と対照群のサンプル数を任意に設定できるので、生存率やオッズの算出に意味はないという資料がありました。ロジスティック回帰式の左辺はlog(p/(1-p))で対数オッズを意味していますので、そもそも回帰することができないということになるでしょうか。
glmに目的変数と説明変数を与えると計算結果自体は出力されるのですが、根本的に前提が間違っているような気がしています。

宜しくお願いします。

  • 症例群と対照群を併せて分析しようとしているのでしょうか?ならば,サンプルサイズの差が結果に影響するのは当たり前でしょう。
    普通は,別々に分析して結果(パラメータ,オッズ比など)の違いを比較・検討するんじゃないですか? -- 河童の屁は,河童にあらず,屁である。? 2018-04-18 (水) 09:19:07

R言語へのエクセルからのデーターのインポート

ムーミン? (2018-04-17 (火) 18:19:45)

よろしくお願いします
下記データーをRコンソール、にクリップボードを介してインポートする方法を教えてください。

R言語のインストールは下記のページを参考にして行いました
https://qiita.com/FukuharaYohei/items/8e0ddd0af11132031355

インポートは下記ページを参考にしましたがうまくいきません
http://monge.tec.fukuoka-u.ac.jp/r_analysis/input_clipbord00.html

月 電気量 電気料金
1月 306  8184
2月 335  9065
3月 341  9311

R言語のバージョンは3.4.4

ウィンドウズ10でエクセルのバージョンは1803です

  • いやー,ビックリしたわ。初級Q&A15へ投稿したんだねえ。
    あなたが,どういう風にしたのか,その結果どのようなエラーメッセージが出たのか,コンソールの出力をそのままコピーして,ここに,ペーストしなさい!!! -- 河童の屁は,河童にあらず,屁である。? 2018-04-17 (火) 20:45:31
  • よろしくお願いします。 -- ムーミン? 2018-04-18 (水) 20:02:06
  • R version 3.4.4 (2018-03-15) -- "Someone to Lean On" -- ムーミン? 2018-04-18 (水) 20:03:18

2つのデータをlmで回帰したときに傾き切片が有意差を持つかどうか

kagakuma? (2018-04-13 (金) 12:00:31)

両対数を取ると明らかに線形関数に乗るデータがあります。
このデータをlm()で線形回帰したところR2が0.95でフィッティングしました。
このデータ群は2つあり、線形回帰したのですが、いずれも違う直線に乗りました。 この結果から傾きと切片が2組でますが、両者の傾きと切片に有意差があるかどうかを確かめる方法はございませんでしょうか。
Rの機能で閉じていると助かりますが、手計算でもできるならやろうと思います。
単純には両者を引き算してそれの傾きと切片を求めるという方法も考えましたが、正当性がわかりません。 よろしくおねがいします。

  • 共分散分析というものです。http://www.snap-tck.com/room04/c01/stat/stat99/stat0204.pdf に重回帰分析との関連が詳しく書かれています。lm でできます。 -- 河童の屁は,河童にあらず,屁である。? 2018-04-13 (金) 14:26:25
  • 河童様:ご返事ありがとうございます。lmでできる、という言葉の意味が取れなかったので重ねての質問になりますが、二群をクラス分けして一本のデータにしたものをanova(aov(value~factor))に掛けるのと等価の処理でしょうか。共分散分析は勉強したはずなのに思いつきませんでした。精進します。 -- kagakuma? 2018-04-13 (金) 17:40:43
  • データを挙げて,結果も示すという分析例があるのですから,まずは,その例題を lm を使ってやってみれば良いと思います。
    ただ,それをここで書いて見せてあげても仕方ないでしょうということ。 -- 河童の屁は,河童にあらず,屁である。? 2018-04-13 (金) 19:09:49

入門機械学習による異常検知p137

ksaki? (2018-04-12 (木) 14:42:13)

お世話になります。

井手剛先生著の
入門 機械学習による異常検知のp137の実行例5.1を真似て、以下のコマンドを入力しましたが、エラーが出ます。
何が間違っているか教えてください。
以下、入力したコマンド

> data <- read.table("R_r1.csv",header=T)
> data
                                                  a.b.c.d.e.f.g.h.i.j.k
1         P2201A1,-4.78,11,1.44,-0.48,3.83,-10.05,15.79,4.31,1.44,18.66
2     P2201A2,-3.35,-21.53,-11,8.13,-6.22,18.66,-2.87,-0.96,2.39,-11.48
3            P2201A3,2.88,3.63,4.69,0.61,1.21,1.36,-1.67,1.06,1.06,1.67
4       P2201A4,-0.48,3.35,-5.26,-6.22,-3.35,-4.78,0.48,-0.48,1.91,3.35
5    P2201H1,-14.35,-5.26,0,-11.96,-7.66,14.35,-11.48,-5.26,14.35,-0.48
6  P2201H2,12.92,13.87,0.96,13.4,-0.96,1.91,-12.44,-13.87,-10.05,-14.35
7           P2201H3,7.18,3.83,3.35,0,-4.31,-8.61,-12.44,-3.35,1.44,5.74
8     P2201H4,-7.18,-16.27,-16.27,-5.26,2.39,8.61,12.92,12.92,6.7,-0.48
9        P2201V1,-6.7,-3.83,-12.92,2.39,2.39,2.87,12.44,2.87,1.91,-8.13
10   P2201V2,-3.03,-18.17,22.71,-27.25,27.25,-15.14,0,-10.6,-18.17,21.2
11      P2201V3,-4.31,-0.48,1.91,-1.44,-5.26,1.44,3.35,-8.61,-6.7,-2.87
12        P2201V4,0.48,2.39,1.91,3.35,-1.44,-7.18,-5.74,-1.91,7.66,7.66
13       P2202A1,-1.44,-0.96,3.35,1.91,-0.48,0.96,-3.35,-2.39,0.48,-6.7
14         P2202A2,-3.83,0,1.44,-6.22,-3.35,0.96,-1.44,0.48,-0.48,-1.91
15      P2202A3,0.61,-1.36,-0.15,0.61,0.76,-2.57,-2.12,-2.57,-0.91,0.76
16        P2202A4,0.45,2.27,-0.45,1.21,0.15,-0.15,0.76,-1.06,-1.21,1.06
17       P2202H1,5.26,-4.31,-0.48,-1.44,-3.83,0.48,0.48,5.74,-0.48,0.48
18         P2202H2,-0.96,-1.91,-0.96,-5.74,0,1.91,-5.74,1.44,-1.91,-6.7
19       P2202H3,-1.67,0.45,-0.45,2.57,2.27,-1.21,-0.61,-0.3,-0.15,0.45
> cc <- c("a","b","c","d","e","f","g","h","i","j","k")
> mask <- is.element(colnames(data),cc)
> Xc <- t(scale(data[,mask]))
> colnames(Xc) <- t(data[,"a"])
Error in `[.data.frame`(data, , "a") :  未定義の列が選ばれました

コマンド、ここまで
〇環境設定

> sessionInfo()
R version 3.3.3 (2017-03-06)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

locale:
[1] LC_COLLATE=Japanese_Japan.932  LC_CTYPE=Japanese_Japan.932   
[3] LC_MONETARY=Japanese_Japan.932 LC_NUMERIC=C                  
[5] LC_TIME=Japanese_Japan.932    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] MASS_7.3-45

loaded via a namespace (and not attached):
[1] tools_3.3.3
  • 入力ファイル "R_r1.csv" は自分で入力したのでしょうか。data の次の出力を見ると,19 行のデータは,"a.b.c.d.e.f.g.h.i.j.k" という名前を持つ文字型の1変数として入力されていますよ。その原因のひとつは,変数名の行が "," ではなく "." で区切られていること。もう一つは,もし read.table を使うのなら,sep="," を加えるか,あるいは read.table ではなく read.csv を使う(ファイル名からいうと,read.csv を使う前提なのだろうけど)。 -- 河童の屁は,河童にあらず,屁である。? 2018-04-12 (木) 18:13:01
    > data = read.csv("R_r1.csv") # read.csv は header = TRUE はデフォルトなので,指定不要
    > data
             a      b      c      d      e     f      g      h      i      j      k
    1  P2201A1  -4.78  11.00   1.44  -0.48  3.83 -10.05  15.79   4.31   1.44  18.66
    2  P2201A2  -3.35 -21.53 -11.00   8.13 -6.22  18.66  -2.87  -0.96   2.39 -11.48
       中略
    18 P2202H2  -0.96  -1.91  -0.96  -5.74  0.00   1.91  -5.74   1.44  -1.91  -6.70
    19 P2202H3  -1.67   0.45  -0.45   2.57  2.27  -1.21  -0.61  -0.30  -0.15   0.45
    > cc = c("b","c","d","e","f","g","h","i","j","k") # 以下の意図はよく分からないけど,"a" は外さないとエラーが出るよ
    > # cc = letters[2:11] とも書けるけど
    > mask = is.element(colnames(data),cc)
    > Xc = t(scale(data[,mask]))
    > # Xc = t(scale(data[, 2:11])) で十分。mask や cc なんか使わなくても良いし
    > colnames(Xc) = t(data[,"a"]) # t は不要
  • 河童様、教えて頂き、感謝の気持ちでいっぱいです。
    "R_r1.csv"は自分で入力しました。テキストで開くと","になっています。が、"."表示になっているのを理解しました。 -- ksaki? 2018-04-12 (木) 18:40:19
  • 「テキストで開くと","になっています」
    そうですか。もしかして,全角のカンマ(あるいは,半角だけど ascii 文字でないカンマ)で入力されていませんか? -- 河童の屁は,河童にあらず,屁である。? 2018-04-12 (木) 18:44:09
  • エクセルに張り付けた測定記録に追加で自分で各行にabc..を追加して作ったcsvなので、何故上段の1行目と2行目で"."と","なっているのか頭の整理ができません。試して後ほどご報告いたします。 -- ksaki? 2018-04-12 (木) 23:46:33


トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Google
WWW を検索 OKADAJP.ORG を検索
Last-modified: 2018-04-24 (火) 18:29:02 (19h)