Hello,
sorry for incomplete code...
with this I read the file and calculate my stuff. I have a plenty of
them 80-300K every 5 Kelvin. I start with 79K, 80K...300K
test<-read.table("T300_both.txt",header=FALSE,sep="")
RH2<-c(RH2,2.5e7*.32e-4/100e-6/5100*(test$V3[c(2)]-test$V3[c(1)]+test
$V3[c(4)]-test$V3[c(3)]))
n2<-c(n2,-1/(RH2[c(46)]*1.60217e-19))
PA2<-c(1.1331/100e-6*.32e-4*(test$V1[c(2)]+test$V1[c(4)]-test
$V1[c(1)]-test$V1[c(3)]))
PB2<-c(1.1331/100e-6*.32e-4*(test$V1[c(6)]+test$V1[c(8)]-test
$V1[c(5)]-test$V1[c(7)]))
P2<-c(P2,(PA2+PB2)/2)
bew2<-c(bew2,-RH2[c(46)]/P2[c(46)])
QA2<-c(QA2,(test$V1[c(2)]-test$V1[c(1)])/(test$V1[c(4)]-test$V1[c(3)]))
QB2<-c(QB2,(test$V1[c(6)]-test$V1[c(5)])/(test$V1[c(8)]-test$V1[c(7)]))
#.......
Temp<-c(79,80,85,90,95,100,106,110,115,120,125,132,135,140,145,151,156,160,165,170,175,180,185,190,195,200,206,210,216,220,225,230,235,240,247,250,255,261,265,270,275,280,285,290,295,300)
#Here I transform my data to make this Arrhenius Plot:
Temp2<-c(1000/Temp)
n2log<-c(log10(n2))
#This function should be plottet as log10(n)~1000/T into the data:
kurve2<-function(T)2^-.5*sqrt(1e17*2*6*I(T^1.5)*(2*pi*0.689*9.109e-31*1.38e-23/43.9e-68)^1.5)*exp(-.047*1.6e-19/(2*1.38e-23*T))
#Here is the plot-function, but I don't know the argument for
#panel.curve()
trellis.device(jpeg,file="Tabh1.jpg",color=FALSE)
n2plot<-xyplot(n2log~Temp2,xlab=list(label="1000/T
(1/K)",fontsize=20),ylab=list(label="log(n/cm^3)",fontsize=20),scales=list(cex=1.5),panel=function(){panel.xyplot(Temp2,n2log);panel.curve(kurve2)})
print(n2plot)
dev.off()
thanks
Markus
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
