[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