Rで箱ひげ図をプロットする

 　Rで箱ひげ図の作成を試みたいと思います。

　

　少しずつ慣れてきた感じはありますが、ネット情報でそのまま導入しても、素人の私にはできないことが多いです。基本を理解していないんだと思います。

　今回は、Qiitaさんのところから参考にしています。

まずは、データを読み込みます（データはでたらめに作成したもの）。

#箱ひげ図　データは縦に並べ、もう一列にSpeciesとして設定する

dat <- read.csv("Book2.csv")

→こんな感じのデータです。

＃ggplot2を立ち上げる

install.packages("ggplot2")

library(ggplot2)

g <- ggplot(Book2)+

  geom_boxplot(aes(x=Species, y=test))+

  geom_jitter(aes(x=Species,y=test,color=Species),width=0.3)

plot(g)

　上記のような箱ひげ図がプロットされました。

　もしカラーをはずしたい場合は、color=Speciesをそのまま削除します。

　もし、灰色カラー（背景）をはずしたい場合は、

theme_classic()+

を加えるようです。すると下記のような図がプロットされます。

　その他にたくさんのおしゃれな図があるようなので、トライしてみたいと思います。

Rで2要因のグラフをプロットする方法

 　2要因のデータからエクセルにて標準偏差や標準誤差を加えた折れ線グラフを作成したとき、どうも2要因のデータ同士が重なって見えずらくなる現象が起きてしまっていました。

　エクセルの手法で改善できるかもしれませんが、Rで簡単にプロットできることがわかりました。

　下記の表のようなデータ（でたらめに作られたもの）を使うとします。

　このデータファイル名をinfecとし、csvファイルに保存します。

　実際のRで作成してみます。

dat <- read.csv("infec.csv") #データ読み込み

＃プロットと平均値や標準偏差を求めるためこの2つを読み込む　

　なければinstall.packagesで対応

library(ggplot2)

library(dplyr)

library(cowplot)

library( Hmisc )

pd = position_dodge( width = 0.15 )

ggplot( infec, aes( x = prepos, y = infection, color = intervention, group = intervention ) ) +

  stat_summary( fun.y = "mean", geom = "line" )+

  stat_summary( fun.y = "mean", geom = "line", size = 1.0, position = pd ) +

  stat_summary( fun.data = mean_cl_normal, geom = "errorbar",

                alpha = 0.9, size = 1.0, width = 0.2, position = pd ) + 

  stat_summary( fun.y = "mean", geom = "point", size = 3.5, position = pd )+

  scale_y_continuous( limits = c( 20, 70 ), breaks =seq(20,70,5), expand = c( 0, 0 ) ) +

  theme_cowplot( 15 ) +

  scale_colour_manual( values = c( "#5B84B1FF", "#FC766AFF" ) ) +

  xlab( "Differences in infection status by town" ) +

  ylab( "Infection" )+ 

  scale_x_discrete(limit=c('before', 'after'))

　グラフの出来上がりです。

　このように、2要因のデータ間で隙間があるので、標準偏差が描かれてもはっきり違いが判ります。

　現時点で、素人の私にとってこれが限界ですが、もう少しグラフの調整ができるように引き続き調べていきたいと思います。

Rを使ってROC曲線を求める方法

　

　ROC曲線を求めるとき、普段使用しているHADには備わっていないため、spssを使用していました。

　spssが嫌いなわけではありませんが、高額なこと、使える場所が限られてしまうことから、やや抵抗感がありました。

　そこで今回、Rを使ったROC曲線を求める方法を探していたところ、統計ERさんのブログでとてもわかりやすく説明されているのを発見しました。

　ただし、この通りに進めてみましたが、dat機能でvariablesが上手く分けられませんでした。もちろん、私の知識の問題がありますが、別の方法で解決できましたので、紹介させて頂きます。現在使用しているRは、「version 4.2.2」です。

　「Book1」のcvsファイルをデータとして使用しています。このデータは、BellCurveさんからダウンロートできます。BellCurveさんのファイル名は「ROC-curve-sample-data」です。

このデータのHeadingには、

検査項目：testA 

疾患の有無：disease　（疾患無しを0、疾患有りを1）

と適当に英語で記載しました。

　現在のRversion 4.2.2機能では、日本語で記載すると上手く読み込めなくなる？ためです。

まず、パッケージをインストールします。

install.packages("Epi")　

library(Epi)

dat = read.csv("Book1.csv", na.strings="") →ここが統計ERさんと異なるところ

統計ERさんの場合→dat <- read.table("ROC-curve-sample-data.txt")

str(dat)

ROC(test=dat$testA, stat=dat$disease, plot="ROC")  #testAが検査の点数で、diseaseが真の分類

ROC(test=dat$testA, stat=dat$disease, plot="sp")

　すると、下記のような図がプロットされます。

　BellCurveさんと同じ結果が得られてますね。

　このEpi パッケージの ROC() 関数の魅力は、AUC面積の算出はもちろんのこと、カットオフ値まで同時に判定してくれるところです。

　spssの場合、Youden Indexなどを求めないと、明確なカットオフ値を得ることができません。そして統計、Rの素人である私にも使いやすい。

　

　その他Rで用いる便利な方法があるか、今後も探していきたいと思います。

※Rのパッケージを論文に引用する必要がある場合、

今回のEpiを用いたとき、

citation("Epi")関数を使います。すると、下記の通り表示されます。

To cite Epi in publications use:

  Bendix Carstensen, Martyn Plummer, Esa Laara, Michael Hills (2022). Epi: A Package for Statistical Analysis in Epidemiology. R

  package version 2.47. URL https://CRAN.R-project.org/package=Epi

If you use Lexis objects/diagrams, please also cite:

  Martyn Plummer, Bendix Carstensen (2011). Lexis: An R Class for Epidemiological Studies with Long-Term Follow-Up. Journal of

  Statistical Software, 38(5), 1-12. URL https://www.jstatsoft.org/v38/i05/.

For use of Lexis objects in multi-state models, please also cite:

  Bendix Carstensen, Martyn Plummer (2011). Using Lexis Objects for Multi-State Models in R. Journal of Statistical Software,

  38(6), 1-18. URL https://www.jstatsoft.org/v38/i06/.

To see these entries in BibTeX format, use 'print(<citation>, bibtex=TRUE)', 'toBibtex(.)', or set

'options(citation.bibtex.max=999)'.

作業療法士という職種について（学校選びについて）

 　個人的な見解ですが、教員と臨床を担っている立場として、作業療法士の数は圧倒的に少ないと思っています。

　

　就職の時期になると、作業療法学生1人に対して6～10施設から就職への依頼がくることもあります。

　作業療法士は、医療・介護領域においても、中心的役割を担う職種です。

　

　つまり、作業療法士は他の職種のことも理解しつつ、リーダー的存在としてチームをひっぱっていくことになります。

　身体機能にかかわる理学療法士以上に、応用的な動作へのアプローチが求められます。食事にかかわる言語聴覚士以上に、嚥下や姿勢、精神面へのアプローチが求められます。

　それだけ作業療法士は必要とされる職業なんですね。

　また作業療法士は、看護師や介護士さんのように夜勤がありません。

　残業も少なく、休日もしっかり確保されています。ゆとりのある生活スタイルを保つという意味では、作業療法士は非常に安定した職業といえるでしょう。

　作業療法士をめざす学校には大きく専門学校、大学に分けられます。

　学校を選ぶ基準はさまざまですが、私がおすすめする見方として以下の3点を挙げます（個人的な意見です）。

1：国家試験合格率

2：教員人数

3：付属病院の有無

1：国家試験合格率

　国家試験合格率は、学生や教員の関係性、そして学校そのものの環境の良さがもろにあらわれると考えることができます。

　一見、偏差値が高い学校でも、なぜか国家試験の合格率は低いことが頻繁にあります。

　つまり、学生や教員が国家試験に対してあまり取り組んでいない可能性があると考えることができます。

　この点において、学校選びの基準として国家試験合格率は重要な見方になるでしょう。

2：教員人数

　当たり前のことですが、教員人数が多いほど、一人当たりの授業数やその他業務に割り当てられる量は減ってきます。

　つまり、教員人数が多いほど、個別で学生にかかわれる時間をもつことができると考えることができます。

　教員は常に精神的に安定した状態のなかで、学生への教育を提供できる環境下にいるということになります。

　学校を選ぶ際は、教員人数を注目することに価値があるといえます。ただし、注意が必要です。教員の数といっても、その教員が作業療法士であるかどうかが重要になってきます。

　作業療法士でない教員は、それに属する学科に深くかかわれないことが多いです。作業療法に関する指導ができない分、実習のサポートも薄くなってしまう可能性があると思います。

3：付属病院の有無

　付属病院があるということは、卒業に必須な実習を身内で受けることができます。

　これは、教員そのものが学生に直接指導することができることに加え、教員と臨床スタッフのやりとりも密に行えることができます。

　担当する教員がいつでも早急に学生をサポートできる点でいえば、付属病院を有する学校を選ぶことは非常に有意義になるでしょう。

　この3つを担保している学校といえば、杏林大学保健学部作業療法学科になります。

　その他付属病院を有している学校には北里大学医療衛生学作業療法学専攻や昭和大学保健医療学部作業療法学科もあります。

　近隣の大学もＨＰで確認しつつ、慎重に学校を選ぶことをおすすめします。

　今回の重要な視点として、大学教員の研究業績は省きましたが、研究業績の高さは教員そのもののレベルを反映すると考えることもできます。教員紹介のところで、研究業績を開示しているので、興味があればぜのぞいてみると良いかと思います。

　作業療法士は心身ともに負担のない職業で、年を重ねても健康な状態でお仕事ができる職業の一つだと思います。

　さまざまな職業と見比べてみて、もし興味があれば、ぜひ作業療法士を目指してみてはいかがでしょうか。

日産軽の100％電気自動車

 　日産の軽が新たに100％電気自動車として来年度に発売されるそうです。

　

　写真はレスポンスさんのホームページより引用

　値段は200万位ですが、国から補助が50万程支給される可能性があり、かなり安価で購入できる品物となっているようです。

　このように、電気自動車も手ごろな値段で購入できる日が近づいてきました。

　私は車に全く興味がありませんので、軽でも何でも使い勝手が良ければ購入してもよいかなとは思っています。

　

　普通の軽自動車よりもさらに大きいみたいですね。

　問題は走行距離ですが、おそらく170㎞程だと言われていますので、遠出をする方は不向きかもしれません。私のように街乗りレベルであれば十分活用できそうです。

　これまでの古いリーフとの走行距離はあまり変わらないかもしれませんが、バッテリー機能がかなり変化しているので、劣化の抑制については期待できそうです。

コロナ対策の一つ換気の重要性について

　コロナ社会において、換気は日常生活に欠かせない用語になってきていると思います。

　まずどこのお店に入っても換気がきちんとなされていないとそれだけでお店のイメージはダウンしてしまうし、入りたくなくなってきます。

　もちろんこれは個人的な見解なので、誰しも同じように解釈すべきであるとは思いません。

　この情報によると、国基準による換気では飛沫をすべて防ぐことができないと言われております。とてもわかりやすい動画ですね。これをみますと、換気が不十分な環境にいると思うだけで何だか気持ち悪くなってきます。

　

　コロナどうこうというより、このような状況下ではとても不潔な感じがして、食事をとること自体かなり抵抗あります。

　理想的な換気を行うためには、最低2ヶ所の窓を開けることが必要だといわれています。1ヵ所でも扇風機やサーキュレーターによってなかの飛沫をきれいに除去することができるんですね（特設サイト新型コロナウイルス）。

　2か所の窓があれば、それぞれの窓にサーキュレーターの向きを外向きと内向きに設置するくことで更なる換気が強化されるようです。

　サーキュレーターは365日私たちの生活に欠かせない存在になりそうなアイテムですね。

アイデンティティ文献検索とKHコーダーによる解釈

　

　先日、アイデンティティと精神的健康の関連について述べました。

　今回、作業療法におけるアイデンティティの国内論文の実態を把握しようと医中誌で検索したところ、該当する論文27件が確認されました。

　先行研究の文献研究を参考にし、27件の抄録をＫＨコーダーを用いて分析してみました。下記の図は共起ネットワークという図になります。

　検索キーワードを軸にしているので、「作業療法」と「アイデンティティ」は四角で囲まれています。

　この図から解釈すると、作業療法士には「職業」をキーワードとした尺度を用いて調査研究がなされていることがわかりました。

　別のグループでは、学生に対するアイデンティティの研究もヒットされているようです。尺度はＰＩ（Professional Identity）を用いて検証がなされています。特に作業療法学生のＰＴは低い状況であることがわかりました。

　また、アンケートの妥当性、有効性を検証しているので、アイデンティティの指標を検討している段階なのかもしれません。

　セルフエフィカシーという用語も確認できました。やはりアイデンティティには、仕事の達成感や自身の可能性を把握するうえで必要な要因なのかもしれません。

　ちょろっと理学療法がでてきています。やはり作業療法を対象としたとき、理学療法もセットで検証されるようです。

　現段階のアイデンティティの指標として、作業療法士には職業アイデンティティ調査が、作業療法学生にはPIが、補足指標としてセルフエフィカシーが挙げられました。

　作業療法士のアイデンティティはどのような状況であるかについては、クリアになっていない、と判断してもよいのでしょうか。質的な研究はヒットされておりませんでした。

　自由回答等で明らかにすることで、より具体的な内容がみえてくるのかもしれません。