群馬大学 | 医学部 | サイトトップ | 医学情報処理演習
第12回に出てきた関数と文の主なものはこちらを参照されたい。
vcdパッケージに含まれているArthritisというデータは，慢性関節リウマチ(RA)患者84人を対象とした二重盲検による臨床試験の結果である。患者をプラセボ群(Placebo)と新しい治療群(Treated)にランダムに割り付け（変数名Treatment），症状に改善がみられたかどうか(変数名Improved）を，改善なし(None)，多少改善(Some)，顕著に改善(Marked)の3段階で評価している。
改善がみられたかどうかを改善なしと改善ありの2段階に再区分して，新しい治療の有効性を，改善ありのリスク比によって示し，さらに治療の有無と改善の有無の関連性を属性相関係数によって評価せよ。
学籍番号・氏名とともに，下のフォームと解釈文を穴埋めして送信せよ。
The data Arthritis included in the vcd package is the result of double blind randomized controlled trial about the 84 rheumatoid arthritis (RA) patients. Randomly assign the patients into 2 treatment groups (variable name: Treatment) of Placebo or Treated, then evaluate the extent of improvement (variable name: Improved) by three levels of None, Some, or Marked.
Relevel the improvement into 2 levels of None and Exist, then calculate the risk ratio as the effectiveness of the new treatment compared to placebo, and evaluate the relationship between the treatment and improvement by association coefficients.
Please write the registry number and name, fill the boxes by adequate characters.
コードは以下の通り。
(The code is shown below.)
library(vcd)
TAB <- xtabs(~+
, data=Arthritis)
TAB2 <- cbind(TAB[,1],TAB[,2]+TAB[,3])
colnames(TAB2) <- c("None","Exist")
TAB2
TreatImproved <- TAB2[2,2]
PlaceboImproved <- TAB2[1,2]
TreatAll <- rowSums(TAB2)[2]
PlaceboAll <- rowSums(TAB2)[1]
library(fmsb)
riskratio(TreatImproved, PlaceboImproved, TreatAll, PlaceboAll)
(TAB2) # calculate association coefficients
以上の結果から，プラセボ群に比べて新治療群では約2.1倍の改善が見られたといえる。この改善は有意水準5％で統計的に有意である。また，属性相関係数の値はすべて0.3から0.4の間であり，治療と改善の間には
（1:弱いが有意な関連がある，2:強くて有意な関連がある，3:関連が無い）といえる。
Based on the above results, the improvement in new treatment group was about 2.1 times higher than placebo group, which was statistically significant at 5% level. Because the all association coefficients ranged 0.3 to 0.4, we can state that there is (1: weak but significant, 2: strong and significant, 3: not significant) relationship between treatment and improvement.
| 項目 | 入力欄 |
|---|---|
| Treatment （＊多かった誤答は，treatment，Treatedなど） |
| Improved （＊多かった誤答は，improved，Improveなど） |
| assocstats （＊誤答は1人だけ，assoctatsというタイプミスがありました） |
| 1 （＊2と判定した人もかなりいましたが，CramerのVが0.3〜0.4では「強い関連」とはいえません。しかし「関連が無い」帰無仮説は有意水準5%で棄却されるので，統計的に有意な関連があるとはいえます。そこで解釈は1となります） |
変数名は正確に書かないと動作しないので注意してください。今回は回答を要する質問・コメントはありませんでした。