source("commands.s")
ls()
l
ls()
read.table("calorimeter.data")
ls()
aa<-read.table("calorimeter.data")
aa
sink()
aa<-read.table("calorimeter.data")
aa
ls()
rm(aa)
rm(xx.start)
ls()
source("commands.s")
ls()
aa
xx.start
ls()
10<11
(10<11)*2
(10<10)*2
time
t
aa$time
?read.table
attach("aa")
attach(aa)
time
attach(xx.start)
ls()
ls(pos=2)
source("commands.s")
attach(as.data.frame(xx.start))
ls(pos=2)
xx.start
tzero
source("commands.s")
xx.start<-c(tzero=200,base1=30,base2=0,amp=30000,kinit=.02,kdecay=.2)
source("commands.s")
source("commands.s")
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=119,base1=24.7,base2=15.1,amp=41910,kinit=.1188,kdecay=..00832)
xx.start<-c(tzero=119,base1=24.7,base2=15.1,amp=41910,kinit=.1188,kdecay=.00832)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=200,base1=30,base2=0,amp=30000,kinit=.02,kdecay=.2)
xx.start<-c(tzero=200,base1=30,base2=0,amp=30000,kinit=.2,kdecay=.02)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start
xx.start<-c(tzero=100,base1=30,base2=0,amp=30000,kinit=.2,kdecay=.02)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=100,base1=30,base2=10,amp=30000,kinit=.2,kdecay=.02)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=130,base1=30,base2=10,amp=30000,kinit=.2,kdecay=.02)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=130,base1=30,base2=0,amp=30000,kinit=.2,kdecay=.02)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=150,base1=30,base2=0,amp=30000,kinit=.2,kdecay=.02)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start<-c(tzero=120,base1=30,base2=-50,amp=30000,kinit=.2,kdecay=.01)
nls(voltage ~ (time < tzero)*(base1 - base2) + base2 + amp*(exp(-kinit*(time>=tzero)*(time-tzero)) - exp(-kdecay*(time>=tzero)*(time-tzero))), aa,start=xx.start,trace=TRUE)
xx.start
source("commands.s")
xx<-0
xx<-ls()
rm(xx)
ls()
detach()
aa<-ls()
rm(list=xx)
rm(list=aa)
ls()(
)
ls()
source("commands.s")
ls()
summary(aa)
summary(aa.out)
:q
aa.out$fitted
out)
fitted.values(aa.out)
source("commands.s")
source("commands.s")
graphics.off()
sink()
source("commands.s")
source("commands.s")
ls()
aa<-ls()
rm(list=aa)
ls()
source("commands.s")
ls()
a.out
aa.out
var(aa.out)
var(aa)
var(aa$voltage)
sqrt(var(aa$voltage))
?nls
?nlsModel
aa.out$deiance()
aa.out$deviance()
aa.out$deviance
aa.out[[deviance]]
aa.out[deviance]
aa.out$$deviance
aa.out$$deviance()
aa.out$deviance[]
ls()
list(aa.out)
?deviance
deviance(aa.out)
residuals(aa.out)
sum(residuals(aa.out))
sum(residuals(aa.out)^2)
sum(residuals(aa.out)^2)/912
sqrt(sum(residuals(aa.out)^2)/912)
sqrt(sum(residuals(aa.out)^2)/906)
q()
