examples/example4.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: Repeated measure
## Data: [WA income data](examples/data/2019_WA_income.csv)

library(rTisane)

pid <- Participant("ID")
age <- numeric(unit=pid, name="AGEP")
edu <- nominal(unit=pid, name="SCHL", cardinality=10)
sex <- nominal(unit=pid, name="SEX", cardinality=2)
year <- Time(name="year", cardinality=5) 
income <- numeric(unit=pid, name="PINCP", numberOfInstances=per(1, year)) # One for each year

# Construct Conceptual Model
cm <- ConceptualModel()
# Specify conceptual relationships
cm <- assume(causes(age, income), cm)
cm <- hypothesize(relates(edu, income), cm)
cm <- hypothesize(whenThen(when=equals(sex, 1.0), then=increases(income)), cm)
conditions <- list(equals(sex, 1.0), increases(edu))
cm <- hypothesize(whenThen(when=conditions, then=increases(income)), cm) # Interaction

# 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, data="examples/data/2019_WA_income.csv")