29364

掲示板

当ホームページに関するお問い合わせ,講習会,研修会などの問い合わせはこちらへどうぞ.統計相談に関しましては,回答を約束できませんことをご了承ください.

HOME | お知らせ(3/8) | 記事検索 | 携帯用URL | フィード | ヘルプ | 環境設定 | ロケットBBS
Rコマンダーダウンロードについて / つしま
 最近,ダウンロードが集中してリンク切れが起こるようになりました.いまサーバーを変えるなどの対策を考え中ですが,最悪の際には配布できない可能性もあります.
 申し訳ございませんが,ご容赦お願いいたします.
No.401 - 2015/04/17(Fri) 13:05:54
クロスオーバー比較試験 / 坂本
対馬先生 御侍史

対馬先生の”東京図書の書籍”および”R”をいつも使用させていただいております。
理学療法士の坂本と申します。
ランダム化クロスオーバー比較試験の統計解析で行き詰っております。
分かりやすく書かれた書籍(事例も掲載されている)はございますでしょうか?
同じ研究デザインの論文を集めてみたのですが、どれが正しいのか分からなくなっている状態です。
持越し効果の検討も含めて説明されている書籍が助かります。
突然のお願いで大変恐縮ですがアドバイスいただけますと幸いです。
何卒よろしくお願い申し上げます。
No.398 - 2015/02/10(Tue) 10:30:44

Re: クロスオーバー比較試験 / 坂本
問題が一応の解決をみました一報させていただきます。
面倒な質問で申し訳ございませんでした。
No.400 - 2015/03/30(Mon) 08:37:52
分割プロット分散分析について / 渡邉淳一 [ Mail ]
対馬先生

はじめまして。理学療法士の渡邉淳一と申します。この度、改変Rコマンダーの分割プロット分散分析について、ご助言をいただきたくメールさせていただきました。宜しくお願い致します。
先日より、改変Rコマンダーを使用させていただいております。大変重宝しておりますが、分割プロット分散分析の結果がどうしても出力されず、原因が分からないでいます。解析の手順も先生のパワーポイントの説明通りに行っているはずなのですが・・。
データは、a群とb群の2群で、それぞれに別の介入を行い、介入前および二日後、七日後、一ヶ月後の影響をみています。
解析の結果、以下のような出力が成されます。
もし、何か思い当たるような問題があるようであれば、ご助言いただけると幸いです。
ご多忙とは存じますが、宜しくお願い致します。

出力ウインドウ

> aaa <- sqlQuery(channel = 1, select * from [Sheet1$])

> names(aaa) <- make.names(names(aaa))

> names(aaa)[c(1)] <- c(".X1value")

> .Responses <- na.omit(with(aaa, cbind(.X1value, a, b, c, d)))

> anovakun(as.data.frame(.Responses),AsB,range(.Responses[,
+ 1])[2]-range(.Responses[,1])[1]+1,ncol(.Responses)-1,s2n=T,s2r = T,auto = T)

> plot.data.sp( aaa )

> remove(.Responses)

警告
aaa contains non-standard variable names:
1value
These have been changed to:
X1value
エラー:  オブジェクト AsBは存在しません
エラー:  有限な ylimの値が必要です
No.389 - 2014/12/29(Mon) 23:55:15

Re: 分割プロット分散分析について / 対馬
 恐らくこれは,変数名の問題かもしれません.変数名をX1valueではなく,スペースを空けずにX1valueでどうでしょうか?もしくは,最新版を再ダウンロードして下さい.
No.392 - 2015/01/08(Thu) 16:06:23

Re: 分割プロット分散分析について / 渡邉淳一 [ Mail ]
ご助言ありがとうございました。
スペースを空けずにX1valueでも出来ませんでしたので、再ダウンロードして試してみたいと思います。
お忙しい中ありがとうございました。
No.394 - 2015/01/08(Thu) 22:40:45

Re: 分割プロット分散分析について / つしま
 これでダメでしたら,もう一度連絡下さい.
No.395 - 2015/01/09(Fri) 09:17:44

Re: 分割プロット分散分析について / 渡邉
再ダウンロードを行い試してみたところ、結果を出力することができました。
お忙しい中、ご迷惑をおかけしました。

ありがとうございました。
No.397 - 2015/01/13(Tue) 17:33:57
直線の追加とpredict関数について / teruwow [ Mail ]
対馬先生 御机下
改変Rcmdrを利用させていただいております。以前にもダウンロードの件で質問させていただきました杉 と申します。
層別(健常者と脳卒中の転倒者と非転倒者の3層)の散布図を書き、一つの層(健常者)だけの回帰直線と信頼区間・予測区間の直線をかこうとしております。subsetでデータセットを作り、lm関数で回帰分析はできるのですが、その直線を散布図に書くことが出来ません。scatterplotではablineやlinesは利用できないのでしょうか。また、信頼区間はconfintで求められたのですが、予測区間がpredictでうまく求められません。結果をsummaryで表示しようとしても"objectがありません"となってしまいます。Rcmdrdではpredictの結果はobjectの対象にならないのでしょうか?
Rcmdrのダウンロードがうまくいっていないのか、Rcmdrではできないことをやろうとしているのかさえも分かりません。
唐突なうえ、わかりにくい質問で大変失礼と承知いたしておりますが、解決案等をご教授頂ければ幸いです。
ご多忙中、申し訳ありません。
No.384 - 2014/08/08(Fri) 16:35:15

Re: 直線の追加とpredict関数について / つしま
 Rguiの方でやった方が良いかもしれません.コマンダーは対応していません.機能拡張できればいいのですが…
No.396 - 2015/01/09(Fri) 09:26:21
共分散分析について / MM
対馬先生
いつも勉強させていただいております。
コントロール群とケース群(疾病の程度で3群)の合計4群の分散分析を行っております。年齢に偏りがあるので調整するために年齢を共変量として共分散分析しました。結果、平行性は仮定できましたが回帰モデルの傾きは有意でありませんでした。従って共分散分析する意味がないので、ANOVAで分析すればよいと思いますが、各群で年齢に有意差が出ている場合はどのように対処・考察をしたら良いのでしょうか?ANOVAでの有意差は年齢差ではないの?と突っ込まれる可能性があります。ちなみにコントロール群で従属変数と年齢の相関を検定すると、R=0.2で有意でした。勉強不足で申し訳ありませんが何卒ご教授頂きたく、お願い申し上げます。
No.385 - 2014/11/10(Mon) 13:21:06

Re: 共分散分析について / 対馬
 年齢の影響が交絡として有意に影響していないだけ,だと思います
No.393 - 2015/01/08(Thu) 16:08:18
大学院について / 森
弘前大学 対馬先生
初めまして。札幌で理学療法士をしています森といいます。大変恐縮ですが、現在職に就きながら、大学院に進学を目指しております。そこで、大学院の環境についてお聞きしたくご連絡させていただきました。
先生の院生の中には北海道の恵庭や函館の病院に勤めながら、通学している方がいらっしゃいます。
その方々はどのように講義や対馬先生のご指導を受け、学位を取得しているのですか?


お忙しいとは存じますが、何卒ご教授お願い致します。m(_ _)m
No.390 - 2015/01/05(Mon) 19:48:50

Re: 大学院について / つしま
 e-mailにてご連絡下さい.
No.391 - 2015/01/05(Mon) 19:49:43
2つの相関係数の差の検定について / 齋藤
対馬先生

 PTの齋藤と申します。
 対馬先生のHPで紹介されている「2つの相関係数の差の検定」を利用させていただいております。この検定について質問です。
 この検定は、「偏相関係数」についても同様に適応してよいのでしょうか?
 例えば、集団A(100名)における性別と年齢を調整した下肢筋力とADL能力の偏相関係数(r=a)と、集団B(50名)における性別と年齢を調整した下肢筋力とADL能力の偏相関係数(r=b)とすると、aとbの差の検定を行う時に、適応しても良いのでしょうか?
 お手すきの時に御返答いただけましたら幸いです。
No.386 - 2014/11/18(Tue) 06:08:14

Re: 2つの相関係数の差の検定について / つしま
 大丈夫です.
No.387 - 2014/11/18(Tue) 08:48:07

Re: 2つの相関係数の差の検定について / 齋藤
>  大丈夫です.

早速のご回答をいただき、深謝いたします。
活用させていただきます。
No.388 - 2014/11/19(Wed) 04:15:56
ROC曲線の結果解釈について / JY
対馬先生へ
いつも拝見させていただいております。
今回、R2.8.1のROC曲線を用いて解析しているのですが、結果、グラフのAUCにNaNと表示されております。これは対象となる症例が不足しているという解釈でよろしいでしょうか?
また、その際のRコマンダーの出力ウィンドウに、$A、$n…とありますが、$Aは何を示しているのかわかりませんでした。御多忙中申し訳ありませんが、ご教授いただけると幸いです。
No.380 - 2014/08/05(Tue) 09:03:38

Re: ROC曲線の結果解釈について / つしま
http://cran.r-project.org/web/packages/verification/verification.pdf
のP35をご覧下さい.
No.381 - 2014/08/05(Tue) 09:46:57

Re: ROC曲線の結果解釈について / つしま
Nanは症例数に関係ないようです.なぜ出るかは不明ですが,面積は代わりに $A として出力されます.
No.382 - 2014/08/05(Tue) 09:51:01

Re: ROC曲線の結果解釈について / JY
早速のお返事ありがとうございました!
No.383 - 2014/08/05(Tue) 12:28:28
kappa値について / はやし
いつも使用させていただいております。
今回は8人の解答者に3つの設問について2択で答えてもらい一致率を見ることを試みましたが、設問の答えを「0または1」というエクセルデータ「1または2」のエクセルデータでκ値が異なってしまいました。何の数字で置き換えるかの違いだけで、答え自体は一緒のはずですが・・・。(一方で1の解答の場合は一方では0、一方で2の解答の場合は一方では解答1です)使用する数字によって答えが変わることってあるのでしょうか?
わかりづらい文章ですみませんが、対策等ご教授いただけましたら幸いです。
No.377 - 2014/07/08(Tue) 13:52:00

Re: kappa値について / つしま
 おそらく,Kappa(Cohen)を参照していると思われますが,これには入力規則の問題があります.
 まず,データは,如何なる時も「1」から始めないとならない点です.0から入れますと,0=無,と判断して計算します(つまり0と1だけだと,カテゴリーは1のみと判断します).
 また,2,3の2カテゴリーにしますと,1のカテゴリーの一致=0と判断して,計算します.ですので,ランクデータにして入力する必要があります.
 次に,計算の問題です.私が想定しているのは,Fleissのκ係数です.計算は,千野先生(http://www.aichi-gakuin.ac.jp/~chino/anova/chapter3/sec3-1.html#Fleiss71)のページをご参照下さい.この計算を,Rコンコルドパッケージのcohen.kappaでは(つまり私の改変Rコマンダー),Kappa(Siegel)で出力しています.Kappa(Siegel)では,上記の問題はありません(つまり0から始めると計算が違うなど).
 この点については, http://www.hs.hirosaki-u.ac.jp/~pteiki/research/stat/S/kappa/ のページもご参照下さい.
No.378 - 2014/07/08(Tue) 23:35:39

Re: kappa値について / はやし
早々のお返事ありがとうございました。
今後ともどうぞよろしくお願い申し上げます。
No.379 - 2014/07/10(Thu) 22:48:36
ICCについて / しまだあさみ [ Mail ]
理学療法士のしまだと申します。改変Rコマンダーを使用させていただいています。
改変Rコマンダーを使用してICCを出したいのですが、エラーになってしまいます。どんな問題点が考えられるのでしょうか‥。初歩的な質問で大変申し訳ありませんがアドバイスいただければと思います。
No.371 - 2014/03/03(Mon) 19:12:21

Re: ICCについて / つしま
 いろいろなエラーがあるので何とも言えません.
 詳しい情報が必要です.
No.376 - 2014/07/06(Sun) 15:45:47
以下のフォームに記事No.と投稿時のパスワードを入力すれば
投稿後に記事の編集や削除が行えます。
100/100件 [ ページ : << 1 2 3 4 5 6 7 8 9 10 >> ]

- HOME - お知らせ(3/8) - 記事検索 - 携帯用URL - フィード - ヘルプ - メール - 環境設定 -

Rocket Board Type-LS (Free) Rocket BBS