0

私はRにあまり詳しくありません。Rを使用して、さまざまなラムダ（1から10）のポアソン分布プロットを作成し、それぞれのプロットを比較として表示していました。

しかし、私はタイトルを追加したいと思います：ラムダに従って自動的にグラフにプロット1の場合は「ラムダ= 1」、プロット2の場合は「ラムダ=2」など。しかし、タイトルを自動的に更新する方法がわかりませんでした。これは私のコードです。10個の異なるグラフを正しく出力できましたが、対応するラムダをタイトルに自動的に更新または追加する方法がわかりません。誰かが私にいくつかのヒントを教えてもらえますか？

また、プロット1〜5のフォントサイズを「小さい」、次にフォントサイズを6〜10と言うこともできますか？

ありがとう

the_data_frame<-data.frame(matrix(ncol=10,nrow=21))
lam<-seq(1,10,1)
lam
x<-seq(0,20,1)
x
for (i in 1:10){
the_data_frame[i]<-exp(-lam[i])*lam[i]**x/gamma(x+1)
}
the_data_frame<-cbind(the_data_frame, x)


par(mfrow=c(5,2))
for (i in 1:10){
plot(the_data_frame[[i]]~the_data_frame[[11]], the_data_frame)
}
4

4 に答える 4

1

問題を単純化できます。ラムダ値に対して 1 つのループを使用して、各反復で毒の式を使用して y の値を計算し、それをプロットします。引数を使用mainして、各プロットのタイトルを追加します。ここではbquote、ラムダ値のプロットマス形式を取得するために使用しています。

たとえば、 lambda の 4 つの値の場合、次のようになります。

x<-seq(0,20,1);lam = c(0.5,1,2,4)
par(mfrow=c(2,2))
lapply(lam,function(lamd){
  y <- exp(-lamd)*lamd*x/gamma(x+1)
  plot(x,y,main=bquote(paste(lambda,'=',.(lamd))),type='l')

})

ここに画像の説明を入力

于 2013-03-14T19:45:45.737 に答える
0

テキストを入力して図のタイトルを更新できるインタラクティブなウィンドウを探していましたか？はいの場合は、tcltkパッケージを探すことをお勧めします。

http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/modalDialog.htmlを参照してください

于 2013-03-16T19:13:28.670 に答える
0

これは役立つかもしれません：

for (i in 1:10){
    plot(the_data_frame[[i]]~the_data_frame[[11]], the_data_frame,
        main=paste("lambda=", i, sep=""))
}
于 2013-03-14T19:22:53.383 に答える
0 
library(ggplot2)
xval <- rep(0:20,10)
lambda <- rep(1:10,21)
yvtal <- exp(- lambda)*lambda**xval/gamma(xval+1)
the_new_data_frame <- data.frame(cbind(xval,lambda,yval))

plot1 <- ggplot(the_new_data_frame, aes(xval, yval)) + geom_line(aes(colour=factor(lambda)))
plot1
plot1 + facet_grid(~lambda)
于 2013-03-14T20:58:38.830 に答える