[Rsiena-help] Possible regression with sienaTimeTest on 2-mode networks

Guido Conaldi conaldi at gmail.com
Thu Jan 13 00:57:33 CET 2011


Dear All,

I was checking whether the latest versions of RSiena/Test had solved the
problem we encountered when using includeTimeDummy() on sienaEffects
objects with 2-mode networks. I found instead a possible regression in
sienaTimeTest().


-1-
This is what I get with the latest RSiena/Test and a sienaFit object
with 3 time periods and one two-mode network (and I tried many different
ones):


> sienaTimeTest(a.compAll.mainEff)
Error in sienaTimeTest(a.compAll.mainEff) : 
You must have at least three time periods to test for non-heterogeneity
across time.


> a.compAll.mainEff
Estimates, standard errors and t-statistics for convergence

                                                  Estimate   Standard   t statistic 
                                                              Error                

Rate parameters: 
  0.1      Rate parameter period 1                 3.5386  ( 0.2790   )             
  0.2      Rate parameter period 2                 5.1183  ( 0.4053   )

...


With RSiena 1.0.11.120 RSienaTest 1.0.12.123 I get the expected output.


-2-
This is the error I get with all versions of RSiena when trying to
include a timedummy in a sienaEffect object with a two-mode network:


>
epi.effects[[7]]<-includeTimeDummy(epi.effects[[7]],cycle4,timeDummy="2")

** caught segfault ***
address 0x10, cause 'memory not mapped'

Traceback:
1: .Call("ChangingCovariates", PACKAGE = pkgname, pData, lapply(f,
unction(x) x$vCovars))
2: initializeFRAN(z, x, data, effects, prevAns, initC, profileData =
profileData, returnDeps = returnDeps)
3: x$FRAN(z, x, INIT = TRUE, ...)
4: robmon(z, x, useCluster, nbrNodes, initC, clusterString,
clusterIter, ...)
5: siena07(epi.model[[i]], data = epi.data[[i]], effects =
epi.effects[[i]], returnDeps = FALSE, batch = TRUE, verbose = TRUE,
nbrNodes = 3, useCluster = TRUE, initC = TRUE)

I can provide my objects in case you need them.

I hope you can help me,
Thanks in advance and a Happy new year to all of you!

Guido



More information about the Rsiena-help mailing list