群馬大学 | 医学部 | サイトトップ | 医学情報処理演習
前回も参照した通り，http://phi.med.gunma-u.ac.jp/medstat/sample02.txtは，ソロモン諸島の首都のとある学校で実施した健診結果を，昨年の実習で入力してもらった後にエラーを訂正したタブ区切りテキスト形式データである。
このファイルをインターネットから直接Rのデータフレームxに読み込み，婚姻状態（変数名はMARITAL）別に，Body Mass Index（変数名はBMI）の記述統計量（サンプルサイズ，平均値，標準偏差，第１四分位，中央値，第３四分位）を求めるコードは以下のように書ける。ただし，わざとBoxで穴をあけてある。
descstat <- function(X) {
XX <- X[!is.na(X)] # removing missing values
Q5S <-(XX)
res <- c((XX),
(XX), sd(XX), Q5S[2], Q5S[3], Q5S[4])
names(res) <- c("N","Mean","SD","Q1","Median","Q3")
return(res) }
x <-("http://phi.med.gunma-u.ac.jp/medstat/sample02.txt")
tapply(x$BMI,x$MARITAL,descstat)
それぞれのBoxに入る正しい関数名を解答せよ。また，求めた記述統計量からいえることを文章で記せ。さらに，もしあれば感想・要望・質問などとともに下のフォームから送信せよ。学籍番号と氏名の入力を忘れないこと。
(The program listed above is to calculate descriptive statistics (sample size, mean, unbiased standard deviation, 1st and 3rd quartiles, and median) of body mass index (BMI) for married and single subjects, separately. Please fill appropriate functions in the boxes. In addition, please explain what the results mean.)