Data_Import_Howland.RMD

#installed and loaded “here” & “palmerpenguins”

#install.packages("here")

#install.packages("palmerpenguins")

#Reading the data

library(here)

## here() starts at /Users/zacharyhowland/OneDrive/intro_quantitative_ecology

library(palmerpenguins)

grazing_data = read.csv(here("data", "grazing_data.csv"))

grazing_dat <- data.frame(grazing_data)

head(grazing_dat)

##   X abundance replicate grass pasture
## 1 1         9         1 short   upper
## 2 2        11         2 short   upper
## 3 3         6         3 short   upper
## 4 4        14         1   med   upper
## 5 5        17         2   med   upper
## 6 6        19         3   med   upper

#By Name

print(grazing_dat$abundance)

##  [1]  9 11  6 14 17 19 28 31 32  7  6  5 14 17 15 44 38 37

#By Position

print(grazing_dat[, 1])

##  [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18

#Second & Fourth Columns At Once

print(grazing_dat[, c(2, 4)])

##    abundance grass
## 1          9 short
## 2         11 short
## 3          6 short
## 4         14   med
## 5         17   med
## 6         19   med
## 7         28  tall
## 8         31  tall
## 9         32  tall
## 10         7 short
## 11         6 short
## 12         5 short
## 13        14   med
## 14        17   med
## 15        15   med
## 16        44  tall
## 17        38  tall
## 18        37  tall

Subset the data into only the penguins that were measured on Torgersen island

require(palmerpenguins)
penguins = data.frame(penguins)
torgersen_penguins <- subset(penguins, island == "Torgersen")

Plot a histogram of penguin body masses on Torgersen Island

torger_penguins = subset(penguins, island == "Torgersen")
hist(
  x = torger_penguins$body_mass_g,
  main = "Body mass of penguins on Torgersen Island",
  xlab = "body mass (g)"
)

# Plot a histogram of sheep abundance from grazing_dat

hist(x = grazing_dat$abundance, main = "Sheep Abundance", xlab = "Sheep Abundance (total ewes)")