問題タブ [nlm]

関数 F1 が与えられます。

F1 は最後の式の結果を返します

最後の方程式の結果がゼロに収束するように、パラメーター x の値を決定したいと考えています。と

結果は得られますが、満足のいくものではありません:

nlm コマンドの引数（gradtol、stepmax、iterlim など）でいろいろいじりましたが、良い結果が得られませんでした。optim、optimize、uniroot も試しましたが、まったく実行できませんでした。

u と x は負の指数関数的展開を示します。x が減少すると、u は指数関数的に増加します。x がゼロの場合、u は有限値になります。x にも上限がありますが、これは不明です。したがって、反復が下限（ゼロ）から始まり、段階的に増加する場合は有望だと思いました。しかし、stepmax の値を減らしても増やしても、結果は良くなりません。

r-communityからのヒントをいただければ幸いです。

どうもありがとうございました。

PS: matlab では、同僚が fsolve(@(x) F1(x,u_target,C1,C2,...),0) を使用しており、正常に動作します。

GenABEL に実装された非加法的 GWAS モデルを使用しています。ドミナントおよびリセッシブ GWAS 非加法的モデルの VIFGC 関数は正常に動作していますが、同じデータを使用した過剰支配継承では、VIFGC_ovedom 関数がエラーを出しています: Error: Error in nlm(GC_VIF_nlm, c(F = F1, K = K1)) :パラメータに値がありません

また、これらの非相加モデルの重要なマーカーによって説明される表現型の分散を計算する方法。

ありがとう、