#R CODE FOR BENJAMINI-HOCHBERG PROCEDURE: FAMILYWISE BY PREDICTOR

#code for creating vectors of p values for families of analyses; families grouped by predictor & outcome
pssrqa<-c(.208,.209,.217) #p values for general stress predicting RQ; actor effects
pssnria<-c(.175,.195,.209) #p values for general stress predicting NRI; actor effects
pssrqp<-c(.290,.293,.295) #p values for general stress predicting RQ; parnter effects
pssnrip<-c(.800,.778,.807) #p values for general stress predicting NRI; partner effects

ismrqa<-c(.009,.009,.010) #p values for internalized stigma predicting RQ; actor effects
ismnria<-c(.002,.002,.002) #p values for internalized stigma predicting NRI; actor effects
ismrqp<-c(.504,.512,.512) #p values for internalized stigma predicting RQ; partner effects
ismnrip<-c(.071, .058,.054) #p values for internalized stigma predicting NRI; partner effects

micrqa<-c(.180,.175,.178) #p values for microaggressions predicting RQ; actor effects
micnria<-c(.00003,.00005,.00007) #p values for microaggressions predicting NRI; actor effects
micrqp<-c(.863,.892,.893) #p values for microaggressions predicting RQ; partner effects
micnrip<-c(.289,.278,.292) #p values for microaggressions predicting NRI; partner effects

vicrqa<-c(.108,.105,.106) #p values for victimization predicting RQ; actor effects
vicnria<-c(.886,.890,.86) #p values for victimization predicting NRI; actor effects
vicrqp <-c(.483,.503,.506) #p values for victimization predicting RQ; partner effects
vicnrip<-c(.747,.755,.880) #p values for victimization predicting NRI; partner effects

sodrqa<-c(.015,.013,.013) #p values for outness predicting RQ; actor effects
sodnria<-c(.141,.134,.145) #p values for outness predicting NRI; actor effects
sodrqp<-c(.999,.976,.982) #p values for outness predicting RQ; partner effects
sodnrip<-c(.727,.706,.666) #p values for outness predicting NRI; partner effects

rqdepa<-c(.555,.266,.448,.380,.302) #p values for RQ predicting depression; actor effects
rqaa<-c(.212,.785,.234,.231,.924) #p values for RQ predicting AUDIT; actor effects
rqca<-c(.861,.328,.928,.945,.436) #p values for RQ predicting CUDIT; actor effects
rqdepp<-c(.897,.894,.717,.636,.944) #p values for RQ predicting depression; partner effects
rqap<-c(.796,.866,.841,.892,.963) #p values for RQ predicting AUDIT; partner effects
rqcp<-c(.823,.650,.878,.731,.770) #p values for RQ predicting CUDIT; partner effects

nridepa <-c(.001,.006,.001,.0002,.023) #p values for NRI predicting depression; actor effects
nriaa<-c(.019,.018,.030,.018,.052) #p values for NRI predicting AUDIT; actor effects
nrica<-c(.620,.889,.539,.595,.660) #p values for NRI predicting CUDIT; actor effects
nridepp<-c(.360,.753,.365,.310,.633) #p values for NRI predicting depression; partner effects
nriap<-c(.861,.954,.895,.810,.993) #p values for NRI predicting AUDIT; partner effects
nricp<-c(.619,.355,.752,.848,.546) #p values for NRI predicting CUDIT; partner effects

msdepa<-c(.007,.867,.544,.694,.621) #p values for stressors predicting depression; actor effects
msaa<-c(.777,.582,.492,.861,.710) #p values for stressors predicting AUDIT; actor effects
msca<-c(.199,.375,.933,.047,.692) #p values for stressors predicting CUDIT; actor effects
msdepp<-c(.865,.097,.783,.394,.220) #p values for stressors predicting depression; partner effects
msap<-c(.362,.676,.745,.679,.491) #p values for stressors predicting AUDIT; partner effects
mscp<-c(.240,.150,.870,.245,.806) #p values for stressors predicting CUDIT; partner effects

#Benjamini-Hochberg corrections to p values
p.adjust(pssrqa,"BH",n=length(pssrqa))
p.adjust(pssnria,"BH",n=length(pssnria))
p.adjust(pssrqp,"BH",n=length(pssrqp))
p.adjust(pssnrip,"BH",n=length(pssnrip))


p.adjust(ismrqa,"BH",n=length(ismrqa))
p.adjust(ismnria,"BH",n=length(ismnria))
p.adjust(ismrqp,"BH",n=length(ismrqp))
p.adjust(ismnrip,"BH",n=length(ismnrip))

p.adjust(micrqa,"BH",n=length(micrqa))
p.adjust(micnria,"BH",n=length(micnria))
p.adjust(micrqp,"BH",n=length(micrqp))
p.adjust(micnrip,"BH",n=length(micnrip))

p.adjust(vicrqa,"BH",n=length(vicrqa))
p.adjust(vicnria,"BH",n=length(vicnria))
p.adjust(vicrqp,"BH",n=length(vicrqp))
p.adjust(vicnrip,"BH",n=length(vicnrip))

p.adjust(sodrqa,"BH",n=length(sodrqa))
p.adjust(sodnria,"BH",n=length(sodnria))
p.adjust(sodrqp,"BH",n=length(sodrqp))
p.adjust(sodnrip,"BH",n=length(sodnrip))

p.adjust(rqdepa,"BH",n=length(rqdepa))
p.adjust(rqaa,"BH",n=length(rqaa))
p.adjust(rqca,"BH",n=length(rqca))
p.adjust(rqdepp,"BH",n=length(rqdepp))
p.adjust(rqap,"BH",n=length(rqap))
p.adjust(rqcp,"BH",n=length(rqcp))

p.adjust(nridepa,"BH",n=length(nridepa))
p.adjust(nriaa,"BH",n=length(nriaa))
p.adjust(nrica,"BH",n=length(nrica))
p.adjust(nridepp,"BH",n=length(nridepp))
p.adjust(nriap,"BH",n=length(nriap))
p.adjust(nricp,"BH",n=length(nricp))

p.adjust(msdepa,"BH",n=length(msdepa))
p.adjust(msaa,"BH",n=length(msaa))
p.adjust(msca,"BH",n=length(msca))
p.adjust(msdepp,"BH",n=length(msdepp))
p.adjust(msap,"BH",n=length(msap))
p.adjust(mscp,"BH",n=length(mscp))