【極秘】Rでlmer()、random slopeなど igu***さん、お世話になります。 mo

Rでlmer()、random slopeなど . . igu***さん、お世話になります。 mo…

Rでlmer()、random slopeなど . . igu***さん、お世話になります。 mo…

Rでlmer()、random slopeなど
igu***さん、お世話になります。
modelは前回お世話になったのと同じで、その続きです。よろしくお願いします。 「>ある生薬を服用すると、身体の表面で温度の分布パターンが出来るようなので、あらかじめそういうパターンが示される位置を6箇所選んだ。
>次にこの生薬を服用する前(1)、服用後30分後(2)、1時間後(3)にこれら6箇所の皮下血流量をドップラー血流計で測定して比較データを得える試行を10回繰り返した。」
これをランダムな数値で置き換えて、Rで試して練習にしています。

> setseed(3)
> area =rep(letters[1:6] 10)
> before=floor(runif(60 400 1200)) #;from 400 – 1200
> after 5=floor(runif(60 100 600)) #;from 100 – 600
> after 10=floor(runif(60 50 200)) #;from 50 – 200
> table all=dataframe(area before after 5 after 10)
> table all
___area__before_____after_5___after_10
1______a____534_____508______157
2______b___1046_____128______173
3______c____707_____501______173
4______d____662_____152_______79
5______e____881_____483______140
6______f____883_____252______114
7______a____499_____484______145
8______b____635_____370______113
(以下60までつづく)
[area_]はfixed effectですが、
[before]値は[area]によって上下するから、これはrandom slopeになるなぁと考えました。
それから、従属変数を(before+after_5+after_10)にしたmodel(A)と、
(A)の従属変数から(+after_10)を外して、従属変数を(before+after_5)にしたmodel(B)とを作って、この二つをANOVAで比較したら(+after_10)の有意差が出るかな、と考えました。
質問1:この考えは正しいですか?
それでまずこれをlmerで演算してみたら、warningをもらいました。

#####
> temp_lmer=lmer((before+after_5+after_10)~area_+(1+before|area_) data = table_all)
Warning messages:
1: In checkConv(attr(opt “derivs”) opt$par ctrl = control$checkConv :
unable to evaluate scaled gradient
2: In checkConv(attr(opt “derivs”) opt$par ctrl = control$checkConv :
Model failed to converge: degenerate Hessian with 1 negative eigenvalues
3: In as lmerModLT(model devfun) :
Model may not have converged with 1 eigenvalue close to zero: -11e-12
#####
質問2:
1項目の「unable to evaluate scaled gradient」というのはrandom slopeをつくれない、ということですか?
質問3:
2項目の「Model failed to converge: degenerate Hessian with 1 negative eigenvalues」の意味をご教示お願いします。
質問4:
3項目の「モデルが固有値に収束しない、、、」とは何を言っているのですか?
でもwarningのあとで結果が一応出ました。

#####
> temp lmer
Linear mixed model fit by REML [‘lmerModLmerTest’]
REML criterion at convergence: 7347086
Random_effects:
_Groups___Name________StdDev_Corr
_area_____(Intercept)_17852_______
__________before________104___100
_Residual_____________16293_______
Number_of_obs:_60__groups:__area___6
Fixed_Effects:
(Intercept)___area_b__area_c___area_d____area_e____area_f
___3821___-1123____1612____-1049____-2006____-2403
convergence code 0; 2 optimizer warnings; 0 lme4 warnings
#####
質問5:
Random_effects:のところで、before値のSdが104で、area_のが178というのは、意味として
「areaごとにbefore値が上下するけれど、それはareaによる影響ではないよ」
ということですか?
代わりにResidualが162なので、「ここに出ていない何かほかの因子が影響しているよ」
ということですか?
質問6:
Fixed Effectsのそれぞれの数値の意味をご教示お願いします。
よろしくお願いします。

>lmer((before+after_5+after_10)~area
??奇妙な入れ方ですね・・...アンサー

>lmer((before+after_5+after_10)~area
??奇妙な入れ方ですね・・
質問1の前ですが
>random slopeになるな
という点は,ランダム効果が,intercept の場合とslopeの場合で,モデルを比較し(AICを見る),どちらかを決めれば良いのです。
>従属変数を(before+after_5+after_10)
違います。
従属変数は時間(time)です。
time が3水準で,データは
time: 1, 2, 3
と入ります。この場合は,順序尺度にしてあるので
factor(time)
とします。
あるいは,before を0にして
time: 0, 5, 10
でもよい。
# 計算例
area<- rep(c("a", "b"), each=3) time<- rep(c(0, 5, 10), 2) y<- c( 534, 508, 157, 1046, 128, 173) dat<- dataframe(area, time, y) dat library(lmerTest) # ランダム傾きモデル mod1<- lmer(y ~ time + (time | area), data=dat) summary(mod1) AIC(mod1) # ランダム切片モデル mod2<- lmer(y ~ time + (1 | area), data=dat) summary(mod2) AIC(mod2) ## 終了 結果は, ランダム傾きモデル AIC 7370576 ランダム切片モデル AIC 6971192 ランダム切片モデルのほうが適合が良く time p値 00622 なので,時間変化に有意傾向ありです。 質問2~4は,変数がデータに比べて多すぎるために確定しないものと思われます。 以下の質問も,変数の入れ方がおかしいので,やり直してから考察してください。


質問タグ:random slope,In checkConv,derivs&quot,0 lme4 warnings,table all,ランダム切片モデル,data.frame

RWBYのvolume3の エメラルド、マーキュリー対 ココ、ヤツハシ戦で流れ…

RWBYのvolume3の エメラルド、マーキュリー対 ココ、ヤツハシ戦で流れ…

RWBYのvolume3の エメラルド、マーキュリー対
ココ、ヤツハシ戦で流れていた曲名がわかりません教えて下さい

アンサー

I’m the one


質問タグ:volume3,マーキュリー,エメラルド,RWBY,ヤツハシ戦,I’m the one

Ricottaの『ワルキューレロマンツェ』ってエロゲーは恋愛ゲームですか…

Ricottaの『ワルキューレロマンツェ』ってエロゲーは恋愛ゲームですか…

Ricottaの『ワルキューレロマンツェ』ってエロゲーは恋愛ゲームですか? 絵が好みなので気になってるんですが、恋愛を目的にしたようなゲームがあまり好きじゃないので迷ってます。
何か他にテーマがあって、そこにヒロインとの恋愛も絡んでくるとか、そういうのなら大丈夫なんですけど。

アンサー

完全な恋愛目的ではないですかね。
ヒロインが騎士、主人公がその騎士を補佐するヴェグライターという立場でタッグを組んで学園の頂点を目指すのが目的でそれを目指す過程で恋愛関係になる感じです。
とは言ってもイチャイチャするシーンも結構あるので半々って感じですかね。


質問タグ:ワルキューレロマンツェ,エロゲー,Ricotta,ヒロイン,恋愛,ヴェグライター,騎士

Republicans who remain silent in the face of Donald Trump’s ugly,…

Republicans who remain silent in the face of Donald Trump’s ugly,…

Republicans who remain silent in the face of Donald Trump’s ugly,
divisive tactics and attacks, take note of what another Republican said 150 years ago:
“To sin by silence when they should protest makes cowards of men”
Abraham Lincoln
日本語ではどういう意味でしょうか?

アンサー

ドナルド・トランプの醜く、分断を起こさせるような戦術と攻撃に直面して、沈黙を守り続ける共和党員たちは、150年前、もう1人の共和党員が言ったことに注目すべきだ。「抗議の声を上げるべき時に、沈黙で罪を犯すことは、人々から臆病者を作り出す」
※難しいですね、参考にしてください。


質問タグ:エイブラハム・リンカーン,日本語,共和党員,沈黙,臆病者,ドナルド・トランプ,共和党員たち

ReLIFEやシュタインズゲートみたいな現実ではありえないけどバトル系…

ReLIFEやシュタインズゲートみたいな現実ではありえないけどバトル系…

ReLIFEやシュタインズゲートみたいな現実ではありえないけどバトル系ではないアニメを教えてください。リゼロや僕だけがいない街などの有名な作品は見尽くしているので、そこまで有名じゃないや つも教えてください

アンサー

現代劇であることに絞らないとわけわかんなくなりますね。
アベノ橋魔法商店街
エイリアン9
乱歩奇譚
輪るピングドラム
地獄少女
四畳半神話大系
残響のテロル
凪のあすから
物語シリーズ
あたりで。


質問タグ:シュタインズゲート,リゼロ,リライフ,ピングドラム,地獄少女,乱歩奇譚,バトル系

refrainという英単語は (…を)控える、断つ、やめる、我慢する という…

refrainという英単語は (…を)控える、断つ、やめる、我慢する という…

refrainという英単語は (…を)控える、断つ、やめる、我慢する という意味があるようですが、 リトルバスターズのrefrainはどのような意味からrefrainという単語が使われているのでしょうか?

アンサー

確かに動詞としてのrefrainはそのような意味ですが、
名詞としてのrefrainは”繰り返し”という意味や(歌の)”サビ”という意味があります。
おそらくリトバスでは前者の”繰り返し”の意味でのrefrainでしょう。


質問タグ:Refrain,リトルバスターズ,英単語,意味,動詞,名詞

RDG、八雲みたいな恋愛要素がある小説を探しています! レッドデータ…

RDG、八雲みたいな恋愛要素がある小説を探しています! レッドデータ…

RDG、八雲みたいな恋愛要素がある小説を探しています! レッドデータガールの相楽と泉子の関係が好きで、似たようなヒーローヒロインが出てくる小説はありませんか?
最初は仲が悪かったけど、だんだんと仲が良くなって恋愛に発展するのが好きです。
恋愛ミステリにおいても男女コンビが好きです。
10代?30代まででお願いします。
好みの作品
・RDG
・八雲
・東雲清一郎は書を書かない
・SMシリーズ
・女学生探偵シリーズ

アンサー

ファンタジーでもミステリーでもないのですが、瀬尾まい子さんの「幸福な食卓」。
ヒロインの佐和子と同学年の少年、大浦は最初は勉強のライバルでしたが、高校に入ると、恋愛関係になり、つきあうようになります。ただ、ハッピーエンドではないので、泣いてしまうかも。


質問タグ:東雲清一郎,SMシリーズ,八雲,RDG,泉子,瀬尾まいこ,ハッピーエンド

RADWIMPSのOne Man Liveについて One Man Liveで使っているエフェクタ…

RADWIMPSのOne Man Liveについて One Man Liveで使っているエフェクタ…

RADWIMPSのOne Man Liveについて One Man Liveで使っているエフェクターはなんでしょうか?

アンサー

ギターのトーンを絞ってフェイザーかオートワウ+リヴァーブ+オーバードライブorブルースドライバーで歪ませます


質問タグ:One Man Live,エフェクター,RADWIMPS,フェイザー,オートワウ+リヴァーブ,トーン,ギター

R15は15歳でも買えるんですか?

R15は15歳でも買えるんですか?

R15は15歳でも買えるんですか?

アンサー

買えます


質問タグ:R15,15歳