読者です 読者をやめる 読者になる 読者になる

Rを通じて統計学を学ぶ備忘録ブログ

SPSSからRに移行したい私のような人向けのR解説ブログ兼学習用備忘録。

ちょっとした備忘録のページ。

その他(R関係)
options(scipen = 5)

##想定クリック数1000回、クリック率1%の時の二項分布。
d <- as.data.frame(dbinom(x = 0:30, size = 1000, prob = 0.01))
d$x <- 0:30
colnames(d) <- c("y","x")
ggplot(data = d, aes(x = factor(x), y = y))+ geom_bar(stat = "identity")

##平均λ=1.43のポアソン分布。7時間で10cvの場合。
d <- as.data.frame(dpois(x = 0:10, lambda = 1.428))
d$x <- 0:10
colnames(d) <- c("y","x")
ggplot(data = d, aes(x = factor(x), y = y))+ geom_bar(stat = "identity")

参照元http://www.mm-lab.jp/analysis/data-analysis-start-with-excel-5th/

【オンライン広告の認知効果】

x <- 1:1000
y <- 5.980469*x^0.265342 + rnorm(x, 0, 5)
d <- data.frame(x,y)
plot(x,y)
fit <- lm(log(y) ~ log(x))
summary(fit)

result.con <- data.frame(x = d$x, conf = exp(predict(fit, interval="confidence", level = 0.90)))

ggplot(d, aes(x = x, y = y)) + geom_point() +
  geom_ribbon(data = result.con, aes(ymin=result.con$conf.lwr, ymax = result.con$conf.upr), fill = "darkgreen", alpha=0.2) +
  geom_line(data = result.con, aes(x = x, y = conf.fit)) 

参照元http://www.videor.co.jp/press/2015/150903.pdf

広告を非表示にする