-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy pathhlaModelFromObj.Rd
61 lines (53 loc) · 1.49 KB
/
hlaModelFromObj.Rd
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
\name{hlaModelFromObj}
\alias{hlaModelFromObj}
\alias{hlaModelToObj}
\title{
Conversion between the in-memory model and the object that can be
saved in a file
}
\description{
Build a model \code{\link{hlaAttrBagClass}} from an object of
\code{\link{hlaAttrBagObj}} which is stored in an R object file, or convert
\code{\link{hlaAttrBagClass}} to \code{\link{hlaAttrBagObj}}.
}
\usage{
hlaModelFromObj(obj)
hlaModelToObj(model)
}
\arguments{
\item{obj}{an object of \code{\link{hlaAttrBagObj}}}
\item{model}{an object of \code{\link{hlaAttrBagClass}}}
}
\value{
\code{hlaModelFromObj} returns \code{hlaAttrBagClass}, and
\code{hlaModelToObj} returns \code{hlaAttrBagObj}.
}
\author{Xiuwen Zheng}
\seealso{
\code{\link{hlaAttrBagging}}
}
\examples{
# make a "hlaAlleleClass" object
hla.id <- "DQB1"
hla <- hlaAllele(HLA_Type_Table$sample.id,
H1 = HLA_Type_Table[, paste(hla.id, ".1", sep="")],
H2 = HLA_Type_Table[, paste(hla.id, ".2", sep="")],
locus=hla.id, assembly="hg19")
# training genotypes
region <- 100 # kb
snpid <- hlaFlankingSNP(HapMap_CEU_Geno$snp.id, HapMap_CEU_Geno$snp.position,
hla.id, region*1000, assembly="hg19")
train.geno <- hlaGenoSubset(HapMap_CEU_Geno,
snp.sel = match(snpid, HapMap_CEU_Geno$snp.id))
# train a HIBAG model
set.seed(1000)
model <- hlaAttrBagging(hla, train.geno, nclassifier=2)
print(model)
mobj <- hlaModelToObj(model)
is(model)
is(mobj)
# close the HIBAG model explicitly
hlaClose(model)
}
\keyword{HLA}
\keyword{genetics}