[quote]Ariel8732 发表于 2016-4-28 20:10

点击<>这个符号，然后把你的代码复制进来你就知道了。

[code]
学习一下

[code]

[code]
set.seed(1234) #设置种子点
p <- runif(10,min=0,max=0.1)
p <- sort(p)
p_adj0 <- p.adjust(p,method = 'bonferroni')
p_adj1 <- p.adjust(p,method = 'holm')
p_adj2 <- p.adjust(p,method = 'BH')
p0 = length(p)*p
p1 = (length(p) - rank(p) + 1) * p
p2 = length(p) * p / rank(p)
pbind = cbind(p,p0,p_adj0,p1,p_adj1,p2,p_adj2)
pbind

[code]
set.seed(1234) #设置种子点
p <- runif(10,min=0,max=0.1)
p <- sort(p)
p_adj0 <- p.adjust(p,method = 'bonferroni')
p_adj1 <- p.adjust(p,method = 'holm')
p_adj2 <- p.adjust(p,method = 'BH')
p0 = length(p)*p
p1 = (length(p) - rank(p) + 1) * p
p2 = length(p) * p / rank(p)
pbind = cbind(p,p0,p_adj0,p1,p_adj1,p2,p_adj2)
pbind
[/嗷code]

[code]
*1、生成二次项
gen SMconerball=Mconerball*Mconerball
*2、线性回归
reg Mdf  Mconerball SMconerball  Msms Mszl Mqj Mqjcgl Mfqdq Mqcryq Myw Myellowcard Mredcard ,r //
*3、utest
utest Mconerball SMconerball  

*2、ologit回归
ologit Mdf  Mconerball SMconerball  Msms Mszl Mqj Mqjcgl Mfqdq Mqcryq Myw Myellowcard Mredcard ,r //
*3、utest
utest Mconerball SMconerball

*3、utest
utest Mconerball SMconerball ,prefix(Mdf)
[code]

［code
这样？

/code］

[code]
这样？
［/code］

复制代码

复制代码