×
[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
# psychパッケージのdescribe.by関数でもいいが…
# グループごとの記述統計
aggregate(iris[1:4], list(iris[,5]), mean, na.rm=TRUE) # 平均
## これでも同じ aggregate(iris[1:4], iris[5], mean, na.rm = T)
aggregate(iris[1:4], list(iris[,5]), median, na.rm=TRUE) # 中央値
aggregate(iris[1:4], list(iris[,5]), sd, na.rm=TRUE) # 標準偏差
aggregate(iris[1:4], list(iris[,5]), max, na.rm=TRUE) # 最大値
aggregate(iris[1:4], list(iris[,5]), min, na.rm=TRUE) # 最小値
aggregate(iris[1:4], list(iris[,5]), length) #データ数、NAなし
aggregate(iris[1:4], list(iris[,5]), function(x) length(na.omit(x))) # データ数、NA抜き
aggregate(iris[1:4], list(iris[,5]), function(x) sum(is.na(x))) # データ数、NAのみ
# 追加
by(iris[1:4], list(a = iris[,5]), summary)
library(psych)
describe(iris)
by(iris[1:4], list(a = iris[,5]), describe)
describe.by(iris[1:4], iris[5])
# グループごとの記述統計
aggregate(iris[1:4], list(iris[,5]), mean, na.rm=TRUE) # 平均
## これでも同じ aggregate(iris[1:4], iris[5], mean, na.rm = T)
aggregate(iris[1:4], list(iris[,5]), median, na.rm=TRUE) # 中央値
aggregate(iris[1:4], list(iris[,5]), sd, na.rm=TRUE) # 標準偏差
aggregate(iris[1:4], list(iris[,5]), max, na.rm=TRUE) # 最大値
aggregate(iris[1:4], list(iris[,5]), min, na.rm=TRUE) # 最小値
aggregate(iris[1:4], list(iris[,5]), length) #データ数、NAなし
aggregate(iris[1:4], list(iris[,5]), function(x) length(na.omit(x))) # データ数、NA抜き
aggregate(iris[1:4], list(iris[,5]), function(x) sum(is.na(x))) # データ数、NAのみ
# 追加
by(iris[1:4], list(a = iris[,5]), summary)
library(psych)
describe(iris)
by(iris[1:4], list(a = iris[,5]), describe)
describe.by(iris[1:4], iris[5])
PR
Comment
Trackback
Trackback URL
Comment form
