違うところに書き込んでたようなので移動しました
453 :132人目の素数さん:2015/03/25(水) 11:42:34.89 ID:XFYYJvdYK
R初心者です
nls.lmの使い方を教えてください
いままでフィッティングはエクセルでやってましたが
誤差評価をつけろと言われてRを始めたところ全然分からなくて困っています
時間に対して観測値が変化するのを普通の式でフィッティングするのはできるのですが
一定の時間の後に式が変化するときの場合分けはどうしたらよいのでしょうか?
具体的に言うと、1回投与に対するフィッティングはできるけど
数回に分けて投与する場合のフィッティングが分かりません
推定値をgetPred(parS, t)という関数で求めるとして、
getPred関数中のif文でtの値による場合分けをして
サブルーチンを呼び出そうと思ったのですが
nls.lm中にはtの値はベクトルで渡されていて
nls.lmが最少化する目的関数の計算中に変化させているt[?]の状態で
サブルーチンに渡すことができず
単純にtをgetPredに渡してif(t<1)とかやると
「条件の長さが2以上なので・・・」
というエラーメッセージがでてしまいます
うまく場合分けする方法はあるでしょうか
454 :453:2015/03/25(水) 11:49:25.36 ID:XFYYJvdYK
補足です
if文ではなく、
getPredの中に場合分けに対応した複数の式を書いて
それぞれの頭に
ceiling(min(1,max(0,t-1)))*
のような0か1になる場合分けの式を書き込んでもnls.lmは分けてくれませんでした
【R言語】統計解析フリーソフトR 第5章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
451132人目の素数さん
2015/03/25(水) 14:04:02.55ID:CH/w1DxE■ このスレッドは過去ログ倉庫に格納されています
ニュース
- ランドセルにくぎ刺される「国に帰れ」など言われ、転校を余儀なくされた海外からの転校生 仙台市教育委員会が「いじめ重大事態」認定★2 [煮卵★]
- 「いいの?前科ついちゃうよ」万引きした女子大学生を脅す 元コンビニ店長の男(54)逮捕 ★2 [煮卵★]
- 【芸能】立川志らく「ヤクザの事務所が隣近所にあれば、強盗だとか闇バイトが来ない」 暴力団の“利点”指摘に違和感 [冬月記者★]
- 《活動休止から365日》国分太一 すべての仕事を失い、TOKIOもスピード解散…いまだ復帰叶わぬも残された“最大の謎” [Ailuropoda melanoleuca★]
- 【テレ東音楽祭】長渕剛、テレ東初出演!『とんぼ』『乾杯』熱唱へ 酒井法子は『碧いうさぎ』歌唱 [Ailuropoda melanoleuca★]
- 【サッカーW杯】24年ぶりV狙うブラジル、ハイチに3-0快勝で今大会初白星!波乗りクーニャが2発、好調ヴィニシウスが1G1A [鉄チーズ烏★]
- 【地上波/DAZNほか】 FIFAワールドカップ2026 総合スレ★106【メキシコ/カナダ/アメリカ】
- 【MLB】ドジャース vs オリオールズ ★2
- 【MLB】ドジャース vs オリオールズ
- はません
- 巨専】
- はません ★3
- 木原「高市首相がG7でボッチだったり無視された事実はない!デマを流すなー😡」 [931948549]
- /ほく(プロ)「炙りあぶ、下さい。」店主「炙りアブ?」
- 【悲報】「フェルメール展」のチケットを取れなかったジャップ民がSNSで大暴れwwwwwwwwwwwwwwwwwwwwwwww [551743856]
- 新幹線の個室でやりたいこと
- 🪖自衛官、退職しまくり😱入隊から5年以内で退職するケースが約5割 [718678614]
- 🏡ガイジ製造工場