R言語でZipfの法則
R言語でジップの法則のヒストグラムを描く。
R help archive:Zipf random number generationより引用。
x0<-1:10000; t<-0.999; p<-((t^x0)/x0)/(-log(1-t)) Y<-sample(x0,5000,replace=TRUE,prob=p) hist(Y,breaks=100)
R言語でジップの法則のヒストグラムを描く。
R help archive:Zipf random number generationより引用。
x0<-1:10000; t<-0.999; p<-((t^x0)/x0)/(-log(1-t)) Y<-sample(x0,5000,replace=TRUE,prob=p) hist(Y,breaks=100)