# 9-8 Exercise 3 - Solar Flares

# The data is held in the file SolarFlares.csv
# This is within the SP09 folder of the current working directory
# The file contains headers

# Load SolarFlares file

solarflares<-read.csv("SP09/SolarFlares.csv",stringsAsFactors = T, header=TRUE)

# Display solarflares

solarflares

# Display headers

names(solarflares)

# Mean number of predicted C-class flares

mean(solarflares$PredictedC)

# Mean number of predicted X-class flares
# for small and large areas

tapply(solarflares$PredictedX,solarflares$Area, mean)

# Create a new column for the total number of flares

totalflares<-solarflares$PredictedC+
  solarflares$PredictedM+
  solarflares$PredictedX

# Bind the total number of flares to solarflares

solarflares<-cbind(solarflares,PredictedFlares=totalflares)

# Display solarflares

solarflares

# Display the contingency table for 
# spot size vs number of predicted flares

split(solarflares$PredictedFlares, solarflares$SpotSize)

# Display the contingency table for 
# area vs number of predicted flares

split(solarflares$PredictedFlares, solarflares$Area)

# Display the contingency table for 
# activity vs area vs number of predicted flares

split(solarflares$PredictedFlares, 
      list(solarflares$Activity,solarflares$Area))