network: This function enumerates conditional sample spaces associated...

Description Usage Arguments Details Value References Examples

View source: R/network.R

Description

This function uses a network algorithm to enumerate conditional sample spaces associated with logistic regression, using a minimal version of the algorithm of \insertCitemehtapatel;textualPHInfiniteEstimates.

Usage

1
2
3
4
5
6
7
8
9
network(
  dm,
  n = NULL,
  resp = NULL,
  conditionon = NULL,
  sst = NULL,
  addint = TRUE,
  verbose = FALSE
)

Arguments

dm

matrix of covariates

n

Vector of number of trials. If null, make them all ones.

resp

vector of successes. Used only to calculate the sufficient statistics, unless sufficient statistics are entered directly. Either resp or sst must be provided.

conditionon

indices of covariate matrix indicating sufficient statistics to be conditioned on.

sst

sufficient statistic vector, if input directly. Otherwise, recomputed from resp.

addint

logical, true if a column of 1s must be added to the covariate matrix.

verbose

logical; if true, print intermediate results.

Details

Examples apply this to data from \insertCitemehtapatel;textualPHInfiniteEstimates, citing \insertCitegoorinetal87;textualPHInfiniteEstimates.

Value

For a successful run, a list with components:

For an unsuccessful run (because of input inconsistencies) NA

References

\insertRef

mehtapatelPHInfiniteEstimates

\insertRef

goorinetal87PHInfiniteEstimates

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
#Columns in table are:
# Lymphocytic Infiltration (1=low, 0=high)
# Sex (1=male, 0=female)
# Any Ostioid Pathology (1=yes, 0=no)
# Number in LI-Sex-AOP group
# Number in LI-Sex-AOP group with disease free interval greater than 3 y
goorin<-data.frame(LI=c(0,0,0,0,1,1,1,1),Sex=c(0,0,1,1,0,0,1,1),
   AOP=c(0,1,0,1,0,1,0,1),N=c(3,2,4,1,5,5,9,17),Y=c(3,2,4,1,5,3,5,6))

out<-network(goorin[,1:3],goorin[,4],conditionon=1:3,resp=goorin[,5])
inference(out)

PHInfiniteEstimates documentation built on July 24, 2021, 1:07 a.m.