瓜生真也(横浜国立大学大学院環境情報学府)
June 30, 2014 @ Yokohama.R#1
多くの研究分野で使用される統計解析環境ソフト
1 + 3
[1] 4
a <- c(1, 2, 3)
a
[1] 1 2 3
b <- c(4, 5, 6)
a + b
[1] 5 7 9
a * b
[1] 4 10 18
mean(b)
[1] 5
10 * pi
[1] 31.42
三種類のあやめの品種のそれぞれからの50の花の蕚(がく)片の長さと幅、花弁の長さと幅
data(iris)
head(iris)
Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1 5.1 3.5 1.4 0.2 setosa
2 4.9 3.0 1.4 0.2 setosa
3 4.7 3.2 1.3 0.2 setosa
4 4.6 3.1 1.5 0.2 setosa
5 5.0 3.6 1.4 0.2 setosa
6 5.4 3.9 1.7 0.4 setosa
plot(
iris$Sepal.Length, iris$Petal.Length,
col=c("royal blue", "green", "orange")[iris$Species],
pch=16, cex=2)
一元配置分散分析… oneway.test()
3種の萼片の長さの平均値に差があるか
oneway.test(iris$Sepal.Length ~ iris$Specie)
One-way analysis of means (not assuming equal variances)
data: iris$Sepal.Length and iris$Specie
F = 138.9, num df = 2.00, denom df = 92.21, p-value < 2.2e-16
P < 0.001で萼片の長さの平均値には3種間で差がある
install.package("lattice", dependencies=TRUE)
library(lattice)
xyplot(
Petal.Length ~ Sepal.Length | Species,
data=iris)
Rの統合開発環境… もっと便利にRを使える
RStudioをインストールしたら…