2nd example

Author: n8thangreen

inst/joss_code.Rmd

@@ -46,7 +46,7 @@ outstandR_gcomp_ml <-
 
 outstandR_gcomp_stan <-
   outstandR(AC_IPD_binY_contX, BC_ALD_binY_contX,
-            refresh = 0,  # turns off the iteration updates.
+            refresh = 0,  # turns off the iteration updates
             strategy = strategy_gcomp_stan(
               formula = lin_form,
               family = binomial(link = "logit")))
@@ -58,3 +58,52 @@ outstandR_mim <-
               formula = lin_form,
               family = binomial(link = "logit")))
 ```
+
+```{r}
+data("AC_IPD_contY_mixedX")
+data("BC_ALD_contY_mixedX")
+
+# individual-level data
+head(AC_IPD_contY_mixedX)
+
+# aggregate-level data
+BC_ALD_contY_mixedX
+
+lin_form <- as.formula("y ~ X1 + X2 + X3 + trt + trt:(X1 + X2 + X4)")
+
+outstandR_maic <-
+  outstandR(AC_IPD_contY_mixedX, BC_ALD_contY_mixedX,
+            strategy = strategy_maic(
+              formula = lin_form,
+              family = continuous(link = "identity")))
+
+str(outstandR_maic)
+
+print(outstandR_maic)
+
+outstandR_stc <-
+  outstandR(AC_IPD_contY_mixedX, BC_ALD_contY_mixedX,
+            strategy = strategy_stc(
+              formula = lin_form,
+              family = continuous(link = "identity")))
+
+outstandR_gcomp_ml <-
+  outstandR(AC_IPD_contY_mixedX, BC_ALD_contY_mixedX,
+            strategy = strategy_gcomp_ml(
+              formula = lin_form,
+              family = continuous(link = "identity")))
+
+outstandR_gcomp_stan <-
+  outstandR(AC_IPD_contY_mixedX, BC_ALD_contY_mixedX,
+            refresh = 0,  # turns off the iteration updates
+            strategy = strategy_gcomp_stan(
+              formula = lin_form,
+              family = continuous(link = "identity")))
+
+outstandR_mim <-
+  outstandR(AC_IPD_contY_mixedX, BC_ALD_contY_mixedX,
+            refresh = 0,
+            strategy = strategy_mim(
+              formula = lin_form,
+              family = continuous(link = "identity")))
+```