続いて、本来の目的である個人ベースの言及数の解析に入ります。こんな感じのスクリプトで解析してみることにします。 data.stat <- function(d, pickup=seq_len(ncol(d))){ d <- d[,pickup] lm.ba <- apply(d, 1, function(v) lm(as.numeric(v) ~ pickup)[[…