異なる日付の臨床分離株から特定された多数の感染のデータフレームがあります。
これまでのところ、作業を開始したい形にデータを編成しました。私は、レポートの記述統計のために一連の表を準備しようとしています。
私はftable
を使用しており、以下を取得します。
onset.types <- ftable(SAB$Onset,SAB$MRSA.Type,year(SAB$Collection.Date))
2005 2006 2007 2008 2009 2010
Community 454 472 512 499 525 512
AUS-2/3-like 28 23 27 29 32 38
EMRSA-15-like 9 4 4 9 8 8
nmMRSA 40 47 53 39 64 60
Other mMRSA 1 3 3 11 5 9
unclassified MRSA 0 2 0 0 1 1
Hospital 163 163 156 164 149 165
AUS-2/3-like 31 33 27 31 29 28
EMRSA-15-like 3 8 5 9 4 3
nmMRSA 10 9 13 17 13 12
Other mMRSA 5 1 6 2 3 10
unclassified MRSA 2 0 1 0 0 0
2つの質問:
1:限界合計をどのように計算しますか
2:パーセンテージを計算し、限界合計で再度カウントする簡単な方法はありますか
私はエピツールを試しましたが、思ったほど便利ではありませんでした。
どうもありがとう。
マージンを追加するには、addmargins()
を使用します
_addmargins(table(state.division, state.region))
state.region
state.division Northeast South North Central West Sum
New England 6 0 0 0 6
Middle Atlantic 3 0 0 0 3
South Atlantic 0 8 0 0 8
East South Central 0 4 0 0 4
West South Central 0 4 0 0 4
East North Central 0 0 5 0 5
West North Central 0 0 7 0 7
Mountain 0 0 0 8 8
Pacific 0 0 0 5 5
Sum 9 16 12 13 50
_
パーセンテージを計算するには、prop.table()
を使用します
_prop.table(table(state.division, state.region))
state.region
state.division Northeast South North Central West
New England 0.12 0.00 0.00 0.00
Middle Atlantic 0.06 0.00 0.00 0.00
South Atlantic 0.00 0.16 0.00 0.00
East South Central 0.00 0.08 0.00 0.00
West South Central 0.00 0.08 0.00 0.00
East North Central 0.00 0.00 0.10 0.00
West North Central 0.00 0.00 0.14 0.00
Mountain 0.00 0.00 0.00 0.16
Pacific 0.00 0.00 0.00 0.10
_
また、調査パッケージのsvytable()
でも機能します。
addmargins(svytable(formula = ~x1+x2, design = df.w))