問題タブ [rpy2]

回帰には rpy2 を使用しています。返されたオブジェクトには、係数、残差、適合値、適合モデルのランクなどを含むリストがあります)。

ただし、フィット オブジェクトで標準誤差 (または R^2) を見つけることができません。R で lm を直接モデルを実行すると、summary コマンドで標準エラーが表示されますが、モデルのデータ フレームで直接アクセスできません。

rpy2を使用してこの情報を抽出するにはどうすればよいですか?

サンプルのpythonコードは

Rを使用してベクトルのヒストグラムを作成しようとしています。これをプロットすると、x軸は5から9になりますが、ヒストグラムの一部は5より前（最小数は4.414002）で、一部は9より後です。同じことが私のy軸にも当てはまり、0から5000になりますが、タレスの「ステーブ」はこの5000より少し上になります。

ベクトルの最小値から開始しないのはなぜですか？念のために言っておきますが、長さが足りないのはバーだけで、グラフは適切なサイズです（したがって、marまたはomaの設定に問題はありません。変更してみましたが、役に立ちませんでした）

以下は私のコードです

編集：

私は問題が何であるかを知り、Rコードでは次のようになります。

しかし、範囲は4〜9の間であるため（この例ではそれより少ない）、20回の休憩をとるのはあまりにも多すぎました。6に休憩を入れると、問題が解決します。

このRスクリプトと同等のことをしたい：

rpy2で。ただし、rpy2.robjects を R としてインポートすると、Rr['!duplicated'] が認識されません (次のように)。

rpy2で!duplicatedを使用するにはどうすればよいですか?

編集：

が動くので作り方募集中！rpy2で作業

トレースバック（コピー貼り付け）は次のとおりです。

どんな手掛かり？

編集：これがコードです

誤解しないでください、それはエラーが発生している行です。別のPythonプロンプトを試してみましたが、import rpy失敗しません。ウェブサイトでのみ失敗しています。なんで？

すべてが開発サーバーで正常に動作します。ただし、httpd と wsgi を実行している運用サーバーにプロジェクトを移動すると、次のエラーが表示されました。

何が起こっているのか分かりますか？ありがとう。

更新 Apache エラー ログまで追跡すると、「バイトコードの不一致」のようなものが見つかりました。GoogleはそれがRに関連していることを発見しました。そのため、Rを再コンパイルして問題を解決しました。

ここにあるRの例と同じことをしたい: http://www.r-tutor.com/r-introduction/data-frame/data-frame-column-slice。彼は、次のような mtcars データベースを利用しています。

次に、インデックス ベクトルを使用して mpg および hp 列を選択します。

この場合、mpg を取得するために、rpy2 で 1 行を取得するのは簡単です。

ただし、2つの列を選択する方法がわかりません

与える

R.vectors にある他のすべてのベクトルと同じように (私はそれらすべてを試しました)

私の質問は、どうすれば rpy2 でインデックス ベクトルを取得できますか、それが存在しない場合、どうすれば 2 つの列を選択できますか?

ありがとう！

私は rpy2 と numpy の新しいユーザーです。RHEL5でR 2.14.1、python 2.7.1、rpy2 2.2.5、numpy 1.5.1を使用しています。

データを numpy 配列に読み込み、rpy2 関数を使用する必要があります。ただし、rpy2 で配列を使用する前に、欠損値をマスクする必要があります。

値のマスキングに問題はありませんが、結果のマスクされた配列で rpy2 を動作させることはできません。マスクされた numpy 配列では numpy2ri 変換が機能しないように見えますか? (以下のエラーを参照)

どうすればこれを機能させることができますか？マスクされた値を無視するように rpy2 に指示することは可能ですか? 後でより高度な統計を行うため、scipy/numpy を直接使用するのではなく、R に固執したいと思います。

ありがとう。

更新: rpy2 はマスクされた numpy 配列を処理できないため、-99 値を numpy NaN 値に変換してみました。どうやら rpy2 は numpy NaN 値を R スタイルの NA 値として認識します。

以下のコードが機能するのは、r.sd() 呼び出しで、rpy2 に NA 値を使用しないように指示できるためです。しかし、最初の NaN 置換は、numpy マスクを適用するよりも明らかに遅くなります。

大規模なnumpy ndarray全体で-99からNaNへの置換を行うより高速な方法をPythonウィザードで教えてもらえますか? または、別のアプローチを提案しますか？

ありがとう。

このRコードと同等のことをしたい：

私がする時

このエラーが発生します：

私は試した

と

どちらも私が望む結果をもたらさない。では、colnamesを使用してデータフレームの列名を変更するにはどうすればよいですか？

最後の 3 つの prcomp() 呼び出しで、次のエラーが発生します。

rpy2 経由で prcomp() の NA 値を処理するために na.action (またはその他の手段) を適切に使用する方法がわかりません。単なる構文の問題であることを願っています。

助けてくれてありがとう。

R、Rpy2、および ggplot2 の既存の知識に妨げられずに、Python から単純なテーブルの散布図を作成することは決して好きではありません。

これをセットアップするために、インストールしたばかりです：

• Ubuntu 11.10 64 ビット
• R version 2.14.2(r-cran ミラーより)
• ggplot2(経由R> install.packages('ggplot2'))
• rpy2-2.2.5(経由easy_install)

これに続いて、ggplot2 を使用してインタラクティブな R セッションからいくつかのサンプル データフレームをプロットできます。

ggplot2ただし、オンラインで見つけた例で見たように単にインポートしようとすると、次のエラーが発生します。

誰が私が間違っているのか教えてもらえますか? 私が言ったように、問題のあるインポートはオンラインの例から来ているので、rpy2を介してgplot2を使用する必要がある他の方法があるかもしれません.

参考までに、上記の問題とは関係ありませんが、インポートが機能するようになったら、プロットしたいデータフレームの例を次に示します (例を見ても問題はありません)。アイデアは、x 軸に長さ、Y 軸にパーセンテージを示す散布図を作成することです。ブール値を使用してドットに色を付け、ファイル (画像または pdf) に保存します。これらの要件は非常に限られているため、代替ソリューションも歓迎されます。