[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