トップへ ｜ テーマトップへ ｜ テーマランキング一覧へ
| タイトル | 日 時 |
|---|---|
成人病３項目グラフ【７年半】：2008/7/24
|
2008/07/24 07:44 |
体脂肪率の管理図【約１ヶ月】：20080723
|
2008/07/23 03:15 |
体重の管理図【約1ヶ月間】：20080723
|
2008/07/23 02:55 |
血圧・脈拍グラフ【2年2ヵ月間】：20080722
|
2008/07/22 06:00 |
朝顔日記【花数グラフNo.5】：20080713
|
2008/07/13 06:37 |
室温と体脂肪率【散布相関図】：20080711
|
2008/07/11 08:37 |
気温と体脂肪率【時系列的相関係数のグラフ】：20080711
|
2008/07/11 08:18 |
朝顔日記【花数グラフNo.3】：20080706
|
2008/07/06 13:48 |
体脂肪率測定について【重回帰分析】：20080703
３５日間の体重、体脂肪率(１日に３回立て続けに測定しその平均をとる)と気温、湿度(安価なデジタル式で測定)のデータから、ゆっくり変化する成分を取り除いて、日ごとの相互の関係を分析した上で、外れ値を取り除いて、体脂肪率の値を推定する式を立てると ...続きを見る |
2008/07/03 05:29 |
朝顔日記【花数グラフ】暫定：20080629
|
2008/06/29 05:54 |
体脂肪率のXバーR管理図【１ヶ月】：20080629
|
2008/06/29 03:55 |
血圧のXバーR管理図【１ヶ月間】：20080629
|
2008/06/29 02:43 |
体重グラフ【28ヵ月間】：20080624
|
2008/06/24 08:58 |
気温と体脂肪率【相関散布図】：20080622
|
2008/06/22 13:14 |
血圧・脈拍【25ヶ月間グラフ】：20080620
|
2008/06/20 09:08 |
血圧・脈拍のばらつき【2年分グラフ】：20080609
|
2008/06/09 01:19 |
体重の変化と脈拍【散布図】：20080605
|
2008/06/05 07:25 |
体重のXバーR管理図【１ヶ月間】：20080601
|
2008/06/01 08:11 |
最低血圧XバーR管理図【1ヵ月】：20080522
|
2008/05/22 08:15 |
血圧・脈拍、ここ２年間【折れ線グラフ】：20080521
|
2008/05/21 02:43 |
血圧・脈拍、連続１時間測定生データ【フリー扱い】：20080514
2008/5/11(日)の、昼食(11:45AM終わり)後、さらに昼寝(13:25PM終わり)の後で、１時間にわたり約１分おきに６０回連続で、血圧と脈拍を測定した結果の生データです。学生さんのレポートの補助的な役割を果たせたら、と思うので、敢えてUPいたします。 静かに座っている状況下でも、この程度は変化があるもののようです。なおTVはつけっぱなしでした。TVの内容に血圧・脈拍が、かなり影響を受けたか、と、きかれれば、答えはYesです。もっと具体的に書けないのは、記録補助員がいなかったのでやむ... ...続きを見る |
2008/05/14 09:36 |
１時間血圧・脈拍連続測定、分析【残差相関図、検定】：20080513
|
2008/05/13 03:25 |
血圧・脈拍の連続測定、１時間【折れ線グラフ】：20080512
|
2008/05/12 06:20 |
体脂肪率のXバーR管理図【１ヶ月】：20080509
|
2008/05/09 13:50 |
体重のXバーR管理図【１ヶ月】：20080509
|
2008/05/09 13:47 |
最低血圧のXバーR管理図【1ヵ月】：20080507
|
2008/05/07 07:39 |
6万アクセス突破【御礼】：20080505
アクセス状況の表示で、2008/5/3から5/4にかけて、拙ブログ開始以来の累計アクセスが6万を突破いたしました。 日頃から、定期的にアクセスしてくださっている皆様に御礼申し上げます。 不定期にアクセスしてくださる方々にも感謝いたします。 なお5万アクセス突破は昨年の12/18でした。 ...続きを見る |
2008/05/05 03:46 |
起床時血圧連続測定のばらつき【箱ヒゲ図】：20080504
|
2008/05/04 06:49 |
血圧・脈拍グラフ【１年11ヶ月間】：20080421
|
2008/04/21 03:59 |
最低血圧のXバーR管理図【１ヶ月】：20080419
|
2008/04/19 08:00 |
体脂肪率のXバーR管理図【１ヶ月】：20080419
|
2008/04/19 07:58 |
体重のXバーR管理図【短期】：20080417
|
2008/04/17 08:19 |
体脂肪率ＸバーＲ管理図【短期(１ヶ月)】：20080412
|
2008/04/12 08:31 |
最低血圧XバーR管理図【短期(１ヶ月)】：20080412
|
2008/04/12 08:05 |
体重のXバーR管理図【短期グラフ】：20080408
|
2008/04/08 08:17 |
HbA1c、γ-GTP、中性脂肪(TG)【長期グラフ(７年間)】：20080406
|
2008/04/06 04:16 |
体重グラフ【長期(9か月)】：20080405
|
2008/04/05 08:45 |
体重のXバーR管理図【グラフ】：20080331
|
2008/03/31 08:13 |
１日の歩数と体重変化【グラフ】：20080331
|
2008/03/31 02:57 |
血圧・脈拍、連続測定の生データ：20080325
血圧・脈拍、連続測定の生データ。 フリーのデータとして公開します。 2008/3/23(日)。 7:45AM〜8:45AM。 1回/分。 ±10秒程度の誤差あり。 食事7:20AM頃終わり。 薬は食事直後。 ---------------------- 最高血圧,最低血圧,脈拍 116,74,88 116,77,88 109,76,91 114,78,89 114,78,90 109,77,90 105,78,89 104,80,91 105,75,90 1... ...続きを見る |
2008/03/25 00:36 |
血圧連続測定時のばらつき【ヒストグラム】：20080324
|
2008/03/24 08:40 |
血圧の連続測定【折れ線グラフ】１時間：20080323
|
2008/03/23 12:53 |
血圧・脈拍の22ヶ月間【長期グラフ】：20080321
|
2008/03/21 04:04 |
体重ＸバーＲ管理図【短期グラフ】
|
2008/03/20 02:16 |
血圧の毎朝測定時のばらつき【箱ヒゲ図】：20080313
|
2008/03/13 07:45 |
体重のＸバーＲ管理図【短期グラフ】：20080307
|
2008/03/07 09:46 |
体重グラフ【長期】：20080306
|
2008/03/06 05:03 |
体重のＸバーＲ管理図【グラフ】：20080229
|
2008/02/29 08:40 |
体脂肪・内臓脂肪の変動【グラフ】：20080227
|
2008/02/27 04:22 |
連続測定時の血圧・脈拍のばらつき【頻度図】：20080225
|
2008/02/25 15:38 |
血圧脈拍１時間連続測定【相関図と検定】：20080225
『等倍』でご覧下さい。 一昨日の夕食後に１分おきに１時間かけて測定した血圧・脈拍のデータをプロットしたものです。最高血圧と最低血圧の値が、ほかから飛び離れていた１回づつ、２回分をとりのぞいてありますのでN=58です。 この中で関連性が特に高いのは最低血圧と脈拍です。以前の同様な分析でも類似の結果が得られていますから、これは、かなり本当なのではないかと思われます。無論、これは、私個人についてのデータなので、学生さんは、自分のデータを採って比べることをお勧めします。 ...続きを見る |
2008/02/25 04:20 |
血圧脈拍１時間連続測定【生データ】
今回公表した１時間にわたる血圧・脈拍の測定値です。フリーのデータとして公表します。学生さんのレポートなどの参考になれば幸いです。なお、測定間隔は『ほぼ』１分です。最大、10秒程度の誤差があるものと思ってください。 ...続きを見る |
2008/02/24 03:36 |
血圧脈拍１時間連続測定【グラフ】
08/2/23の夕食後１時間〜２時間にかけての１時間に亘る血圧・脈拍の測定値のグラフです。血圧も脈拍も徐々に上昇していることがわかります。私は、食後に眠たくなる傾向があるのですが、血圧が下がっているためではないかと思われます。 今度のチャンスには食事の直後から計ってみたいと思います。なお、今回も、生データと各々の値の相関図も公表の予定です。 ...続きを見る |
2008/02/24 03:34 |
体重管理図【ＸバーＲ】：20080222
『等倍』でご覧下さい。 体重のＸバーＲ管理図です。 極めて異常だった2/9,2/10のデータは、除去しました。 70.1Kg付近に安定してきた、と13日に書きましたが、17日頃から一段と体重が減少したことが分かります。今後もリバウンドせぬよう、食べ物と運動に注意したいと思います。 ...続きを見る |
2008/02/22 08:45 |
血圧・脈拍21ヶ月間【グラフ】：20080219
『等倍』でご覧下さい。 ほぼ、21ヵ月間に亘る血圧の変動です。毎日、起床時、７回連続で測定したものの、平均値と、最も高かった値と低かった値の幅(範囲：レンジ【R】)です。2006/9/27に、測定姿勢と血圧計を変更していますので、厳密には、その前後の比較はできないのですが……。一時低下した最高血圧が、再度上がってきていますが、部屋の気温は下げていないので、別な理由を考えなくてはならないでしょう。 下の段は、上の、毎日の血圧の変動のグラフの、７日移動平均をとって、動きを滑らかにしてみたものです... ...続きを見る |
2008/02/19 05:29 |
体重XバーR管理図【グラフ】：20080213
『等倍』でご覧下さい。 体重のＸバーＲ管理図です。 2/9,2/10の異常は、水の飲み過ぎと下剤の服用忘れ(下剤は薬の副作用を打ち消すために飲んでいる)のためです。 70.1Kg付近に安定してきているのがわかります。リバウンドせぬよう、運動にも努めたいと思います。 ...続きを見る |
2008/02/13 17:30 |
起床時血圧７回連続測定【グラフ、箱ヒゲ図】：20080213
『等倍』でご覧下さい。 2007/12/5〜2008/2/13にわたって、毎日起床時、７回連続で測定した、最高・最低血圧と脈拍が、毎日の平均値からどのくらいずれていたかを示した「箱ヒゲ図」です。 データの量がすくないためか、これまでみられた、１回目(V1,V8,V15)の値のずれが、特に大きいという傾向が見られません。 この傾向が続くなら、睡眠時の室温を温かくしたことの影響か？などと言えるかもしれませんが、本当のところ、よくわかりません。 また２ヶ月くらい経ったら、再度、データを集計して... ...続きを見る |
2008/02/13 04:58 |
今までのアクセス数【グラフ】：20080211
詳しくは図内コメントをご参照ください。 ...続きを見る |
2008/02/11 04:18 |
血圧・脈拍の連続測定でのバラツキ【グラフ】
『等倍』でご覧下さい。 2008/1/6の早朝、１時間の血圧の変動を計った時の、データのばらつきが、どれくらいあるかを示した図【頻度図、ヒストグラム】です。 最高血圧で±10、最低血圧で±8、脈拍で±7、位のずれは覚悟しないといけないようですね。 ...続きを見る |
2008/01/28 03:46 |
血圧の変動【約20ヵ月間】：20080118
『等倍』でご覧ください ほぼ、20ヵ月間に亘る血圧の変動です。毎日、起床時、７回連続で測定したものの、平均値と、最も高かった値と低かった値の幅(範囲：レンジ【R】)です。2006/9/27に、測定姿勢と血圧計を変更していますので、厳密には、その前後の比較はできないのですが……。ここのところ部屋の暖房の設定温度を１度上げたところ、血圧に影響があったもようです。詳しくは平均の平均を示した、下のグラフをご覧下さい。 ...続きを見る |
2008/01/18 01:55 |
１時間血圧・脈拍【生データ】：20080112
先日公表した１時間にわたる血圧・脈拍の測定値です。フリーのデータとして公表します。学生さんのレポートなどの参考になれば幸いです。なお、測定間隔は『ほぼ』１分です。最大、10秒程度の誤差があるものと思ってください。 ...続きを見る |
2008/01/12 05:28 |
血圧と脈拍の関係、連続測定時【グラフと検定】
『等倍』でご覧下さい。 １時間にわたって連続測定した血圧と脈拍について、その関係を探ったグラフです。脈拍と関係が深いのは、最高血圧から最低血圧を引いた数字であることがわかりました。見た目では、最低血圧と脈拍の関係も、同じくらい強そうですが、例外の点が３つほどあり相関係数の上昇を妨げています。これらは、いずれも脈拍が高いのですが、その割に最低血圧が低くなっています。これらは、連続測定の終わりごろに集中して現われており、何らかの測定エラーであることが否定できません。これら３つを取り除くと最低血圧と... ...続きを見る |
2008/01/08 16:37 |
血圧・脈拍の連続測定【グラフ】
『等倍』でご覧下さい。 今朝早く、１時間かけて１分おきの、最高・最低血圧、脈拍を、連続的に測定したグラフです。最高血圧が徐々に低下し、脈拍が、後半、徐々に上昇しています。『TR』とついているのが、大まかな変化を、カルマンフィルタ、と呼ばれる手法で取り出したものです。最低血圧には、変化が認められませんでした。連続測定で、この程度ばらつくのですから、毎日１回の測定では若干不安になる気持ち、分かっていただけますか(笑)？ ...続きを見る |
2008/01/06 08:59 |
体重グラフ【ダイエット、管理限界線入り】：20080103
『等倍』でご覧下さい。 管理限界の目標を示した、ダイエットグラフ(体重)です。 2007/12/30を除いた体重の数字で、２．５σ(シグマ)の管理限界線を引いています。大晦日から新年にかけて体重の減少が足踏みしていることが、一目瞭然です。 特に、今日(2008/1/3)の朝の体重は、上方管理限界を超えていますから、今日明日は、しっかりと食事制限しないと、えらいことになりそうです。昼飯までは、しっかり食べてしまいましたから、夕食が、寂しいものになるのは避けられないでしょう。それで、明日の朝、... ...続きを見る |
2008/01/03 15:17 |
前の記事の訂正……
最高・最低血圧の差と、脈拍との相関係数は、 仰向けで、+0.075位ですから、差と脈拍は正の相関関係があることになります。 統計的に有意かどうか、ギリギリのところです。 座位では、-0.12〜-0.14位で、負の相関関係があることになります。 こちらは、かなり積極的に主張できる有意さです。 ...続きを見る |
2007/12/27 12:49 |
最高・最低血圧の差と脈拍の関係
起床時に、７回立て続けに血圧・脈拍を測定していると、なんとなく実感されるのですが、最高・最低血圧の差が大きいときは脈拍が少ないような気がしていました。 ...続きを見る |
2007/12/27 04:50 |
血圧、１年７ヶ月間【グラフ】
『等倍』でご覧下さい。 計り始めてから１年７ヶ月間の血圧の動きです。 『平均』の方は起床時７回連続測定の平均値ですから、１回しか測らなければこの２．６倍くらい、ばらつくはずです。 2006/9/27頃に、測定姿勢変更〜仰向け→座って、と、血圧計を新しくしていますから、それから後だけ見ると、秋から冬にかけて、血圧が漸増しているように見えますが、年周期があるかどうかは来年の春、暖かくなったころ、血圧が低下しないと、なんとも言えません(笑)。 ...続きを見る |
2007/12/22 23:39 |
アクセスの多かった記事(2007/12/21朝現在)
５万アクセスを記念して、今までの記事の中で、アクセスの多かったもの(300アクセス以上)を、列挙しておきます。 ------------------------------------------------ 写真の集合体を作れるフリーソフト 2006/6/28 329アクセス ...続きを見る |
2007/12/21 08:23 |
血液検査の結果【グラフ】
『等倍』でご覧下さい。 ここ７年間の、 HbA1c(血糖の長期的な指標である検査項目)、 γ-GTP(肝臓の傷められ具合の指標である検査項目)、 中性脂肪(カロリーの摂り過ぎと日頃の運動の熱心さの指標)、 のグラフです。 ...続きを見る |
2007/12/21 01:13 |
１年１０ヶ月間の体重変化【グラフ】
『等倍』でご覧下さい。 約１年１０か月間にわたる、体重の変化です。 2005年の10月頃は81キロあったので、ＢＭＩで25付近、体重で72.3キロ付近に抑えていられれば、上出来なのでしょう。負け惜しみかな(笑)。 ...続きを見る |
2007/12/12 14:56 |
体脂肪重量が減少に転じた【グラフ】
『等倍』でご覧下さい。 ...続きを見る |
2007/12/06 11:20 |
７回連続測定時の血圧・脈拍の変動
『等倍』でご覧下さい。 朝起きた時に、７回連続して測定した、最高・最低血圧・脈拍が、その日の平均値からどの程度隔たりがあったかを見ている「箱ヒゲ図」です。 V1〜V7が最高血圧、V8〜V14が最低血圧、V15〜V21が脈拍です。 最高血圧は１回目が、最低血圧は１〜２回目が、脈拍は１〜２回目と７回目が、平均からずれているようです。有意差と言えるかどうかは、最高血圧の１回目を除いては疑問があるのですが、大勢のデータを処理する場合、５回計って、その３，４回目の平均をとるくらいの用心がいるかもしれ... ...続きを見る |
2007/12/04 06:30 |
統計(直前の記事のデータを見る“Ｒ”code)
この寸前の記事をファイルにしたら、“Ｒ”を起動して、以下の「Ｒ code」を「Ｒ console」に、コピー and ペーストしてみてください。次の記事の箱ヒゲ図の横並びが出るはずです。 ---------------------------------------- x<-read.table("c:\\BP\\bpp",sep=",") y<-t(x) z<-matrix(rep(0,21*268),nrow=21,ncol=268) for(i in 1:268){ a<-y[... ...続きを見る |
2007/12/04 06:29 |
統計(血圧データ【異常値除去済み】)
起床時に、立て続けに７回測定した、最高血圧、最低血圧、脈拍について、同じ日の中での変動がどうであるかを、“Ｒ”を使って見るためのデータです。 (c:\)ドライブの(BP)ディレクトリ「(c:\BP\)」に、「bpp」というファイル名で、点線より下の、数字とコンマを登録して下さい。 ---------------------------------------------- 109,106,108,109,107,104,98,73,71,77,74,78,72,75,73,82,85,84... ...続きを見る |
2007/12/04 06:27 |
体重の変化：増加し続け！：グラフ
『等倍』でご覧下さい。 2007/7/1からの体重の変化です。 ...続きを見る |
2007/11/30 11:07 |
“Ｒ”の参考書
“Ｒ”の参考書は、いろいろあるのだけれど、私が統計数理研究所の公開講座で“Ｒ”を習った時の先生である、金 明哲さんの「Ｒによるデータサイエンス」(森北出版)は、統計学の基礎＋αくらいの知識のあるパソコンのユーザー（PC経験中程度）には向いている本です。ただ若干誤植が目立つのですが、ホームページの方でフォローしているそうです。ただなにせ高い(3600+税)けれど、“Ｒ”の本は軒並み高いからなぁ(嘆息)。 ...続きを見る |
2007/11/15 08:08 |
統計_07_005（Ｒ code）
この前のデータを「ｃ」ドライブの下の「BP」ディレクトリに「bp_day_av」という名前のファイルにした上で、点線以下を、点線の区切りごとに「Ｒ code」ウィンドウにペーストします。【点線は入力しない】 -------------------------------- #データの読み込み x<-read.table("c:\\BP\\bp_day_av",header=T) #散布図(相関図)の出力 plot(x) --------------------------------... ...続きを見る |
2007/11/14 13:19 |
血圧・脈拍の相関図にもちいたデータ（フリーのデータとします）
先の血圧・脈拍の相関図 に用いたデータ（平均値 にしたもの）です。“Ｒ” で、いろいろな分析をして みたい方のための「フリー」 のデータとして公開します。 最初が(2007/4/1)です。 起床時に７回連続測定 した値によるものです。 ...続きを見る |
2007/11/14 08:23 |
最高・最低血圧、脈拍の分布図：グラフ
『等倍』でご覧下さい。 2007/4/1〜2007/11/14 のその日、立て続けに７回計ったデータの平均値による、最高血圧、最低血圧、脈拍、の関連性の図です。 欠測の日があるのでｎ＝223、です。 ...続きを見る |
2007/11/14 05:30 |
アクセス数の曜日依存性の有無
毎日のアクセス・レポートで、アクセス数の平均値に大きな曜日差があるかどうか検討したのですが『大きな』差は見られませんでした。ただ、日曜日に少なく月曜日に多いことが分かりましたが、アクセス数としては２割くらいの差しかありませんでした。 ...続きを見る |
2007/11/12 09:21 |
統計_07_003(血圧等の連続測定における変動)
先日の血圧等のファイルを読み込んで、 各日の何回目に計ったかによって 分けた『箱ひげ図』を出力するものです。 最高血圧(3つに分かれた左側)の 1回目が高めになっていることが分ります。 これは昨年と、ほぼ同様の結果です。 血圧を測っている方で、時間に余裕のある方は、 自分でも、似たような結果が出るかどうか、 2週間くらいにわたって、毎日、立て続けに 7回測定するのも面白いかと思います。 ------------------------------------ x<-read... ...続きを見る |
2007/11/07 08:57 |
このブログの統計“Ｒ”用codeはフリー扱いでいいです
タイトルの通りです。ただ常識的に扱って下さいね。 ...続きを見る |
2007/11/05 09:53 |
統計_07_002(血圧等、グラフ化)
点線より下を『Ｒ code 窓』に コピー and ペーストして実行すると、 (“c:\BP\”)に作った(“bp_k7.txt”)が あれば、最高血圧・最低血圧・脈拍の データを『１日連続７回測定の平均』をとって 『７日移動平均』して『標準化』したグラフが 表示されます。カラーです。 『標準化』してやることによって、様々な項目の 平均値やバラツキの大きさが異なるデータを、 類似の基準で眺める事ができるようになります。 ただ常識的には３項目くらいが限界でしょう。 -----... ...続きを見る |
2007/11/05 09:27 |
統計_07_001(血圧・脈拍、生データ)
点線より下を(“c:\BP\”)というディレクトリに(“bp_k7.txt”) という名前で、数値文字ファイルとする。 ...続きを見る |
2007/11/05 09:26 |
各種測定項目の分類(体重、etc)
2007/8/1〜8/25の ...続きを見る |
2007/08/25 11:12 |
アクセス時刻の分布
『等倍』でご覧ください。 このブログへの、アクセス時刻の分布です。2006/6/17〜2007/3/16についてのアクセスレポートから作成したものです。 このブログへのアクセスは、健康と朝顔の両日記のUPのタイミングを除けば、特別な点はないでしょうから、午前５時〜午前８時の凹みは、どこでも起こっていることかもしれませんね。 ...続きを見る |
2007/03/17 11:59 |
wt2用、Ｒ_code(バグ取り後)
x1)&(sitac_old0){sitac1){uefg1)&(uec_old>0)){lines(c((i-1),i),c(yuec[(i-1)],yuec[i]),col=3,lwd=2)} if(uec1)&(sitac_old2.7){ sitacovfg[i]0){sitac1){uefg1)&(uec_old>0)){lines(c((i-1),i),c(yuec[(i-1)],yuec[i]),col=3,lwd=2) if (abs(yuec[i])>2.7){ uec... ...続きを見る |
2007/02/21 05:31 |
BPmax用、Ｒ_code(バグ取り後)
x1)&(sitac_old0){sitac1){uefg1)&(uec_old>0)){lines(c((i-1),i),c(yuec[(i-1)],yuec[i]),col=3,lwd=2)} if(uec1)&(sitac_old2.7){ sitacovfg[i]0){sitac1){uefg1)&(uec_old>0)){lines(c((i-1),i),c(yuec[(i-1)],yuec[i]),col=3,lwd=2) if (abs(yuec[i])>2.7){ uec... ...続きを見る |
2007/02/21 05:30 |
バグのお詫びと体脂肪重量の管理図(decision limit cusum 法)
以前公開した『Ｒ code』に、えらいバグがあったので、今日(2007/2/21)中に改訂版をUPする予定です。 『date』が予約語だった上に、未定義のまま使っていました。 ワークエリアが全く空白の状態で作動させるなどすると、エラーとなったはずです。 まことに申し訳ありませんでした。 ...続きを見る |
2007/02/21 05:10 |
下４つの記事のサンプルグラフ：解説あり
『等倍』でご覧下さい。 下４つの記事のデータと『Ｒ code』で作った『decision limit cusum 管理図』です。 『Ｒ code』の１行目にある、読み込むファイル名を変えることで、かなり広範に適用できると思います。 使い心地が悪い点を指摘していただければ、改善版を作りますので、コメントしてください。 カラーによる表現を使っていますのでその点、御了承ください。 黒の細い実線が、管理すべきデータ。 緑のやや太い実線が上方CuSum。 青のやや太い実線が下方CuSum。 ... ...続きを見る |
2007/02/14 15:09 |
BPmax(最高血圧データ：R code用)
点線から下を（"c:\kenko\"）というディレクトリに （"BPmax"）という名前で、数字のファイルとして 登録してください。 ...続きを見る |
2007/02/14 15:06 |
wt2(体重データ：R code用)
点線から下を（"c:\kenko\"）というディレクトリに （"wt2"）という名前で、数字のファイルとして 登録してください。 ...続きを見る |
2007/02/14 15:06 |
朝顔花数。グラフを滑らかにしようとする試み。
2006年の朝顔の花数のグラフ(細い線)と、それを滑らかにしようとした線(太い線)です。 6/24(最初)から9/1頃までは周期的に増減しているようですが、統計的に、その傾向を取り出すことは出来ませんでした。主に、初めの頃、花数０(ゼロ)が何回も出ているからです。 ...続きを見る |
2007/01/29 03:44 |
統計_028(_027、グラフ見本)
下(統計_027)の記事で作ったグラフです。 動作確認のために『等倍』で、ご覧下さい。 グラフの下の説明は『ペイント』で付けたもので、 『R_code_009』には含まれていません。 ...続きを見る |
2006/12/28 02:54 |
統計_027(R_code_009)：平滑化と標準化：白黒
点線より下を『Ｒ code 窓』に コピー and ペーストして実行すると、 以前に(“c:\BP\”)に作った(“k_dx7.txt”)が 残してあれば、最高血圧・最低血圧・脈拍の データを『１日連続７回測定の平均』をとって 『７日移動平均』して『標準化』したグラフが 表示されます。白黒です。 『標準化』してやることによって、様々な項目の 平均値やバラツキの大きさが異なるデータを、 類似の基準で眺める事ができるようになります。 ただ常識的には３項目くらいが限界でしょう。 ... ...続きを見る |
2006/12/28 02:52 |
統計_026(_024、025、のグラフ見本)
『等倍』でご覧下さい。 “統計_024(最高・最低血圧：生データのうち公式記録分)” “統計_025(R_code_008)：白黒、カラー編はなし” で作った『濃度･等高線表示の図』です。 データと『R_code_008』が正常作動したかどうかの確認のためにご覧下さい。 ...続きを見る |
2006/12/25 05:33 |
統計_025(R_code_008)：白黒、カラー編はなし
x<-read.table("c:\\BP\\max_and_min_nama.txt",sep=",") a<-matrix(c(rep(0,200*200)),nrow=200,ncol=200) for (i in 1:1526){q<-x[i,1] r<-x[i,2] a[q,r]<-a[q,r]+1 } b<-matrix(c(rep(0,60*45)),nrow=60,ncol=45) for (i in 81:140){for (j in 56:100){q<-i-8... ...続きを見る |
2006/12/25 05:32 |
統計_024(最高・最低血圧：生データのうち公式記録分)
私の最高血圧・最低血圧の公式記録された生データ。 (“c:\BP\”)なるディレクトリに (“max_and_min_nama.txt”)という名前で 点線より下を、数値文字ファイルとする。 そのファイルは、次の(上の)記事 “統計_025(R_code_008)：白黒、カラー編はなし” の「Ｒ code」で濃度表示＋等高線表示の図と する事ができる。 ...続きを見る |
2006/12/25 05:31 |
統計_023(グラフ見本：血圧･脈拍の長期変動)
前の記事によって作成したグラフです。 白黒とカラー版の２つを載せます。 『等倍』でご覧下さい。 「Ｒ」がきちんと走ったかのチェックに ご利用下さい。 ...続きを見る |
2006/12/15 02:45 |
統計_022(R_code_007)：カラー
x<-read.table("c:\\BP\\max_min_pulse.txt",sep=",") d<-as.Date(x[,1]) a<-matrix(c(rep(NA,206*2)),nrow=206,ncol=2) b<-matrix(c(rep(NA,206*2)),nrow=206,ncol=2) c<-matrix(c(rep(NA,206*2)),nrow=206,ncol=2) a[,1]<-x[,2] b[,1]<-x[,3] c[,1]<-x[,4] me... ...続きを見る |
2006/12/15 02:43 |
統計_021(R_code_006)：白黒
x<-read.table("c:\\BP\\max_min_pulse.txt",sep=",") d<-as.Date(x[,1]) a<-matrix(c(rep(NA,206*2)),nrow=206,ncol=2) b<-matrix(c(rep(NA,206*2)),nrow=206,ncol=2) c<-matrix(c(rep(NA,206*2)),nrow=206,ncol=2) a[,1]<-x[,2] b[,1]<-x[,3] c[,1]<-x[,4] me... ...続きを見る |
2006/12/15 02:42 |
統計_020(血圧等_04)
日付つきの最高血圧・最低血圧・脈拍の各日平均値。 (“c:\BP\”)なるディレクトリに (“max_min_pulse.txt”)という名前で 点線より下を、数値文字ファイルとする。 これは、次の(上の)記事“統計_021(R_code_006)：白黒” あるいは“統計_022(R_code_007)：カラー” の「Ｒ code」でグラフ化する事ができる。 ...続きを見る |
2006/12/15 02:34 |
統計_019(統計_009(R_code_001))訂正済み
median.na<-function(x){ y<-na.omit(x) n<-length(y) if (n==0) a<-NA else a<-median(y[0:n]) a } a<-NA b<-13 c<-c(1,2,3,4,5) d<-c(1,2,3,NA,4,5) e<-c(NA,NA,NA,NA,NA) median.na(a) median.na(b) median.na(c) median.na(d) median.na(e) mean.na... ...続きを見る |
2006/12/14 02:09 |
統計_018(統計_008(matplot))訂正済み
この前の「統計_007(体重)」という記事のデータの登録が終わったら、「Ｒ console」プロンプト『左端の(＞)と(＋)』に対して、点線より下のように入力していくと、プロンプトのない行が出力され、最終的にグラフィックウィンドウがひらき３色の実線のグラフが描かれるはずです。赤が７日移動中央値、緑が７日移動平均です。黒は元の体重のデータ。注意深く見るとグラフの黒い実線はところどころ飛んでいることが分かります。 単に本当に作動するかどうか見たい方は、この後の「統計_009(R_code_001)... ...続きを見る |
2006/12/14 02:09 |
間違いの訂正のお知らせ
このブログの『統計_008』と『統計_009』に、重大な誤りがありました。 NAを含むデータの平均を計算すべき関数の定義・その他が中央値をとるものと入れ替わっていました。お詫びして、この上の『統計_018』と『統計_019』に投稿しなおしておきます。 ...続きを見る |
2006/12/14 02:08 |
統計_017(R_code_005)
#白黒・点線 x<-read.table("c:\\BP\\waip.txt",sep=",",header=T) k="平滑化した体重・血圧・脈拍" matplot(x,type="l",lty=1:4,col=1,ylim=c(-3,6),main=k,lwd=1.5) legend(1,6,colnames(x),col=1,lty=1:4,lwd=1.5) ...続きを見る |
2006/12/10 13:23 |
統計_016(血圧等_03)
点線以下は、2006年12月9日の『体重・血圧・脈拍』のグラフを“Ｒ”で作るためのデータ。（“c:\BP\”）に（“waip.txt”）という名前で数字テキストファイル（ヘッダ付き）とする。 ...続きを見る |
2006/12/10 13:23 |
統計_015(R_code_004)
ww<-function(a){y<-na.omit(a) n<-length(y) if (n==0) a<-NA else a<-y[0:n] a} x<-read.table("c:\\BP\\k_dx7.txt",sep=",") y<-as.matrix(x) a<-y[1:196,1] a j<-ww(a) j qq<-NA qq<-matrix(c(rep(0,7*31)),nrow=31,ncol=7) qq for(h in 1:7){ a<-y[1... ...続きを見る |
2006/12/04 02:04 |
統計_014(R_code_003)
x<-read.table("c:\\BP\\k_dx7.txt",sep=",") boxplot(x[,1:7]) ...続きを見る |
2006/12/04 02:03 |
統計_013(血圧等_02)
点線より下を(“c:\BP\”)というディレクトリに(“k_dx7.txt”) という名前で、数値文字ファイルとする。 ...続きを見る |
2006/12/04 02:03 |
体重のグラフ・トレンドの抽出
８月中旬から１１月中旬までの、毎日の体重と、それを本格的な方法で平滑化した『トレンド』の変化です。『等倍』で見てください。 ...続きを見る |
2006/11/19 15:12 |
統計_011(日次平均)
まず前の記事「統計_010(血圧等)」のデータ部分を「BP」にテキストファイル「ketuatu.txt」として登録した上で・・・。 以下の点線より下の行で、先頭に「＞」と「＋」にある行を「Ｒ console」に入力してください。 行先頭の「＞」と「＋」は入力の必要はありません。 行先頭に「#」がある行は「注意書き」ですので入力の必要はありません。 前回と同じく動作だけ見たい方は次の(上の)記事「統計_012(R_code_002)」を「Ｒ console」にコピー・アンド・ペーストしてく... ...続きを見る |
2006/11/13 03:01 |
統計_012(R_code_002)
mean.na<-function(x){ y<-na.omit(x) n<-length(y) if (n==0) a<-NA else a<-mean(y[0:n]) a } u<-function(i,j){a<-abc[i:(i+6),j] a } abc<-read.table("c:\\BP\\ketuatu.txt") gg<-matrix(c(rep(0,175*3)),nrow=175,ncol=3) for(i in 1:175){k<-(i-1)*7+... ...続きを見る |
2006/11/13 03:01 |
統計_010(血圧等)
点線以下は、起床時に７回連続して測定した最高血圧・最低血圧・脈拍の生データです。 今回はこれの毎日の平均のグラフを描きます。 まず、このデータ(点線より下)を(c:\BP\ketuatu.txt)なるファイルにしてください。 -------------------- 122 90 78 111 83 76 109 87 76 108 85 74 103 82 77 107 82 75 113 85 79 118 84 76 112 80 71 120 89 70 11... ...続きを見る |
2006/11/13 03:00 |
統計_009(R_code_001)
median.na<-function(x){ y<-na.omit(x) n<-length(y) if (n==0) a<-NA else a<-median(y[0:n]) a } a<-NA b<-13 c<-c(1,2,3,4,5) d<-c(1,2,3,NA,4,5) e<-c(NA,NA,NA,NA,NA) median.na(a) median.na(b) median.na(c) median.na(d) median.na(e) media... ...続きを見る |
2006/11/06 03:19 |
統計_008(matplot)
この前の「統計_007(体重)」という記事のデータの登録が終わったら、「Ｒ console」プロンプト『左端の(＞)と(＋)』に対して、点線より下のように入力していくと、プロンプトのない行が出力され、最終的にグラフィックウィンドウがひらき３色の実線のグラフが描かれるはずです。赤が７日移動中央値、緑が７日移動平均です。黒は元の体重のデータ。注意深く見るとグラフの黒い実線はところどころ飛んでいることが分かります。 単に本当に作動するかどうか見たい方は、この後の「統計_009(R_code_001)... ...続きを見る |
2006/11/06 03:18 |
統計_007(体重)
点線より下の数値とNAからなるデータを(c:\BP\)に(taizyu.txt)という名前のファイルとします。 次の記事で、欠測値のある時系列データの「移動中央値(メディアン)」と「移動平均」をグラフ化する方法を示します。 ----------------------------------------------- 75.4 75.3 75.3 74.7 74.8 74.7 74.8 74.8 74.8 74.7 74.7 74.9 74.6 74.1 74 ... ...続きを見る |
2006/11/06 03:17 |
統計_004(Ｒ)
「max.txt」、「min.txt」、「p.txt」の３つのファイルが出来たら、適当なサイトで『Ｒ 統計』で検索して、統計用フリーソフトウェア『Ｒ』をダウンロードしてインストールする。 さまざまなバージョンがあるが、私がここで使っているのは『バージョン 2.20.51006.0』である。 ただし、日本語向けの変更が加わっている部分変更版である。日本語機能は使えたほうが便利であるので、日本語向けの処置が施してある版をダウンロードするにこしたことはない。 ダウンロード先について、長期にわたっ... ...続きを見る |
2006/11/01 02:28 |
統計_005(ｂｏｘｐｌｏｔ)
先の３つのファイルが出来『Ｒ』のダウンロードが済めば、とにかく実際、使ってみよう。 ...続きを見る |
2006/11/01 02:28 |
統計_006(ｍｅａｎ)
まず「max.txt」を、read.table を使って max に読み込む。 ...続きを見る |
2006/11/01 02:28 |
統計_001(ｍａｘ)
記事最下部の数値データ文字列（コンマを含む）を、 マウスで始めから終わりまで印をつけ（ドラッグし）コピーし、 ...続きを見る |
2006/11/01 02:27 |
統計_002(ｍｉｎ)
82,92,73,64,79,65,81,78,70,78,78,76,72,67,78,90 88,89,75,68,80,68,77,71,71,77,73,74,71,67,74,95 89,84,73,72,78,68,73,74,71,75,82,70,74,68,75,86 83,85,74,70,78,67,78,75,69,77,81,65,75,67,77,93 88,85,78,73,76,67,81,75,72,73,74,68,75,67,81,91 9... ...続きを見る |
2006/11/01 02:27 |
統計_003(ｐ)
87,81,75,79,67,74,72,80,86,71,71,75,69,70,76 ,76 82,82,74,71,66,79,68,78,78,74,70,77,68,77,80 ,76 81,78,75,71,63,90,71,76,85,83,71,76,74,76,77 ,76 87,77,74,87,66,72,77,73,76,78,75,72,71,72,77 ,80 81,81,78,75,66,78,73,78,95,78,75,77,71,84,80... ...続きを見る |
2006/11/01 02:27 |
トップへ ｜ テーマトップへ ｜ テーマランキング一覧へ
