w<-t(matrix(b[1:6],2,3)) rho<-b[7] X<-cbind(1,seq(1, 20, 1),min(birthY)) U<-cbind(X%*%w,0) ival<-log(exp(U[,1]/rho)+exp(U[,2]/rho)) #クラスター選択確率 CL1<-exp(rho*ival)/(exp(rho*ival)+exp(U[,3])) CL2<-exp(U[,3])/(exp(rho*ival)+exp(U[,3])…

person-year型の架空データ(その1)をベースに、 仮定するデータ発生メカニズムを変更就業選択に階層構造を仮定する ①就業するか否か ②就業するとしたら、正規雇用か、非正規雇用か ＊単純化のために、ここでは就業形態を正規、非正規の2つにしているその他の…

パーソン・イヤー型の架空データをRを使用し作成架空データなので、 とりあえずざっくりと以下のような状況を仮定してみる再就業の分析 ・正規雇用からの離職を20-29歳の間に経験した女性が分析対象 ・分析対象となるイベントは「再就業」 ・就業状態は「正…