[Rsiena-help] Error in Print01Report for multigroup analysis

[marco tonellato]

Thank you for your help Tom.

I'm preparing the data object (it's not too big) and sending it to you.

All the best,
Marco

On Wed, Feb 25, 2015 at 4:56 PM, Tom Snijders <
Tom.Snijders at nuffield.ox.ac.uk> wrote:

>  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/b8dc08e2/attachment-0001.html>