10秒ごと更新，2分の移動平均？ RT @Mihoko_Nojiri: ２分まつとかいてあったよね。RT @h_okumura: ソフトバンクの線量計付き携帯の動画 http://t.co/a0KM0x09 10秒もたたないうちに0.05μSv/hとぴったり表示。ちょっと怪しい。

posted at 18:29:05

あ，デモ用でしたか＞あの動画 RT @surblue: @Mihoko_Nojiri @h_okumura 発表会場のものはデモ用で、実際は通常で2分計測、最長40分まで計測できるそうです。

posted at 18:31:39

割れ瓶 RT @Mihoko_Nojiri: ホロビッツがきてひどい演奏をしていった時の吉田秀和の評論ってのは今でも覚えてる。

posted at 11:58:42

そうなんです。福島大のは酷かった RT @Mihoko_Nojiri: 福島大のって全く定量的な情報がないやつじゃない。RT @h_okumura: 放射線教育、自前の副読本つくる http://t.co/yfrrLCAl 「普通の大人がネットや本などの情報を活用し」とあるが…

posted at 10:14:21

納得です RT @Mihoko_Nojiri: 表でもコメント。別に dis ってるわけではなく、これをやったら劇的に結果が改善するという leaf　さんの期待がまちがっているということを説明しています。@h_okumura http://t.co/kFb29cKs

posted at 19:23:58

おかーさんが口とんがらせたマンガが出てきたんですね^^ RT @Mihoko_Nojiri: これ秘密なんだけどさー、娘のそーごがくしゅうの課題がが放射線になってさー、それで「放射線」「測り方」でぐぐって、最初にでてきた漫画をクリックしたらあちゃーになったんだってー。

posted at 21:58:22

文科省への問合せ感謝！ RT @ichimiyar: @hayano @Mihoko_Nojiri @h_okumura 話をする中で、Sv/h, Gy/h両方測れるTCS-171で切り替えて測ると2割違うという事もご存知でしたし、JIS Z4511で校正されたSv/h単位の…

posted at 20:12:02

大入り満員でした @Mihoko_Nojiri @hayano @AkiraOkumura

posted at 18:37:08

これは国じゃなくて福島県なんでしょうね RT @heart8255: 黒猫先生達も知らないって…公開されていなかったということでしょうか？これってNHKがやる事じゃなくて国がやるべきことでは？？RT @Mihoko_Nojiri これは知らなかった。。。…

posted at 08:53:15

.@Mihoko_Nojiri コピペできましたか，よかった。少し前に全データを落としたもの（UTF-8）→ http://t.co/aZh1zRWO

posted at 22:38:24

東海アマ！ RT @Mihoko_Nojiri: 再掲しておくけど、「山崎パンから検出」してしまったデータはここ。http://t.co/XukS8eOe　KH-EL25 売ってるのもここ。検証　→　http://t.co/HIcqBJEI

posted at 09:13:06

パン，やさしおの計算 http://t.co/4yo24r40 と同様にしてみれば数十Bq/kg程度ですか。ただγが出るのは壊変の1割 RT @g4mishit: @Mihoko_Nojiri …日本食品標準成分表2010によると100gあたりKは97mgとのこと…

posted at 10:40:12

うちもFAXは間違いしか届かない RT @Mihoko_Nojiri: うちに一番くるファックスは日○T○○○○○宛。

posted at 19:41:18

@Mihoko_Nojiri 便益分析は経済学者では？

posted at 16:24:24

どこかに物理学者でもあるという経済学者がおられた気がします。何とか話法の人 RT @Mihoko_Nojiri: じゃあ私が欲しいのは科学に明るい経済学者だな。RT @h_okumura: ノビー^^; RT @buvery: 経済学者と言えば、ノビー。除染など、全部無駄と言われ

posted at 16:34:41

@Mihoko_Nojiri 検出限界10Bq/kgとすれば誤差はその1/3で3Bq/kgほどになりそうですが（0〜10Bq/kgのレンジで誤差が一定として）

posted at 16:44:31

積算線量計ならいいですね。あるいはポケガのようにずっと平均をとり±誤差を付ける RT @Mihoko_Nojiri: 子供にもたせるならいいかな。RT @h_okumura: 40cpm/μSv/h程度… エアカウンターSの使用感 http://t.co/qgpS8R1t

posted at 14:33:01

何らかの平滑化をしているなら分散からcpmわかりませんがどうなのかな RT @Mihoko_Nojiri: まあ今ぽちって、来週分散から有効CPM 出すから、みなさんはそれみてから買えばどうでしょうかー＞エアカウンターS

posted at 14:37:30

相関も調べればいいですね RT @Mihoko_Nojiri: おばさんはそういうの得意だからRT @h_okumura: 何らかの平滑化をしているなら分散からcpmわかりませんがどうなのかな RT @Mihoko_Nojiri: まあ今ぽちって、来週分散から有効CPM 出すから

posted at 14:41:24

約40cpm/μSv/hはここから取らせていただきました RT @warafuni: @Mihoko_Nojiri @h_okumura 横から失礼します。エアカウンターSの測定実験（http://t.co/sl0YcPdf）はご覧になられたでしょうか？（下部：感度推定）…

posted at 14:52:19

エアカウンターSの @richeehauser さんのデータ，以前ポアソンとの比較に使わせていただきました http://t.co/X9XHVOpZ これからも約40cpm/μSv/hが確認できます @Mihoko_Nojiri

posted at 15:03:09

いえ分散共分散からcpm推定たいへん興味あります RT @Mihoko_Nojiri: ぽちるまでもなかったかRT @h_okumura: エアカウンターSの @richeehauser さんのデータ，以前ポアソンとの比較に… http://t.co/X9XHVOpZ

posted at 15:06:01

いいなぁそれくらいの線源ほしい RT @Mihoko_Nojiri: 40cpm /(μSv/h) では0.1 の環境で２分はかっても8カウントいかないので大変だが, 手持ちのお道具で5μSv/h はやれるから、安定して数字がw @h_okumura

posted at 15:12:31

そうか。なら2分ごとにリセットすれば分散だけでcpmがわかる RT @richeehauser: @Mihoko_Nojiri 赤ランプが消えるまでの約２分間は約40cpm/(μSv/h)で値が更新されます。PCでシミュレータを作って確認しました。総カウント数÷所要時間…

posted at 16:04:31

@ystricera さんも三重県でしたか RT @Mihoko_Nojiri: 三重県のほうには足をむけてねれんのう。。たかよしさんにビール券を送りたい。@ystricera

posted at 20:02:40

@_@; 歳がわかる RT @Mihoko_Nojiri: だってキノコなんだもの、セシウムでちゃう♪アタック〜アタック〜♪…

posted at 15:18:42

先月はお元気そうなお姿を見ました RT @Mihoko_Nojiri: そういえば一松信先生はお元気であろうか、基礎物理学研究所にいた時は、一般講演とかに時々おいでになっていて、ご案内しようとする、いえいえ大丈夫みたいな、最後におみかけしたのは、７年前くらいかな。

posted at 18:15:00

物理以前の統計学の話で，要は母集団のSDがσ RT @morino_kumasan: ふむ。物理学者とつきあう前は、sigmaなんて言葉使わんかった。うちらの分野じゃ、まだSDの方が通じるかな。 RT @Mihoko_Nojiri: … http://t.co/fkUn2wDc

posted at 08:42:35

第5話！マンガでも男前の勝川先生 RT @Mihoko_Nojiri: あ、コミックりゅうにみそさんの新しいマンガでました。今回は勝川先生　@katukawa が主役でNa I の食品検査の話　 僕と日本が震えた日 - http://t.co/MtBzA1EQ

posted at 21:13:50

昨日雪崩で3人死亡したところ，最大7μSv/h！ RT @Mihoko_Nojiri: 玉川温泉の岩盤浴って地熱じゃなくて放射線なんだけど。あ、ガイガーカウンターの音が聞こえる。http://t.co/fejWqxKv

posted at 08:39:00

ありがとうございます。分解能ないならCs134+137を出すのは良いと思います。測定器の内蔵時計で事故からの日数で比を求めて別々に出すことを測定器レベルで暗黙にするならそう断らないと，福島起源でないものを測られたら困るのでは？ @mw_mw_mw @Mihoko_Nojiri

posted at 11:28:58

Cs137が負になったらCs134:Cs137=1:1を仮定して正にするなんて余計なおせっかいと思うのですが RT @kaztsuda: @Mihoko_Nojiri …和田さんがメーカーのエンジニアに聞いた情報では https://t.co/2Pz7cCWK

posted at 20:54:44

時計でも入っているんでしょうか RT @Mihoko_Nojiri: Cs 134 がでたら、Cs 137 を 0.8: 1 を仮定してだすってのはありという気もしないでもない。。。RT @h_okumura: Cs137が負になったらCs134:Cs137=1:1を仮定して正…

posted at 20:57:20

子どものころガンマといえば→ http://t.co/ezkbBs1i RT @Mihoko_Nojiri: ガンマ来てるな。RT @ystricera: 2号機格納容器内部工業用内視鏡映像一つ目 http://t.co/lF0UMdir

posted at 13:22:16

minimum detectable activityが書いてありますね RT @Mihoko_Nojiri: AT1316 …誤差も検出限界もかいてないけど、なんじゃい。RT @MMsoaring: … http://t.co/jlWLHsZi

posted at 22:16:37

http://t.co/sBUuIwhp の最初の図で黒い線分の上の面積を求めるだけです RT @Mihoko_Nojiri: 緩募　コベル法が何をやっているのか one tweet で説明

posted at 22:02:20

仕事にお使いになるならぜひMac版のRを:-) RT @Mihoko_Nojiri: @h_okumura これは仕事にも使えそうです。

posted at 08:06:03

@Mihoko_Nojiri あれ変ですね。wbc.Rをやっておられるのですか？

posted at 11:48:48

@Mihoko_Nojiri wbc.Rの20行目に rng = 200:850 とかあるでしょうか？ どこでrngが見つからないんでしょうね。不思議。MacですかWindowsですか？

posted at 11:51:55

@Mihoko_Nojiri きっとMac→Winで行末コードが化けて行がつながってしまったのかもしれませんね

posted at 12:10:56

@Mihoko_Nojiri そうでしたか。変数名を適当につけているので脳内置換されてしまったのだと思います。三つ山を順にm1,m2,m3と呼ぶかC-134がm1,m2でC-137がm3にするかの違いだと思います

posted at 12:48:43

@Mihoko_Nojiri 手抜きでご迷惑をおかけしています。http://t.co/Or6RuW23 や http://t.co/Spt6MvsM の例もご参照ください

posted at 13:02:21

@Mihoko_Nojiri sp = 1200 * sp とか？

posted at 13:08:52

@Mihoko_Nojiri text(x座標,y座標,メッセージ)です。メッセージはcat(m1,"keV")のように単位を結合するとわかりやすいかも

posted at 13:35:13

@Mihoko_Nojiri すいません，catじゃなくpasteでした^^; text(50,3,paste(format(m1,digits=2), "keV"))みたいに

posted at 16:41:41

@Mihoko_Nojiri あ，fit=gnls(...)した値でしたらa=coef(fit)してa['m1']です。wbc.Rの最後のほう参照

posted at 16:46:22

今朝撮ったものです。最近毎日どんどん若くなってます^^; RT @Mihoko_Nojiri: @h_okumura おおお、了解です。奥村さんアイコン若返りましたね。

posted at 16:48:08

@Mihoko_Nojiri ベクトルの長さの違うものが混じっているのでしょう。レンジrngを変えてほかと辻褄があっていないとか

posted at 17:46:05

@Mihoko_Nojiri 例えばfitが更新されていないとか？

posted at 17:48:48

@Mihoko_Nojiri いやリセットは必要ありません。どの行でエラーになるかわかればいいのですが。例えばテキストエディタからRコンソールにコピペすればエラー行ですぐ文句を言われるのでどの行が問題かわかりますよね

posted at 17:50:48

@Mihoko_Nojiri メールで送っていただければデバッグします

posted at 17:52:59

@Mihoko_Nojiri 例えばfit=gnls(...)が収束しないでfitに周回遅れのものが残っているということはありそうですね

posted at 18:06:06

@Mihoko_Nojiri 微妙なところで収束していたのかも。nlsTol=1e-3くらいにすれば収束しやすくなります。収束したりしなかったりするのはモデルの選択を間違えているのかもしれませんが

posted at 18:08:43

@Mihoko_Nojiri bin数には依存しないはずですが不思議ですね

posted at 18:10:23

@Mihoko_Nojiri 了解です

posted at 18:11:23

@Mihoko_Nojiri いまお送りしました。データの読み方のところから添削しちゃいました。

posted at 18:31:46

@Mihoko_Nojiri あ，ファイル名のところを一時的に書き換えたのがそのままでした

posted at 18:36:17

@Mihoko_Nojiri CSVの頭と尻尾にゴミがあるのでskipとnrowsで対応，それから1行に3通りのデータがあるのでそのどれを使うかを指定しました

posted at 18:37:57

@Mihoko_Nojiri 走ってますが，レンジをどう設定するとこけますか？

posted at 18:42:48

@Mihoko_Nojiri gnls()の中も若干書き換えました。いまから食事に…

posted at 18:45:24

@Mihoko_Nojiri れれれ，どうしてコピペとファイル読み込みで違うんでしょう？やっぱり行末コードが化けている？

posted at 19:33:49

@Mihoko_Nojiri はてはて，コンピュータは気難し屋です

posted at 19:52:29

@Mihoko_Nojiri 誤差の出し方はさきほどメールしました

posted at 21:12:00