n <- 1000 # ランダムに脱落する prob <- c(0.5,0.5) censor <- sample(0:1,n,replace=TRUE,prob=prob) # 指数乱数で死亡時刻が決まる #times <- cumsum(rexp(n)) times <- sort(rexp(n)) x <- ((n-1):0)/(n:1) v <- rep(1,length(x)) v[which(censor == 1)] …

「仮説Xと仮説Yとのどちらなのかを判定しようとして実験を観察していたら、XでもYでもなく仮説Zなのではないか、と思い至る」ことがあるように思うが、その思考過程とは、どういうことか、という話題が出ました。 そんなことを考えるための例。 今、1か2が出…

いくつかの基準カテゴリがあって、それらに遠近関係があるという。多次元だと視覚化しにくいので、２次元のしかも円周上に配置する。その上で、たくさんあるその他もろもろの点を、その基準カテゴリに帰属する「割合」を基に、その２次元円内にプロットする…