Fuel Cell 1

getwd() 

## [1] "/Users/michaelsczahor/Desktop/EISData"

setwd("/Users/michaelsczahor/Desktop/EISData")
fc1pola=read.csv("FC1POLA.csv")
fc2pola=read.csv("FC2POLA.csv")
fc1pola$Cell="1"
fc2pola$Cell="2"
FCPOLA=rbind(fc1pola,fc2pola)
library(ggplot2)

head(subset(FCPOLA,Cell="1"))

##     U1v  U2v   U3v   U4v   U5v Ustackv    Ia Japercm Timeh Cell
## 1 0.619 0.61 0.611 0.606 0.602   3.048 100.3   1.003     0    1
## 2 0.619 0.61 0.611 0.606 0.602   3.048 100.2   1.002     0    1
## 3 0.620 0.61 0.611 0.606 0.603   3.050 100.2   1.002     0    1
## 4 0.619 0.61 0.611 0.606 0.602   3.048 100.2   1.002     0    1
## 5 0.619 0.61 0.611 0.606 0.602   3.048 100.2   1.002     0    1
## 6 0.619 0.61 0.611 0.607 0.603   3.050 100.1   1.001     0    1

c1pola=subset(FCPOLA,Cell="1")
c1pola$Timeh=as.factor(c1pola$Timeh)
uu=which(c1pola$Ia<25)
c1pola=c1pola[-uu,]
qplot(Japercm,Ustackv,data=c1pola,colour=Timeh,geom="line")

aa=subset(c1pola,Timeh==0)
bb=subset(c1pola,Timeh==48)
cc=subset(c1pola,Timeh==185)
dd=subset(c1pola,Timeh==348)
ee=subset(c1pola,Timeh=515)
ff=subset(c1pola,Timeh==658)
gg=subset(c1pola,Timeh==823)
hh=subset(c1pola,Timeh==991)

m1=lm(Ustackv~Japercm,data=aa)
m2=lm(Ustackv~Japercm,data=bb)
m3=lm(Ustackv~Japercm,data=cc)
m4=lm(Ustackv~Japercm,data=dd)
m5=lm(Ustackv~Japercm,data=ee)
m6=lm(Ustackv~Japercm,data=ff)
m7=lm(Ustackv~Japercm,data=gg)
m8=lm(Ustackv~Japercm,data=hh)

t0=predict(m1,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t48=predict(m2,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t185=predict(m3,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t348=predict(m4,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t515=predict(m5,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t658=predict(m6,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t823=predict(m7,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
t991=predict(m8,data.frame(Japercm=c(.3,.4,.5,.6,.7,.8,.9,1)))
fc1pol=rbind(t0,t48,t185,t348,t515,t658,t823,t991)
colnames(fc1pol)=c(".3",".4",".5",".6",".7",".8",".9","1")
Timehours=c(0,48,185,348,515,658,823,991)

plot(Timehours,fc1pol[,1],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .3 amps/cm^2")

plot(Timehours,fc1pol[,2],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .4 amps/cm^2")

plot(Timehours,fc1pol[,3],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .5 amps/cm^2")

plot(Timehours,fc1pol[,4],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .6 amps/cm^2")

plot(Timehours,fc1pol[,5],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .7 amps/cm^2")

plot(Timehours,fc1pol[,6],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .8 amps/cm^2")

plot(Timehours,fc1pol[,7],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at .9 amps/cm^2")

plot(Timehours,fc1pol[,8],ylab="Voltage",xlab="Time in Hours",main="FC1 Polarization Data at  1 amp/cm^2")