examples/example6.R

In emjun/tisaner: Authoring Statistical Models from Conceptual Models

## Example 5 (from summative eval Income dataset)
## Conceptual relationships: Assume Causes, Assume Relates, Hypothesize whenThen,
## Data collection: Hierarchical nesting
## Data: [WA income data](examples/data/2019_WA_income.csv)

library(rTisane)

pid <- Participant("ID")

age <- continuous(unit=pid, name="AGEP")
edu <- categories(unit=pid, name="SCHL", cardinality=10)
sex <- categories(unit=pid, name="SEX", cardinality=2)
income <- continuous(unit=pid, name="PINCP")

# Construct Conceptual Model
cm <- ConceptualModel() %>%
# Specify conceptual relationships
#   assume(causes(age, income)) %>%
  hypothesize(relates(edu, income)) %>%
  hypothesize(relates(sex, income, when=equals(sex, 1.0), then=increases(income))) %>% 
  interacts(sex, edu, dv=income)

# Education is IV, Income is DV
# Need to disambiguate relationship between Education and Income, between Sex and Income in conceptual model
query(conceptualModel=cm, iv=edu, dv=income) # No data