22:48いつも拝見させていただいています。ネット検索などをしましたが,見つけられなかったので質問させていただきました。在統計ソフトRで,glmmML()を使い,モデル選択を行いたく思っています。logLik()を使うと有効なメソッドがないとメッセージが出て,計算することができません。(glm()では問題なく計算できました)。R上で,glmmMLで得られた結果に対し,logLikをどのように計算すれば良いのでしょうか?ご存知の方がいればお教え願えないでしょうか?文献やネット上で見つけられず困っております。23:23追試できるプログラム(必要ならデータも)を呈示しておくとよいと思うのですけど。現象が再現できる最短のプログラム最小のデータ。あなたがやったやり方でよいのかどうか,あなたの書いた文章からはわからないですよねぇ。倍がデビアンスだったかと)私はこの分析手法については知りませんので,それであなたが求めたいと思っているものが得られるかどうかわかりません。やるというのでは駄目なんですか。わかってなくてすみません。23:54すみませんでした。データと関数を記しました。id)というデータが50サンプル分あります。Nは子の数の上限で,yは1ペア当たりの子の数,xは体サイズ,sは性,idは個体番号になります。体サイズや性を使って,子の数を説明できるか?個体差はどのくらいあるのかが知りたいことです。これでご判断いただけるでしょうか?ご意見いただけたらうれしいです。00:00申し訳ありません。途中で送信してしまいました。書くためのアルゴリズムを知りたいと言うことなんでしょうか?はい,そうなんです。手計算ではなく,その計算をR上でする方法を探しています。言葉足らずですみませんでした。どのようなアドバイスでもありがたく思います。00:038273ですが,ちょっとやってみては書き足していたので,まだあなたが読んでいない部分ができているかも。00:20ありがとうございます。devianceは出してくれるのですが,devianceは出ないのかと思っていました。一度やりなおしてみます。