nniiicc

Publisher	Journal Title	Journal Alias	ISSN	Publications / Year	ISI Category	Total Citaitons 2009	Impact Factor	5 Year Impact Factor	Immediacy Index	Articles	Cited Half Life	EigenFactor Score	Article Influence Score	Journal has affiliation	Affiliation	Peer Reviewed	Subscription Model	Sub Code	Location (url)	Standard publication policy	Date of plan	Data required for peer review	Required peer review text	Policy requested authors to share data  for publication	Policy required authors to share data for publication	Request / Require text	Policy request / require code	Policy has instructions how to submit data?	Policy has instructions for how to submit data? Text	Policy has instructions where to deposit shared data ?	Policy has instructions where to deposit shared data ? Text	Type of place to deposit? (Journal, Repository etc)	Policy requires accession number?	Policy requires format for shared data ?	Required format text	Policy has instructions for how to cite data?	How to cite data text	Policy has instructions for 

nniiicc

#JournalStats

journdat = read.csv("filename", as.is=TRUE, na.strings="N/A")
journdat$simple.var = journdat$Policy.request...require.code ## Renames field to simple.var
journdat$simple.var[journdat$simple.var>1] = 1 ## Turns all fields >1 into 1 for sake of simple variable of request / required code

##Categories

journdat$is.Environmental.Sciences = 0
journdat$is.Environmental.Sciences[grep("*Environmental Sciences*", journdat$ISI.Category)] = 1 ## way of making EnvSci as a category into a varaible in order to add it to our plot

nniiicc

#JournalStats

journdat = read.csv("filename", as.is=TRUE, na.strings="N/A")
journdat$simple.var = journdat$Policy.request...require.code ## Renames field to simple.var
journdat$simple.var[journdat$simple.var>1] = 1 ## Turns all fields >1 into 1 for sake of simple variable of request / required code

##Categories

journdat$is.Environmental.Sciences = 0
journdat$is.Environmental.Sciences[grep("*Environmental Sciences*", journdat$ISI.Category)] = 1 ## way of making EnvSci as a category into a varaible in order to add it to our plot

nniiicc

#JournalStats

journdat = read.csv("filename", as.is=TRUE, na.strings="N/A")
journdat$simple.var = journdat$Policy.request...require.code ## Renames field to simple.var
journdat$simple.var[journdat$simple.var>1] = 1 ## Turns all fields >1 into 1 for sake of simple variable of request / required code

##Categories

journdat$is.Environmental.Sciences = 0
journdat$is.Environmental.Sciences[grep("*Environmental Sciences*", journdat$ISI.Category)] = 1 ## way of making EnvSci as a category into a varaible in order to add it to our plot

nniiicc

journdat = read.csv("filename", as.is=TRUE, na.strings="N/A")

journdat$requires = 0
journdat$requires[journdat$Policy.request...require.code > 1] = 1

journdat$is.Environmental.Sciences = 0
journdat$is.Environmental.Sciences[grep("*Environmental Sciences*", journdat$ISI.Category)] = 1 ## way of making EnvSci as a category into a varaible in order to add it to our plot
journdat$is.Ecology = 0
journdat$is.Ecology[grep("*Ecology*", journdat$ISI.Category)] = 1 ## ditto from above, remember journals have more than one category classification
journdat$is.Evolutionary.Biology = 0

nniiicc

dotchart2.CIs = 
function (data, labels, groups = NULL, gdata = NA, horizontal = TRUE, 
    pch = 16, xlab = "", ylab = "", auxdata, auxgdata = NULL, 
    auxtitle, lty = if (.R.) 1 else 2, lines = TRUE, dotsize = 0.8, 
    cex = par("cex"), cex.labels = cex, cex.group.labels = cex.labels * 
        1.25, sort. = TRUE, add = FALSE, dotfont = par("font"), 
    groupfont = 2, reset.par = add, xaxis = TRUE, width.factor = 1.1, 
    lcolor = if (.R.) "gray" else par("col"), ...) 
{
    if (.R. && !add) {

nniiicc

journdat = read.csv("/Users/nicholasweber/Desktop/Finished DataSets/JournalData.csv", as.is=TRUE, na.strings="N/A")

journdat$affil.var = journdat$Affiliation.Code
journdat$affil.var[journdat$affil.var>1] = 1 

journdat$is.Environmental.Sciences = 0
journdat$is.Environmental.Sciences[grep("*Environmental Sciences*", journdat$ISI.Category)] = 1 ## way of making EnvSci as a category into a varaible in order to add it to our plot
journdat$is.Ecology = 0
journdat$is.Ecology[grep("*Ecology*", journdat$ISI.Category)] = 1 ## ditto from above, remember journals have more than one category classification
journdat$is.Evolutionary.Biology = 0

nniiicc

journdat = read.csv("/Users/JournalData.csv", as.is=TRUE, na.strings="N/A")

journdat$required.sharing = 0
journdat$required.sharing[journdat$Policy.request...require.code > 1] = 1


journdat$is.Environmental.Sciences = 0
journdat$is.Environmental.Sciences[grep("*Environmental Sciences*", journdat$ISI.Category)] = 1 
journdat$is.Ecology = 0
journdat$is.Ecology[grep("*Ecology*", journdat$ISI.Category)] = 1 ## ditto from above, remember journals have more than one category classification

nniiicc

mydata <- read.csv("~/Desktop/JournalData.csv")
attach(mydata)
names(mydata)

how.cite.data.var = Cite.data.code
how.cite.data.var [how.cite.data.var>1] = 1

mylogit<- glm(how.cite.data.var~Impact.Factor+Policy.request...require.code, family=binomial(link="logit"), na.action=na.pass)

summary(mylogit)

nniiicc

filename = "/Users/nicholasweber/Desktop/JournalDat.csv"
mydata = read.csv(filename)
attach(mydata)
requests = ifelse(Policy.request...require.code > 0, 1, 0)
ImFa = Impact.Factor
ImFa[ImFa==0] = NA
hist(ImFa)
summary(ImFa)
SomeOA = ifelse(Subscription.Model == "Sub", 0, 1)
table(SomeOA)