required_plot

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
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


summary(requires ~ 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(requires ~ 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[journdat$simple.var > 0,])

plot(response)