Мне нужно написать собственную функцию, чтобы нарисовать функцию плотности биномиального распределения и, следовательно, провести соответствующий график при n = 20 и p = 0,1,0.2, ..., 0,9. Также мне нужны комментарии на графиках.R, график биномиального распределения
Я пробовал это;
graph <- function(n,p){
x <- dbinom(0:n,size=n,prob=p)
return(barplot(x,names.arg=0:n))
}
graph(20,0.1)
graph(20,0.2)
graph(20,0.3)
graph(20,0.4)
graph(20,0.5)
graph(20,0.6)
graph(20,0.7)
graph(20,0.8)
graph(20,0.9)
#OR
graph(20,scan())
Мой первый вопрос: есть ли способ, так что мне не нужно, чтобы записать несколько раз линия graph(20,p)
за исключением использования scan()
?
Мой второй вопрос:
Я хочу, чтобы увидеть график в одном устройстве или хотите ударить ENTER
, чтобы увидеть следующий график. Я написал
par(mfcol=c(2,5))
graph(20,0.1)
graph(20,0.2)
graph(20,0.3)
graph(20,0.4)
graph(20,0.5)
graph(20,0.6)
graph(20,0.7)
graph(20,0.8)
graph(20,0.9)
, но график слишком крошечный. Как я могу красиво представить графики с помощью строки заголовка n = 20 и p = значение, которое я использовал для рисования графика? [Хотя это можно сделать, написав mtext()
после вызова функции , но при этом мне нужно написать аналогичную линии несколько раз. Поэтому я хочу сделать это, включая функцию . ]
Мой последний вопрос:
О комментариях. Графики показывают, что, как вероятность успеха, p увеличивается, график стремится вправо, т. Е. График искажается справа.
Есть ли способ прокомментировать график, используя program
?
Рассмотрите цикл или 'lapply', чтобы пройти через несколько значений p. – Thomas