I am using R on unix version 63.0
I am doing an image plot of the following data file:
===============================lag1 lag2 cif2d
0.000 0.000 NaN
0.000 1.000 0.500000
0.000 2.000 0.489831
0.000 3.000 0.492986
0.000 4.000 0.493409
0.000 5.000 0.492727
0.000 6.000 0.494485
1.000 0.000 0.500000
1.000 1.000 NaN
1.000 2.000 0.495098
1.000 3.000 0.489831
1.000 4.000 0.492986
1.000 5.000 0.493409
1.000 6.000 0.495455
2.000 0.000 0.489831
2.000 1.000 0.495098
2.000 2.000 NaN
2.000 3.000 0.495098
2.000 4.000 0.449153
2.000 5.000 0.283567
2.000 6.000 0.193032
3.000 0.000 0.492986
3.000 1.000 0.489831
3.000 2.000 0.495098
3.000 3.000 NaN
3.000 4.000 0.436275
3.000 5.000 0.284746
3.000 6.000 0.272545
4.000 0.000 0.493409
4.000 1.000 0.492986
4.000 2.000 0.449153
4.000 3.000 0.436275
4.000 4.000 NaN
4.000 5.000 0.436275
4.000 6.000 0.437288
5.000 0.000 0.492727
5.000 1.000 0.493409
5.000 2.000 0.283567
5.000 3.000 0.284746
5.000 4.000 0.436275
5.000 5.000 NaN
5.000 6.000 0.490196
6.000 0.000 0.494485
6.000 1.000 0.495455
6.000 2.000 0.193032
6.000 3.000 0.272545
6.000 4.000 0.437288
6.000 5.000 0.490196
6.000 6.000 NaN
=========================
I read them into a data frame and do the plot.
data<-read.table("~/r/rt/data/unif/6cbcif2d.out",header=TRUE)
x<-unique(data$lag1)
y<-unique(data$lag2)
z<-matrix(data$cif2d,length(y),length(x))
image(x,y,z, col=gray(0:19/19),xlab="Lag 1 (s)", ylab="Lag 2
(s)")
Then I did:> z
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 17 16 8 11 12 10 13
[2,] 16 17 14 8 11 12 15
[3,] 8 14 17 14 7 3 1
[4,] 11 8 14 17 5 4 2
[5,] 12 11 7 5 17 5 6
[6,] 10 12 3 4 5 17 9
[7,] 13 15 1 2 6 9 17
I find this quite strange, since I expected to find numbers from the cif2d
column in the data table above. Is there something funny going on with
matrix()?
Thanks very much for any help.
Bill Simpson
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
