summary.lmFAB: Summarizing Linear Model Fits with FAB Inference

Description Usage Arguments Details Value Examples

Description

summary method for class lmFAB

Usage

1
2
3
## S3 method for class 'lmFAB'
summary(object, correlation = FALSE,
  symbolic.cor = FALSE, ...)

Arguments

object

an object of class lmFAB

correlation

see summary.lm

symbolic.cor

see summary.lm

...

see summary.lm

Details

A mod of summary.lm that shows FAB p-values in table

Value

A list of summary statistics of the fitted linear model

Examples

 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
# n observations, p FAB variables, q=2 control variables 

n<-100 ; p<-25 

# X is design matrix for params of interest
# beta is vector of true parameter values 
# v a variable in the linking model - used to share info across betas

v<-rnorm(p) ; beta<-(2 - 2*v + rnorm(p))/3 ; X<-matrix(rnorm(n*p),n,p)/8

# control coefficients and variables  
alpha1<-.5 ; alpha2<- -.5
w1<-rnorm(n)/8
w2<-rnorm(n)/8

# simulate data 
lp<-1 + alpha1*w1 + alpha2*w2 + X%*%beta 
y<-rnorm(n,lp) 

# fit model
fit<-lmFAB(y~w1+w2,X,~v)

fit$FABpv
fit$FABci 
summary(fit) # look at p-value column 

FABInference documentation built on Jan. 9, 2020, 5:08 p.m.