R语言生成随机数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助大神：生成混合偏正态分布的随机数。
英文参考里有程序但是自己不太明白，而且没有运行结果。
下面是程序：

### Example 2

## Do not run

# set.seed(150520)

# data(enzyme)

# x <- enzyme

# Enzyme1.out <- MixNRMI1(x, Alpha = 1, Beta = 0.007, Gama = 0.5,

# distr.k = 2, distr.p0 = 2, asigma = 1, bsigma = 1, Meps=0.005,

# Nit = 5000, Pbi = 0.2)

# The output of this run is already loaded in the package

# To show results run the following

# Data

data(enzyme)

x <- enzyme

data(Enzyme1.out)

attach(Enzyme1.out)

# Plotting density estimate + 95% credible interval

m <- ncol(qx)

ymax <- max(qx[,m])

par(mfrow=c(1,1))

hist(x,probability=TRUE,breaks=20,col=grey(.9),ylim=c(0,ymax))

lines(xx,qx[,1],lwd=2)

lines(xx,qx[,2],lty=3,col=4)

lines(xx,qx[,m],lty=3,col=4)

# Plotting number of clusters

par(mfrow=c(2,1))

plot(R,type="l",main="Trace of R")

hist(R,breaks=min(R-0.5):max(R+0.5),probability=TRUE)

# Plotting sigma par(mfrow=c(2,1))

plot(S,type="l",main="Trace of sigma")

hist(S,nclass=20,probability=TRUE,main="Histogram of sigma")

# Plotting u

par(mfrow=c(2,1))

plot(U,type="l",main="Trace of U")

hist(U,nclass=20,probability=TRUE,main="Histogram of U")

# Plotting cpo

par(mfrow=c(2,1))

plot(cpo,main="Scatter plot of CPO's")

boxplot(cpo,horizontal=TRUE,main="Boxplot of CPO's")

print(paste('Average log(CPO)=',round(mean(log(cpo)),4)))

print(paste('Median log(CPO)=',round(median(log(cpo)),4)))

detach()

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：生成随机数 R语言 随机数 example ExamP 正态分布 英文 程序 而且

首先说哪里不明白，还有enzyme在哪个package里面，找到后library(pkg)

bbslover 发表于 2016-5-27 22:37
首先说哪里不明白，还有enzyme在哪个package里面，找到后library(pkg)

r语言我是不懂的，然后现在想生成一组混合偏正态的随机数，关于这个分布的文献里有这个程序生成随机数的，我想直接拿来用。其实你说的我也不懂，哭丧脸

enzyme应该是某个包里自带的数据集吧，楼主百度一下。另，最好说一下你看的是哪本文献。

enzyme应该是某个包里自带的数据集吧，楼主百度一下。另，最好说一下你看的是哪本文献。

install.packages("BNPdensity")
library(BNPdensity)

之后，运行你的代码就可以了。

bbslover 发表于 2016-5-27 22:55
install.packages("BNPdensity")
library(BNPdensity)

library(BNPdensity)
之后跳出CRAN mirror 该选择哪一个呢有好多选择

bbslover 发表于 2016-5-27 22:55
install.packages("BNPdensity")
library(BNPdensity)

library(BNPdensity)
之后跳出CRAN mirror 该选择哪一个呢有好多选择

只R初见 发表于 2016-5-27 22:55
enzyme应该是某个包里自带的数据集吧，楼主百度一下。另，最好说一下你看的是哪本文献。

BNPdensity   就是这个的pdf文件 在这里找的程序

bbslover 发表于 2016-5-27 22:55
install.packages("BNPdensity")
library(BNPdensity)

请问程序包应该怎么放，

> install.packages("BNPdensity")
install.packages("BNPdensity")里有警告:
  'lib="C:/PROGRA~2/R/R-211~1.1/library"'不可写
--- 在此連線階段时请选用CRAN的鏡子 ---
错误于contrib.url(repos, type) : 试试在不设定镜像的情况下用CRAN
>  library(BNPdensity)
错误于library(BNPdensity) : 不存在叫'BNPdensity'这个名字的程辑包