-
Notifications
You must be signed in to change notification settings - Fork 21
Description
Hi,
I want to compare network with the setting of mgm. I don't know if it is right because I didn't find the option about parameter 'estimator'.
This is the code of a easy example:
library(NetworkComparisonTest)
library(mgm)
set.seed(123)
data1_test <- matrix(rnorm(1005), ncol=5)
data2_test <- matrix(rnorm(1005), ncol=5)
data1_test <- as.matrix(data1_test)
data2_test <- as.matrix(data2_test)
nct_result_mgm <- NetworkComparisonTest::NCT(data1_test, data2_test,
it = 1000,
estimator = "mgm",
estimatorArgs = list(type = rep("g", ncol(data1_test)),
level = rep(1, ncol(data1_test)))
)
This is the error:
|----------------------------------------------------------------------------------------| 100%
|----------------------------------------------------------------------------------------| 100%
Note that the sign of parameter estimates is stored separately; see ?mgmError in array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), :
'data' must be of a vector type, was 'NULL'
Anyway, thanks for your kindly help.