[Rsiena-help] Error in Print01Report for multigroup analysis

[marco tonellato]

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/8889b9ba/attachment.html>