iterGlm2nt {mseq} | R Documentation |
Fit the Poisson linear model by iteratively fitting glm and gene expression levels.
iterGlm2nt(data, thrd = 0.01, max_iter = 10)
data |
the data frame generated by expData2nt.R |
thrd |
the threshould of improvement of R squared |
max_iter |
maximum number of iterations |
an glm structure. Summarize it will list the coefficients and their p-values.
data(g1_part) # for real data, please use read.csv, like g1 <- read.csv("g1.csv") data <- expData2nt(g1_part, 2, 1) data.glm <- iterGlm2nt(data) summary(data.glm)