methylMnM (version 1.10.0)

pmultinom: Call C programs to R.

Description

Call C programs to R for calculate p-value of each bin with multinomial distribution.

Usage

pmultinom(T, SIZE,length, P1, P2, P3, P4, C1, C2, pvalue=rep(0,length(T)))

Arguments

T

Statistic.

SIZE

The sum of MeDIP-seq count and MRE-seq count of each bin of two samples.

length

The number of bins.

P1

The probability of MeDIP-seq of sample 1 in multinomial distribution.

P2

The probability of MeDIP-seq of sample 2 in multinomial distribution.

P3

The probability of MRE-seq of sample 1 in multinomial distribution.

P4

The probability of MRE-seq of sample 2 in multinomial distribution.

C1

A constant to balance MeDIP-seq of sample 1 and sample 2.

C2

A constant to balance MRE-seq of sample 1 and sample 2.

pvalue

p-values of windows.

Value

p-value.

Examples

T<-4
SIZE<-200
p<-c(0.25,0.25,0.25,0.25)
c1<-1
c2<-1
length<-1
f<-pmultinom(T, SIZE,length, p[1], p[2], p[3], p[4], c1,
 c2, pvalue=rep(0,length(T)))