Panakwate Research Page

By: Cleve Davis

Welcome to the Panakwate Research Page. This site is dedicated to showcasing some of my research and development. As such, everything here should be considered underdevelopment. Please email me if you have any suggestions or ideas of improvement, especially with R coding.  I also want to thank Dr. Brian Dennis and Dr. Dylan Beaudette who have helped me with the programming. Also, check out my www.waamuuza.com page. This has video clips of movies I made for educational purposes and/or just for fun. A link to literature cited can be found here.

Be sure to cite this page!

Thanks!

Soil-Climate Diagram for Moscow, Idaho area

If you remember my old climate diagram take a look the one below! This one is also based upon real data. Here is a link to the code.

Peyote Sales and Button R Graphic (Dawn Dee Williams 2011)

Link to data.

Nobody does artistic drawings no more, and I think we need to revitalize artistic rendering. I also wish I had some dollars to hire some real artists back home! I'm always in too big of rush to make the drawings good.

 

 

R Bar Plot

#Shannon-Wiener Diversity Indices
mydata <- read.table("C:/Users/Cleve/Documents/R Code/Landscape Ecology/SW_Diversity.csv", header=T, sep=",")
attach (mydata)
par(mfrow=c(1,2))
barplot(Early, main="Shannon-Wiener Diversity Index (Early Settlement)",cex.main=.65, col="red", ylim=c(0,1.2), cex.names=.7)
barplot(Late, main="Shannon-Wiener Diversity Index (Late Settlement)", cex.main=.65, col="green", ylim=c(0,1.2), cex.names=.7)

#Shannon-Wiener Evenness Indices

mydata <- read.table("C:/Users/Cleve/Documents/R Code/Landscape Ecology/SW_Evenness.csv", header=T, sep=",")
attach (mydata)
par(mfrow=c(1,2))
barplot(Early, main="Shannon-Wiener Evenness Index (Early Settlement)", cex.main=.65, col="red", ylim=c(0,1), cex.names=.7)
barplot(Late, main="Shannon-Wiener Evenness Index (Late Settlement)", cex.main=.65, col="green", ylim=c(0,1), cex.names=.7)

#Patches

mydata <- read.table("C:/Users/Cleve/Documents/R Code/Landscape Ecology/Patches_Early.csv", header=T, sep=",")
attach (mydata)
par(mfrow=c(1,2))
barplot(Early, main="Number of Patches (Early Settlement)", cex.main=1, names.arg=c("Forest", "Grassland", "Agriculture"), col="red", ylim=c(0,20), cex.names=.7)
mydata1 <- read.table("C:/Users/Cleve/Documents/R Code/Landscape Ecology/Patches_Late.csv", header=T, sep=",")
attach (mydata1)

R Pie Plot

Landscape Ecology Pie Charts
#Early Proportion of Pixels
slices <- c(59, 28, 13)
par(mfrow=c(2,1))
lbls <- c("Forest", "Grassland", "Agriculture")
pct <- round(slices/sum(slices)*100)
lbls <- paste(lbls, pct) # add percents to labels
lbls <- paste(lbls,"%",sep="") # ad % to labels
pie(slices,labels = lbls, col=rainbow(length(lbls)), main="Proportion of Pixels (Early Settlement)")
slices <- c(39, 32, 29)
slices1=slices
lbls <- c("Forest", "Grassland", "Agriculture")
pct <- round(slices1/sum(slices1)*100)
lbls <- paste(lbls, pct) # add percents to labels
lbls <- paste(lbls,"%",sep="") # ad % to labels
pie(slices1,labels = lbls, col=rainbow(length(lbls)), main="Proportion of Pixels (Late Settlement)")