hisemi-package: Hierarchical semiparametric regression model to a large...

Description Details Author(s) References See Also Examples

Description

The package includes functions for fitting hierarchical semiparametric regression model to a large number of parametric test statistics

Details

Package: hisemi
Type: Package
Version: 1.1-0
Date: 2017-07-09
License: GPL version 2 or newer
LazyLoad: yes

Author(s)

Long Qu rtistician@gmail.com

Maintainer: Long Qu rtistician@gmail.com

References

Long Qu, Dan Nettleton, Jack Dekkers (2012). A hierarchical semiparametric model for incorporating inter-gene relationship information for analysis of genomic data. Biometrics, 68(4):1168-1177

See Also

pi0-package

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
## simulate some fake data
G=100  ## for demonstration only. Normally, G should be much larger
sdncp=1.3
n1=n2=5
df=n1+n2-2
set.seed(54457704)
x=runif(G,1,G)
f=function(x)sin(x*pi/1000)+1
Pi.i=1/(1+exp(f(x)))
Z.i=rbinom(G,1,1-Pi.i)
t0.i=rt(G,df)
ncp.i=rnorm(G,0,sdncp)
t1.i=rt(G,df, ncp.i)
t.i=ifelse(Z.i==0,t0.i,t1.i)

## fit model
(plfit=penLik.EMNewton(t.i, x, df, spar=10^seq(0,8,length=30),plotit=FALSE))
(plfit0=scaledTMix.null(t.i, df))

## Not run: 
plot(plfit)
plot(t.i, plfit$lfdr, pch='.')
lines(sort(t.i), plfit0$lfdr[order(t.i)], col=2, lwd=3)

## End(Not run)

gitlongor/hisemi documentation built on May 17, 2019, 5:28 a.m.