A function to fit a block recursive chain graph model in R following the outlines of Cox & Wermuth 1996 and Caputo et al 1997; some adaptions made by us

1 2 3 4 |

`fmla` |
A formula listing the block structure and the variables in each block. It needs be of the form Block1Vars ~ Block2Vars ~ Block3Vars and the variables per block are separated by +. So a formula y1+y2 ~ x1+x2 ~ z1~z2 refers to the model wehere y1 and y2 are in block 1 and are purely dependent variables, x1 and x2 are in block 2 and are intermediate (i.e. predictors for block 1 and dependents for block 3) and z1 and z3 are in block 3 and are purely explanatory. Internally, the var.frame is built from the formula and the vartypes is automatically detected (but that is crude). |

`data` |
the data frame containing the variables in var.frame or formula |

`var.frame` |
If no formula is given a variable frame with the following structure: First column are the variable names, second column are the variable types (one of categorical, ordinal, continuous, binary, count, odcount (overdispersed count), gamma or invgaussian, see also |

`vartype` |
(optional) the types of variables in data and formula; if formula is given but vartype is not, the function attempts to detect the type of variable and will select one of ordinal, categorical and continuous. The order of types in vartype is expected to be the order of variables in formula (from left to right) or in var.frame from top to bottom |

`automatch` |
flag whether the function should match the variable type to the data object (TRUE if so); will override user specifications in var.frame; defaults to FALSE |

`pen` |
the penalty applied for the numbers of parameters to the information criterion for variable selection; defaults to log(n) (BIC); 2 would be AIC |

`signif` |
the significance level used in the tests for inclusion of higher order effects; defaults to 0.01 |

`contrasts` |
which contrasts to use for ordinal and nominal predictors; defaults to treatment contrasts all around |

`scalepredictors` |
whtehr and how to scale the predictors (no=no scaling, center=mean centering, scale=standardization (mean=0 and var=1) and normalize=scaling to [0,1]) |

`silent` |
flags whether output should be printed (if TRUE, no output; the default) |

`restr` |
a matrix of the same size as the adjacency matrix in the model so as to supress specific dependencies fitted in the algorithm (0 means not fitted) |

`adjfile` |
if TRUE the adjacency matrix is written to a file |

`...` |
additional arguments (not used currently) |

an object of class cw which contains a list of fitted models, the variable frame and the fitted adjacency matrix

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | ```
data(cmc)
#this is a 5 block model with age being
#purely explanatory (block 5) and nrChild and
#contraceptive being purely the targets (block 1)
#Using a formula and autdetection for models
#vartype is automatically detected which is crude for non-factors
res_cmc2<-coxwer(contraceptive + nrChild ~ #block 1 - purely targets
mediaExp ~ #block 2
solIndex ~ #block 3
wifeRel + wifeWork + husbOcc + wifeEdu + husbEdu ~ #block 4
age, #block 5 - purely explanatory
data=cmc)
print(res_cmc2) #Prints adjacency matrix
summary(res_cmc2,target="nrChild") #model path to "contraceptive" as the target variable
predict(res_cmc2,target="nrChild") #predict method
## Not run:
#Using a formula and specifying the vartype
#so the algorithm can use better models for the metric/continuous variables
#(here Poisson model for nrChild instead of OLS)
res_cmc3<-coxwer(contraceptive + nrChild ~
mediaExp ~
solIndex ~
wifeRel + wifeWork + husbOcc + wifeEdu + husbEdu ~
age,
vartype=c("cate","count","bin","ord","bin","bin","ord","ord","ord","metric"),
data=cmc)
summary(res_cmc3,target="nrChild")
#Using a var.frame
#cmc_prep<-prep_coxwer(cmc)
data(cmc_prep)
cmc_prep
res_cmc <- coxwer(var.frame=cmc_prep, data=cmc)
## End(Not run)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.