Journal Stats

#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
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
journdat$is.Evolutionary.Biology[grep("*Evolutionary Biology*", journdat$ISI.Category)] = 1
journdat$is.Analytical.Chemistry = 0
journdat$is.Analytical.Chemistry[grep("*Analytic Chemistry*", journdat$ISI.Category)] = 1
journdat$is.Biology = 0
journdat$is.Biology[grep("*Biology*", journdat$ISI.Category)] = 1
journdat$is.Geosciences.Multidisciplinary = 0
journdat$is.Geosciences.Multidisciplinary[grep("*Geosciences Multidisciplinary*", journdat$ISI.Category)] = 1
journdat$is.Biodiversity.Conservation = 0
journdat$is.Biodiversity.Conservation[grep("*Biodiversity Conservation*", journdat$ISI.Category)] = 1

##Publishers
journdat$is.Elsevier = 0
journdat$is.Elsevier[grep("*Elsevier*", journdat$Publisher)] = 1
journdat$is.Springer = 0
journdat$is.Springer[grep("*Springer*", journdat$Publisher)] = 1
journdat$is.Wiley = 0
journdat$is.Wiley[grep("*Wiley*", journdat$Publisher)] = 1
journdat$is.Wiley.Blackwell = 0
journdat$is.Wiley.Blackwell[grep("*Wiley.Blackwell*", journdat$Publisher)] = 1
journdat$is.Taylor.and.Francis = 0
journdat$is.Taylor.and.Francis[grep("*Taylor.and.Francis*", journdat$Publisher)] = 1

##Summary/Response/Plot
summary(simple.var ~ Impact.Factor + Articles + Data.required.for.peer.review + Annual.publications + Policy.has.instructions.for.how.to.cite.data. + Type.of.place.to.deposit +  Subscription.Model + is.Environmental.Sciences + is.Ecology + is.Evolutionary.Biology + is.Analytical.Chemistry  + is.Biology  + is.Geosciences.Multidisciplinary + is.Biodiversity.Conservation, journdat)

response = summary(simple.var ~ Impact.Factor + Articles + Data.required.for.peer.review + Annual.publications + Policy.has.instructions.for.how.to.cite.data. + Type.of.place.to.deposit +  Subscription.Model + is.Environmental.Sciences + is.Ecology + is.Evolutionary.Biology + is.Analytical.Chemistry  + is.Biology  + is.Geosciences.Multidisciplinary + is.Biodiversity.Conservation, journdat)

plot(response)

##ConfInt

binconf(table(journdat$simple.var)[2], sum(table(journdat$simple.var))) ## Confidence Interval for req code