Skip to content

Instantly share code, notes, and snippets.

@Oreotrephes
Last active January 8, 2016 15:33
Show Gist options
  • Select an option

  • Save Oreotrephes/3e18ea830c0620da2d49 to your computer and use it in GitHub Desktop.

Select an option

Save Oreotrephes/3e18ea830c0620da2d49 to your computer and use it in GitHub Desktop.
market.compare<-function(market.name,df,uses){
this.market.uses<-df[df$Market==market.name,uses]
other.market.uses<-df[df$Market!=market.name,uses]
this.market.unique<-setdiff(this.market.uses,other.market.uses)
this.market.total<-unique(this.market.uses)
percent.unique<-length(this.market.unique)/length(this.market.total)
percent.unique
}
market.compare("Calle Santa Cruz",sample,"Medicinal.Use.categories")
#BUT NOT
sample %>%
group_by(Scientific.name) %>%
summarize(
CSC=length(market.compare("Calle Santa Cruz",sample,"Medicinal.Use.categories")))
##OH THE KLUDGE
species<-c()
n.markets<-c()
Calle.Santa.Cruz<-c()
Cohoni<-c()
Cota.Cota<-c()
El.Alto<-c()
Final.Bajo.Seguencoma<-c()
Mercado.Rodriguez<-c()
for(i in levels(as.factor(RBU$Scientific.name))){
subs<-subset(RBU,Scientific.name==i)
n<-length(unique(subs$Market))
CSC<-market.compare("Calle Santa Cruz",subs,"Medicinal.Use.categories")
CH<-market.compare("Cohoni",subs,"Medicinal.Use.categories")
CC<-market.compare("Cota Cota",subs,"Medicinal.Use.categories")
EA<-market.compare("El Alto" ,subs,"Medicinal.Use.categories")
FBS<-market.compare("Final Bajo Seguencoma",subs,"Medicinal.Use.categories")
MR<-market.compare("Mercado Rodriguez",subs,"Medicinal.Use.categories")
species<-c(species,i)
n.markets<-c(n.markets,n)
Calle.Santa.Cruz<-c(Calle.Santa.Cruz,CSC)
Cohoni<-c(Cohoni,CH)
Cota.Cota<-c(Cota.Cota,CC)
El.Alto<-c(El.Alto,EA)
Final.Bajo.Seguencoma<-c(Final.Bajo.Seguencoma,FBS)
Mercado.Rodriguez<-c(Mercado.Rodriguez,MR)
}
data.frame(cbind(species,n.markets,Calle.Santa.Cruz,Cohoni,Cota.Cota,El.Alto,Final.Bajo.Seguencoma,Mercado.Rodriguez))->output
write.table(output,"differentusessharedspecies.csv",sep=",")
We can make this file beautiful and searchable if this error is corrected: It looks like row 2 should actually have 3 columns, instead of 4 in line 1.
"Scientific.name","Medicinal.Use.categories","Market"
"832","Sambucus australis ","Piel y tejido subcutáneo","Mercado Rodriguez"
"833","Sambucus peruviana","Sistema músculo esquelético","Calle Santa Cruz"
"834","Sambucus peruviana","Piel y tejido subcutáneo","Calle Santa Cruz"
"835","Sambucus peruviana","Síntomas sin específicar y enfermedades generales","Calle Santa Cruz"
"836","Sambucus peruviana","Piel y tejido subcutáneo","Calle Santa Cruz"
"837","Sambucus peruviana","Síntomas sin específicar y enfermedades generales","Mercado Rodriguez"
"838","Chenopodium quinoa","Sistema músculo esquelético","Calle Santa Cruz"
"839","Dysphania ambrosioides","Sistema digestivo","Mercado Rodriguez"
"840","Dysphania ambrosioides","Sistema digestivo","Mercado Rodriguez"
"841","Dysphania ambrosioides","Sistema digestivo","Mercado Rodriguez"
"842","Dysphania ambrosioides","Síntomas sin específicar y enfermedades generales","Mercado Rodriguez"
"843","Dysphania ambrosioides","Sistema digestivo","Mercado Rodriguez"
"844","Dysphania ambrosioides","Sistema digestivo","El Alto"
"845","Dysphania ambrosioides","Sistema digestivo","El Alto"
"846","Dysphania ambrosioides","Sistema digestivo","Calle Santa Cruz"
"847","Dysphania ambrosioides","Sistema digestivo","Calle Santa Cruz"
"848","Dysphania ambrosioides","Sistema digestivo","Calle Santa Cruz"
"849","Dysphania ambrosioides","Sistema digestivo","Calle Santa Cruz"
"850","Dysphania ambrosioides","Sistema digestivo","Mercado Rodriguez"
"851","Allium sp.","Sangre y sistema circulatorio","Calle Santa Cruz"
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment