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
|
rcModelPLMd <- function(y,group.labels){
if (!is.matrix(y))
stop("argument should be matrix")
if (length(group.labels) != ncol(y)){
stop("group labels is of incorrect length")
}
if (!is.factor(group.labels)){
group.labels <- as.factor(group.labels)
}
if (any(table(group.labels) < 2)){
stop("Must be at least two arrays in each group")
}
group.int <- as.integer(group.labels) -1
PsiCode <- 0
PsiK <- 1.345
.Call("R_plmd_model",y,PsiCode,PsiK,as.integer(group.int),as.integer(length(unique(group.labels))),PACKAGE="preprocessCore")
}
|
