errors when using RSDA2SymbolicDA function from symbolicDA package

[simonpedrog]

I have issues with RSDA2SymbolicDA. It seems to work flawlessly with examples such as 'oils' data from RSDA package, but i could not make it work with my data. I have tried using a symbolic object created with RSDA classic.to.sym function, and also loading a sym.data.table from a .csv.
Im using:
rstudio v1.2.5001
R v3.6.1
symbolicDA v0.6.2
RSDA v2.0.3

library("RSDA")
library("symbolicDA")

data <- data.frame(
                   Edad = c(9L, 5L, 71L, 69L, 56L, 22L, 26L, 75L, 35L, 69L,
                            80L, 64L, 24L, 22L, 66L, 61L, 30L, 28L, 40L, 17L),
   RelaciondeParentesco = as.factor(c("Hijo/Hijastro/a", "Hijo/Hijastro/a",
                                      "Jefe/a", "Cónyuge/Pareja", "Jefe/a",
                                      "Hijo/Hijastro/a", "Jefe/a", "Jefe/a",
                                      "Jefe/a", "Jefe/a", "Cónyuge/Pareja", "Jefe/a",
                                      "Hijo/Hijastro/a", "Hijo/Hijastro/a",
                                      "Jefe/a", "Cónyuge/Pareja", "Hijo/Hijastro/a",
                                      "Hijo/Hijastro/a", "Jefe/a",
                                      "Hijo/Hijastro/a")),
                   Sexo = as.factor(c("Mujer", "Mujer", "Varon", "Mujer",
                                      "Mujer", "Mujer", "Mujer", "Mujer",
                                      "Varon", "Mujer", "Varon", "Varon", "Varon",
                                      "Varon", "Mujer", "Varon", "Mujer", "Varon",
                                      "Mujer", "Varon")),
            EstadoCivil = as.factor(c("Soltero", "Soltero", "Casado", "Casado",
                                      "Separado/divorciado", "Soltero",
                                      "Soltero", "Soltero", "Soltero", "Casado",
                                      "Casado", "Separado/divorciado", "Soltero",
                                      "Soltero", "Casado", "Casado", "Soltero",
                                      "Soltero", "Separado/divorciado", "Soltero")),
         NivelEducativo = as.factor(c("SinInstruccion", "SinInstruccion",
                                      "UnivCompleta", "UnivCompleta",
                                      "UnivCompleta", "UnivIncompleta", "UnivIncompleta",
                                      "SinInstruccion", "UnivCompleta",
                                      "UnivCompleta", "UnivCompleta", "UnivIncompleta",
                                      "UnivCompleta", "UnivIncompleta",
                                      "UnivCompleta", "Secundaria", "UnivIncompleta",
                                      "UnivIncompleta", "Secundaria", "Secundaria")),
           EstadoSocial = as.factor(c("MenorDe10Años", "MenorDe10Años",
                                      "Ocupado", "Ocupado", "Inactivo",
                                      "Ocupado", "Ocupado", "Inactivo", "Ocupado",
                                      "Inactivo", "Inactivo", "Inactivo",
                                      "Desocupado", "Inactivo", "Ocupado", "Ocupado",
                                      "Inactivo", "Ocupado", "Ocupado", "Ocupado")),
       ActividadLaboral = as.factor(c("NoCorresponde", "NoCorresponde",
                                      "Obrero/Empleado", "CuentaPropia",
                                      "NoCorresponde", "Obrero/Empleado",
                                      "Obrero/Empleado", "NoCorresponde", "Obrero/Empleado",
                                      "NoCorresponde", "NoCorresponde",
                                      "NoCorresponde", "NoCorresponde", "NoCorresponde",
                                      "Obrero/Empleado", "CuentaPropia",
                                      "NoCorresponde", "Obrero/Empleado",
                                      "Obrero/Empleado", "Obrero/Empleado")),
     RazondeInactividad = as.factor(c("Estudiante", "MenorDe6Años",
                                      "NoCorresponde", "NoCorresponde",
                                      "Jubilado/Pensionado", "NoCorresponde",
                                      "NoCorresponde", "Jubilado/Pensionado",
                                      "NoCorresponde", "Jubilado/Pensionado",
                                      "Jubilado/Pensionado", "Jubilado/Pensionado",
                                      "NoCorresponde", "Estudiante", "NoCorresponde",
                                      "NoCorresponde", "Estudiante", "NoCorresponde",
                                      "NoCorresponde", "NoCorresponde")),
           BuscaTrabajo = as.factor(c("NoCorresponde", "NoCorresponde",
                                      "NoCorresponde", "NoCorresponde", "No",
                                      "NoCorresponde", "NoCorresponde", "No",
                                      "NoCorresponde", "No", "No", "No",
                                      "NoCorresponde", "No", "NoCorresponde",
                                      "NoCorresponde", "No", "NoCorresponde", "NoCorresponde",
                                      "NoCorresponde")),
    Ing_trabajo_estatal = as.factor(c("No", "No", "Si", "No", "No", "No", "No",
                                      "No", "No", "No", "No", "No", "No", "No",
                                      "No", "No", "No", "No", "No", "No")),
    Ing_trabajo_privado = as.factor(c("No", "No", "No", "No", "No", "Si", "Si",
                                      "No", "No", "No", "No", "No", "No", "No",
                                      "No", "Si", "No", "Si", "Si", "Si")),
         Ing_jubilación = as.factor(c("No", "No", "No", "No", "Si", "No", "No",
                                      "Si", "No", "Si", "Si", "No", "No", "No",
                                      "No", "No", "No", "No", "No", "No")),
           Ing_subsidio = as.factor(c("No", "No", "No", "No", "No", "No", "No",
                                      "No", "No", "No", "No", "No", "No", "No",
                                      "No", "No", "No", "No", "No", "No"))
)

variables = c('RelaciondeParentesco'
              ,'Sexo'
              , 'Edad'
              , 'EstadoCivil'
              , 'NivelEducativo'
              ,'EstadoSocial'
              ,'ActividadLaboral', 'RazondeInactividad', 'BuscaTrabajo'
              ,'Ing_trabajo_estatal', 'Ing_trabajo_privado', 'Ing_jubilación', 'Ing_subsidio'
             )
variables.types = c('RelaciondeParentesco' = type.modal(),
                    'Sexo' = type.modal(),
                    'Edad' = type.interval(),
                    'EstadoCivil' = type.modal(),
                    'NivelEducativo' = type.modal(),
                    'ActividadLaboral' = type.modal(),
                    'RazondeInactividad' = type.modal(),
                    'BuscaTrabajo' = type.modal(),
                    'Ing_trabajo_estatal' = type.modal(),
                    'Ing_trabajo_privado' = type.modal(),
                    'Ing_jubilación' = type.modal(),
                    'Ing_subsidio' = type.modal()
                    )
concept = 'EstadoSocial'
sym.table <- classic.to.sym(data = data,
                            concept = concept,
                            variables = variables,
                            variables.types = variables.types
                            )

#The problem is in the next line
symDA.table <- RSDA2SymbolicDA(sym.table)

#>[1] "tt"
#>Error in detailsNMA[detailsNMNo, 3] <- nr : replacement has length zero

[andresrcs]

@simonpedrog

I have tried uploading the data so you can reproduce, but i get 'format not authorized'.

There is no need to upload the file, follow this guide to make a proper reproducible example

FAQ: How to do a minimal reproducible example ( reprex ) for beginners

A minimal reproducible example consists of the following items: A minimal dataset, necessary to reproduce the issue The minimal runnable code necessary to reproduce the issue, which can be run on the given dataset, and including the necessary information on the used packages. Let's quickly go over each one of these with examples: Minimal Dataset (Sample Data) You need to provide a data frame that is small enough to be (reasonably) pasted on a post, but big enough to reproduce your issue. Let's say, as an example, that you are working with the iris data frame head(iris) #> Sepal.Length Sepal.Width Petal.Length Petal.Width Species #> 1 5.1 3.5 1.4 0.…

[simonpedrog]

Thanks, i have changed it. I hope it is ok now.