5

x、a、 の3 つのパラメーターに依存するべき乗関数を描画したいと思いますgamma。関数は次のようになります。

powerlaw <- function(x, a, gamma){
   a*(x**(-gamma))
}

これをプロットしたいのですが、Rに選択した範囲を使用するように指示aしながら、指定する方法がわかりません。私はこれを試しました：gammax

require(ggplot2)
qplot(c(1,10), stat="function", fun=powerlaw(x, a=1, gamma=1), geom="line")

しかし、それは言います

Error in (x^(-gamma)): x is missing

もちろん、次のコードは と を修正することで機能aしgammaます。

powerlaw1 <- function(x){
   1*(x**(-1))
}
qplot(c(1,10), stat="function", fun=powerlaw1, geom="line")

何か案は？

4

2 に答える 2

3

引数を個別に指定する必要があります。

qplot(x=c(1,10), stat="function", 
      fun=powerlaw, geom="line", 
      arg=list(a=1, gamma=1))

詳細については、を参照?stat_functionしてください。

于 2012-11-01T10:32:48.750 に答える
2

data.frame適切なを返す関数を作成するだけですggplot2：

power_data = function(x, a, gamma) {
   return(data.frame(x = x, y = a * (x**(-gamma))))
}

> power_data(1:10, 1, 1)                                           
    x         y                                                    
1   1 1.0000000                                                    
2   2 0.5000000                                                    
3   3 0.3333333                                                    
4   4 0.2500000                                                    
5   5 0.2000000                                                    
6   6 0.1666667                                                    
7   7 0.1428571                                                    
8   8 0.1250000                                                    
9   9 0.1111111
10 10 0.1000000

プロットを作成します（より狭い間隔のxシリーズを使用して、より滑らかな線を取得することに注意してください）：

dat = power_data(seq(1,10,0.01), 1, 1)
qplot(dat$x, dat$y, geom = "line")

ここに画像の説明を入力

于 2012-11-01T10:37:44.333 に答える