85 lines
3.0 KiB
R
85 lines
3.0 KiB
R
|
|
## replace the below paths with your actual path to dropbox
|
|
## e.g., '/home/kaylea/Dropbox/Apps/Overleaf/WhateverICalledMyProject/figures/'
|
|
figPath = '../InOverleaf/figures/'
|
|
## e.g., '/home/kaylea/Dropbox/Apps/Overleaf/WhateverICalledMyProject/knitr_rdata/'
|
|
knitrPath = '../InOverleaf/knitr_rdata/'
|
|
source('lib-00-utils.R') # enables the remember command
|
|
set.seed(424242) # always!
|
|
|
|
library('texreg')
|
|
library('ggplot2')
|
|
|
|
magicNumberOne <- rnorm(1, 0, 1)
|
|
magicNumberOne
|
|
magicNumberTwo <- rnorm(1, 0, 1)
|
|
magicNumberTwo
|
|
|
|
Importance <- data.frame(value=rbeta(100,3,6)+0.1, source='Importance')
|
|
Quality <- data.frame(value=cumsum(rbeta(10,1,4)/10), source='Quality')
|
|
sampleDF <- rbind(Importance, Quality)
|
|
head(sampleDF)
|
|
head(Quality)
|
|
head(Importance)
|
|
|
|
|
|
g <- ggplot(sampleDF, aes(value, color=source)) +
|
|
geom_density(data=Importance, alpha=0.25, size=1) +
|
|
stat_ecdf(data=Quality, aes(x=value), geom='step') +
|
|
theme_classic() +
|
|
xlab("Time") +
|
|
ylab("Standard Deviation Units") +
|
|
theme(axis.ticks.x=element_blank(), axis.text.x=element_blank(),
|
|
axis.ticks.y=element_blank(), axis.text.y=element_blank(),
|
|
axis.title.x=element_text(size=8),
|
|
axis.title.y=element_text(size=8),
|
|
legend.position="bottom", legend.title=element_blank(),
|
|
legend.key.size= unit(.1, 'cm'), legend.key.height = unit(.1, 'cm'),
|
|
legend.margin=margin(0,0,0,0),
|
|
legend.key.width = unit(.1, 'cm'), legend.text = element_text(size=8)) +
|
|
annotate("text", size=3, label="Risk", x=.32, y =1.4) +
|
|
geom_segment(aes(x=.275, y=1.3, xend=.4, yend=1.3), color='black', size=0.13) +
|
|
geom_segment(aes(x=.275, y=1.33, xend=.275, yend=1.27), color='black', size=0.14) +
|
|
geom_segment(aes(x=.4, y=1.33, xend=.4, yend=1.27), color='black', size=0.14)
|
|
|
|
g
|
|
|
|
png(paste0(figPath, "simCurve.png"), units = "in", res = 300, width = 3, height = 3)
|
|
g
|
|
dev.off()
|
|
|
|
m1 <- lm(Importance ~ Quality, data=sampleDF)
|
|
|
|
summary(m1)
|
|
|
|
m2 <- lm(Importance ~ poly(Quality, degree=2, raw=2), data=sampleDF)
|
|
|
|
summary(m2)
|
|
|
|
con <- textConnection("model.texreg", "w") #remembered
|
|
sink(con, split=TRUE, type="output")
|
|
texreg(list(m1,m2), omit.coef = 'factor', stars=NULL, digits=2,
|
|
custom.model.names=c('M1 Description', 'M2 Description'),
|
|
#custom.coef.names=c('(Intercept)', 'Quality', 'Quality Squared'), # run without specifying first, just to be sure these are right!
|
|
use.packages=FALSE, table=FALSE, ci.force = TRUE)
|
|
sink()
|
|
close(con);rm(con)
|
|
|
|
m1.coefs <- summary(m1)$coefficients ## for easy reference within the latex document
|
|
m2.coefs <- summary(m1)$coefficients ## for easy reference within the latex document
|
|
|
|
|
|
|
|
|
|
if (!nosave) {
|
|
r <- list()
|
|
remember(magicNumberOne) #put everything you want to have access to inside latex here
|
|
remember(magicNumberTwo) #put everything you want to have access to inside latex here
|
|
remember(model.texreg)
|
|
remember(m1.coefs)
|
|
remember(m2.coefs)
|
|
save(r, file=paste0(knitrPath, "main.RData"), version=2) ## this creates the file, and we need to specify the version
|
|
print("Remembrances complete")
|
|
}
|
|
|