[R] Arrhenius Plot 2 with lattice

[joris meys]

Please provide minimal and self-contained code. Without being able to
read in the dataset, it's impossible to check the code.

Cheers
Joris

On Thu, Nov 26, 2009 at 9:13 AM, Markus Häge <markus_haege at gmx.de> wrote:
> 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 at 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.
>