[R] matrix weirdness
Bill Simpson
wsimpson at gcal.ac.uk
Thu May 6 12:38:11 CEST 1999
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list