問題タブ [s]

R と S の違いは何ですか?

S-Plus から R を呼び出す良い方法について何か提案はありますか? 理想的には、R にコードを渡してデータを取得するだけで、それらを統合するために複雑すぎるものを書く必要はありません。

Omegahat のRinSパッケージに精通していることを付け加えておきますが、使用したことはありません。私は、Tibco が引き継ぐ前に、Insightful が環境を統合する努力をしていたという印象を受けました。

編集: RinS は Windows では動作しないことがわかりました。最も簡単な解決策は、Rscript を使用することであることがわかりました。これは S-Plus からsystem()コマンドで呼び出すことができます。たとえば、簡単なスクリプトを次に示します。

S-Plus から呼び出します。

次に、結果をテキスト ファイルに保存し、スクリプトの実行後に S-Plus にインポートします。

R / S / 非線形回帰に関連する問題があり、私は R プログラマーではないので、助けが必要です。

tt と td の 2 つの配列があります。

非線形関数の最小二乗和が最小になるように、パラメーター a、b、c を見つける必要があります。

これを行う方法がわかりません。nls()関数を試しnls2()ましたが、運がありませんでした...

前もって感謝します。

編集：

私のデータ:

以下の回答の方法を使用すると、ランダム データの値は問題ありませんが、使用しているデータはモデル メッセージの評価時に生成された欠損値または無限大を返します。

すぐにデータを提供できなくて申し訳ありません。

R / S-PLUSで、カテゴリ列ごとにグループ化された複数の要約統計量を1回のショットで取得しようとしています。いくつかの関数を見つけましたが、それらはすべて、のように、呼び出しごとに1つの統計を実行しますaggregate()。

私が探しているのは、平均、最小、最大、標準など、同じグループの複数の統計を1回の呼び出しで取得することですが、それは実行可能ですか？

ここに書き直された質問:

私はいくつかの進歩を遂げましたが、Rから奇妙な動作を得ています...

これが私が始めているxtsです

私が作りたいもの：

私はこのような関数を書きました：

これにより、thisDate が正常に追加されます (コマンドラインでコードを段階的に実行します)。しかし、nextDate を追加するビットはそれを上書きしているようです! また、予想外の NA の行が表示されるようです。まだこれに取り組んでいます...

最初の列に「タイトルなし」を入れて、実際にはベクトル/マトリックスの一部ではなく、xts 日付インデックスであることを示しました。

余分な行を削除するのは、上書きの問題をまだ解決しておらず、実験していたためです。それは最終的な答えにある必要はありませんが、私が現在行っているところです。

そして最後に、この結果を調べて nextDate を日付に変換しようとすると....

ということで、ちょっとモヤモヤ…。

以下の元の質問:

私はRに時系列を持っています(私は速く学んでいますが、明らかに十分に速くはありません!)

各行の新しい列の NEXT 行の日付インデックスを使用して、その派生物を作成したいと思います。

Value に対しては (Lag を使用して) 行うのは非常に簡単ですが、date-index 自体に対してはそうではありません。

おそらく、さまざまなルックアップなどを使用してそれを行う方法を理解できますが、面倒です。他のフィールドで一致させるか、「R に忠実」とは感じられない行番号をいじる必要があります。

それを行うための素晴らしく、きちんとした、エレガントな方法はありますか?

私は「D'OH！」に行くと確信しています。誰かが答えを出したらすぐに！しかし、これまでのところ、このサイトで日付インデックスの遅れについての答えを見つけていません。

私がこれをやりたい理由は、次に、日付の各ペアを連続して使用して、別のシリーズを調査したいからです。したがって、これを行うためのより良い方法があるかもしれません。

各要素の「ペソ」の値をテストする必要があります (以下の複製コードを参照)。因子が「ペソ」の全体合計の 50% に達するかどうか、各因子の値を新しいオブジェクト「結果」に貼り付ける必要があります。それ以外の場合、R は「ペソ」の集計値が最も低い因子を評価し、検討する必要があります。集計「ペソ」の次の列の係数をもう一度。基本的に、このプロセスは、スコアが最も低い因子を次の因子に置き換えます。要因が 50% のしきい値を超えるまで、このプロセスを繰り返す必要があります。だから私の質問は、どこから始めればいいですか?

次のようになるはずです

したがって、因子 8 はより大きなシェア 4/(4+2+3) = 0.4444444 を得ましたが、ラウンド a で 50% のしきい値に達しませんでした。したがって、さらに何かが必要です。集計を繰り返しますが、最初のラウンドで最も低い集計値を取得したため、列 'a' の係数 9 ではなく、列 'b' の係数 7 を考慮します。

このチュートリアルを使用して、ggplot2 を使用してヒートマップをプロットします

http://learnr.wordpress.com/2010/01/26/ggplot2-quick-heatmap-plotting/

今度は、同じ長さの onother ファクターを使用して垂直方向のファクターを注文したいと思います。

誰かがそれがどのように機能するか知っていますか?

ありがとう