🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

関数

関数(ファンクション・メソッド・サブルーチンとも呼ばれる)は、はプログラムのコードの一部であり、ある特定のタスクを処理するように設計されたものです。

最適化

最適化とはメソッドやデザインの最適な処理方法を選択することです。パフォーマンスの向上を目指す為に行われます。プログラミングにおける最適化は、アルゴリズムのスピードアップや、要求されるリソースを減らすことなどを指します。

Q&A

解決済

1回答

1132閲覧

Rのunirootの結果をoptimizeで同じ結果を算出したい

51sep

総合スコア22

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

関数

関数(ファンクション・メソッド・サブルーチンとも呼ばれる)は、はプログラムのコードの一部であり、ある特定のタスクを処理するように設計されたものです。

最適化

最適化とはメソッドやデザインの最適な処理方法を選択することです。パフォーマンスの向上を目指す為に行われます。プログラミングにおける最適化は、アルゴリズムのスピードアップや、要求されるリソースを減らすことなどを指します。

0グッド

0クリップ

投稿2019/10/11 20:15

unirootの結果をoptimizeで同じ結果を算出できるか確認したいと思っています。
unirootは0になる解を求め、Optimizeは最大化(または最小化)して最適解を求める関数ですが、
以下のunirootの0になる解(0.0036)をOpitimezeで制約をかけて算出することができるのでしょうか。
(以下のOptimizeでは‘6.924’の項があってもなくても、maximumは0.19999となる)

#uniroot uniroot(function(x) 2*log(x)+8*log(1-x)+6.924, c(0.01,0.1)) #$`root` #[1] 0.03638433 #$f.root #[1] 0.0002636439 #$iter #[1] 5 #$init.it #[1] NA #$estim.prec #[1] 6.103516e-05 #Optimize optimize(function(x) 2*log(x)+8*log(1-x)+6.924, c(0.01, 1), maximum=T) #$`maximum` #[1] 0.1999945 #$objective #[1] 1.919976

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

うまく動くかわかりませんが、1つの案としては、関数をabs()で囲って絶対値にしてしまった上でoptimizeで最小値を求める方法があります。
絶対値をとれば最小値=根(あれば)なので。

それではダメというのであれば、optimizeで根を探すというのはできないと思います。

投稿2019/10/11 23:21

gnbrganchan

総合スコア438

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

51sep

2019/10/12 00:02

できました!。gnbrganchanさん、ご回答をどうもありがとうございました。 以下に実行結果を記載しておきます。 おまけですが、nlm関数でもできました。 皆様にもお役に立ちますように。(注意:いずれの関数も初期値により算出値が変動します) optimize(function(x) abs(2*log(x)+8*log(1-x)+6.924), c(0.01, 0.1)) #$`minimum` #[1] 0.03637479 #$objective #[1] 0.0001814986 nlm(function(x)abs((2*log(x)+8*log(1-x)+6.924)), 0.1) #$`minimum` #[1] 5.91723e-07 #$estimate #[1] 0.03637866 #$gradient #[1] 45.49106 #$code #[1] 2 #$iterations #[1] 9
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問