lourymigliorelli — Oct 16, 2013, 5:41 PM
env <- read.csv("rs_env_plot.csv")
Slog <- read.csv("Site_log.csv", header=TRUE)
install.packages("ggplot2")
Error: trying to use CRAN without setting a mirror
require(ggplot2)
Loading required package: ggplot2
install.packages("knitr")
Error: trying to use CRAN without setting a mirror
require(knitr)
install.packages("vegan")
Error: trying to use CRAN without setting a mirror
require(vegan)
Loading required package: vegan
Loading required package: permute
This is vegan 2.0-8
head(Slog)
Man1Rank Man1Log KauRank KauLog Man2Rank Man2Log HonRank HonLog KalRank
1 1 2.2095 1 1.748 1 1.9494 1 1.9085 1
2 2 1.5682 2 1.699 2 1.3617 2 1.3424 2
3 3 1.5051 3 0.000 3 1.2553 3 1.1761 3
4 4 0.8451 4 0.000 4 1.7782 4 0.9542 4
5 5 0.8451 5 0.000 5 0.7782 5 0.9542 5
6 6 0.6990 6 0.000 6 0.7782 6 0.9031 6
KalLog KamRank KamLog Lan1Rank Lan1Log Lan2Rank Lan2Log KaneRank KaneLog
1 2.3856 1 1.3979 1 1.9395 1 1.8195 1 1.7324
2 1.3424 2 0.6990 2 1.7160 2 1.6435 2 1.6812
3 1.3222 3 0.6021 3 1.4472 3 1.0414 3 1.4472
4 1.1139 4 0.3010 4 0.7782 4 0.6021 4 0.6021
5 0.8451 5 0.3010 NA NA 5 0.3010 5 0.0000
6 0.7782 6 0.3010 NA NA 6 0.3010 NA NA
KannRank KannLog LihRank LihLog WaiRank WaiLog NualRank NualLog AwaaRank
1 1 1.146 1 2.149 1 1.7559 1 1.820 1
2 2 1.000 2 1.633 2 1.5563 2 1.785 2
3 3 0.301 3 0.000 3 1.4472 3 1.748 3
4 4 0.301 NA NA 4 1.2041 4 1.431 4
5 5 0.301 NA NA 5 1.1461 5 1.380 5
6 6 0.301 NA NA 6 0.8451 6 1.114 6
AwaaLog MilRank MilLog
1 1.623 1 2.1139
2 1.580 2 1.9494
3 1.568 3 1.7076
4 1.362 4 1.1761
5 1.301 5 0.9542
6 1.255 6 0.9031
ggplot(Slog, aes(x=Slog$Man1Rank, y=Slog$Man1Log)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.2288418", x = 7.5, y = 2)
Warning: Removed 11 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$KauRank, y=Slog$KauLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.3086609", x = 5, y = 1.5)
Warning: Removed 14 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$Man2Rank, y=Slog$Man2Log)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.1485190", x = 10, y = 1.7)
Warning: Removed 7 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$HonRank, y=Slog$HonLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity")+
annotate("text", label = "Slope = -0.1455042", x = 8, y = 1.8)
Warning: Removed 9 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$KalRank, y=Slog$KalLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.1243887", x = 8, y = 2.2)
Warning: Removed 8 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$KamRank, y=Slog$KamLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.1889546", x = 5, y = 1.2)
Warning: Removed 14 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$Lan1Rank, y=Slog$Lan1Log)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.3752949", x = 3.5, y = 1.8)
Warning: Removed 17 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$Lan2Rank, y=Slog$Lan2Log)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.2114740", x = 6, y = 1.6)
Warning: Removed 11 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$KaneRank, y=Slog$KaneLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.4543969", x = 4.5, y = 1.6)
Warning: Removed 16 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$KannRank, y=Slog$KannLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.1727259", x = 6, y = 1)
Warning: Removed 14 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$WaiRank, y=Slog$WaiLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.2093021", x = 7.5, y = 1.6)
Warning: Removed 11 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$NualRank, y=Slog$NualLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.1620015", x = 10, y = 1.8)
Warning: Removed 8 rows containing missing values (position_stack).
ggplot(Slog, aes(x=Slog$AwaaRank, y=Slog$AwaaLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.0916270", x = 18, y = 1.5)
ggplot(Slog, aes(x=Slog$MilRank, y=Slog$MilLog)) +
ylab("Log10(Abundance)") +
xlab("Abundance rank") + geom_histogram(stat="identity") +
annotate("text", label = "Slope = -0.2300050", x = 7.5, y = 2)
Warning: Removed 12 rows containing missing values (position_stack).
par(mfrow=c(3,4))