2007-08-01から1ヶ月間の記事一覧

以下，修正・追記（9/12）

Rメモ

どうも，Multiple R-SquaredやAdjusted R-squaredが大きすぎる。 内容的には同じだが，若干違う定式化で推定してみる。 reg =1 summary(lm(regdata2[[reg]]$y ~ regdata2[[reg]]$X[,2:3])) ##出力 Call: lm(formula = regdata2[[reg]]$y ~ regdata2[[reg]]$X…

■

Rメモ

##################店舗別回帰分析 仕切り直し####### #bayesmを読み込む。 library(bayesm) #cheeseデータを読み込む。 data(cheese) ##回帰分析用のデータ retailer=levels(cheese$RETAILER) nreg=length(retailer) regdata2=NULL for (reg in 1:nreg) { n…

■

Rメモ

#再ロード load("regdata.Rout") #店舗別価格の記述統計の計算 regPriceSummary=matrix(rep(0,440),nrow=88) for (reg in 1:88){ regPriceSummary[reg,1] = mean(regdata[[reg]]$X[,2]) regPriceSummary[reg,2] = median(regdata[[reg]]$X[,2]) regPriceSumm…

■

Rメモ

#bayesmを読み込む。 library(bayesm) #cheeseデータを読み込む。 data(cheese) #定数項をXに含めないバージョン（変数は対数変換前） retailer=levels(cheese$RETAILER) nreg=length(retailer) regdata=NULL for (reg in 1:nreg) { #小売店舗の名前を保存 n…