[Rsiena-help] Error in Print01Report for multigroup analysis

[Tom Snijders]

Dear Marco,

If you send me the data object used (or is it too large to send?) I can try to have a look.

Cheers,
Tom


================================================================
Tom A.B. Snijders
Professor of Statistics and Methodology
Dept. of Sociology, University of Groningen
Emeritus Fellow, Nuffield College, University of Oxford
Associate Member, Dept. of Statistics, University of Oxford
http://www.stats.ox.ac.uk/~snijders/

From: rsiena-help-bounces at lists.r-forge.r-project.org [mailto:rsiena-help-bounces at lists.r-forge.r-project.org] On Behalf Of marco tonellato
Sent: 25 February 2015 14:35
To: rsiena-help at lists.r-forge.r-project.org
Subject: [Rsiena-help] Error in Print01Report for multigroup analysis

Dear All,
while importing data in Rsiena (revision 278) for a test multigroup meta-analysis (7 groups, network size between 14 and 39, 1 dependent network and 1 constant dyadic covariate) I encountered the following error during the Print01Report routine.

Error in format(round(atts$meanp[j], 3), nsmall = 3, width = 10) :
  error in evaluating the argument 'x' in selecting a method for function 'format': Error in round(atts$meanp[j], 3) :
  non-numeric argument to mathematical function
Here is the code I used:

groups <- list()
groups[[1]] <- list(advice=array(c(g1t1_advice_adjacency, g1t2_advice_adjacency), dim=c(dim(g1t1_advice_adjacency),2)), energy=g1t1_energy_adjacency)
groups[[2]] <- list(advice=array(c(g2t1_advice_adjacency, g2t2_advice_adjacency), dim=c(dim(g2t1_advice_adjacency),2)), energy=g2t1_energy_adjacency)
groups[[3]] <- list(advice=array(c(g3t1_advice_adjacency, g3t2_advice_adjacency), dim=c(dim(g3t1_advice_adjacency),2)), energy=g3t1_energy_adjacency)
groups[[4]] <- list(advice=array(c(g4t1_advice_adjacency, g4t2_advice_adjacency), dim=c(dim(g4t1_advice_adjacency),2)), energy=g4t1_energy_adjacency)
groups[[5]] <- list(advice=array(c(g5t1_advice_adjacency, g5t2_advice_adjacency), dim=c(dim(g5t1_advice_adjacency),2)), energy=g5t1_energy_adjacency)
groups[[6]] <- list(advice=array(c(g6t1_advice_adjacency, g6t2_advice_adjacency), dim=c(dim(g6t1_advice_adjacency),2)), energy=g6t1_energy_adjacency)
groups[[7]] <- list(advice=array(c(g7t1_advice_adjacency, g7t2_advice_adjacency), dim=c(dim(g7t1_advice_adjacency),2)), energy=g7t1_energy_adjacency)

advicelist <- list()
for (i in 1:7){advicelist[[i]] <- sienaDependent(groups[[i]]$advice)}

energylist <- list()
for (i in 1:7){energylist[[i]] <- coDyadCovar(groups[[i]]$energy)}

groups_data_list <- list()
for (i in 1:7){groups_data_list[[i]] <- sienaDataCreate(advice = advicelist[[i]], energy = energylist[[i]])}

groups_data <- sienaGroupCreate(groups_data_list)
summary(groups_data)
groups_effects <- getEffects(groups_data)
print01Report(groups_data, modelname = 'groups_descriptives')

Do you have any clues on how to fix this?
Thank you.
Marco
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/rsiena-help/attachments/20150225/78423022/attachment.html>