cvborkulo / NetworkComparisonTest

Statistical comparison of two networks with respect to three invariance measures
21 stars 21 forks source link

problem occurs when set estimator = 'mgm' #51

Open Daiyan-project opened 4 weeks ago

Daiyan-project commented 4 weeks ago

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.