Thanks for all the helpful responses. I include the data file and the
syntax file for reference. Again, if I use the fields function, as is, I
get the message:
Error in svd(tempM) : error 159 in dsvdc
using traceback, I get:
> traceback()
4: stop(paste("error ", z$info, " in dsvdc"))
3: svd(tempM)
2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp, scale.type =
scale.type,
outputcall = Tpscall, p = p, ...)
1: Tps(bvolcap, bdsm)
>
if I change the occurrence of svd in the fields package to La.svd, I get
the error message
> bout <- Tps( bvolcap, bdsm)
Error in "[<-"(*tmp*, (nt + 1):np, (nt + 1):np, value = temp$v) :
number of items to replace is not a multiple of replacement length
>
traceback() then gives:
> traceback()
2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp, scale.type =
scale.type,
outputcall = Tpscall, p = p, ...)
1: Tps(bvolcap, bdsm)
>
Sorry to be so useless, I am a R (S) newbie. Any advice is greatly appreciated.
-------------- next part --------------
data <- read.table("splus-data-black.dat")
bdsmt <- data[1]
bdsm <- bdsmt[bdsmt > -777]
bvolcapt <- data[2]
bvolcap <- bvolcapt[bvolcapt > -777]
bout <- Tps( bvolcap, bdsm)
-------------- next part --------------
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
2 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
2 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 1.5
0 3
1 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3
0 3.5
0 5
1 5
0 5
0 5
0 5
1 5
0 6
1 6.5
0 6.5
0 6.5
0 6.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
1 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
1 7.5
1 7.5
1 7.5
1 7.5
0 7.5
1 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 7.5
0 8.5
0 8.5
0 10
0 10
0 10
0 10
0 10
0 11
1 11
0 11
0 11
1 11
0 12.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
1 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
1 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 13.5
0 14.5
0 14.5
0 14.5
0 14.5
0 14.5
0 14.5
0 14.5
1 17
0 17
0 17.5
0 17.5
1 18.5
0 20.5
1 20.5
0 20.5
1 22
1 22
2 22
0 23
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 25
0 26.5
5 26.5
1 30
0 31
2 31
1 31
2 31
0 31
0 31
1 31
1 36
0 36
0 36
1 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
1 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
2 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 36
0 37.5
0 39
0 39.5
0 39.5
0 39.5
0 39.5
0 42.5
0 43
0 43
0 44.5
0 45
0 45
0 45
1 46.5
0 51.5
0 53.5
0 53.5
0 53.5
0 53.5
2 55
0 59.5
0 59.5
3 61
0 61
1 65.5
2 65.5
0 67.5
0 67.5
0 67.5
0 67.5
0 67.5
0 68.5
0 73
0 73.5
1 74
0 82.5
0 85
0 87
0 91.5
0 97.5
0 97.5
0 97.5
3 97.5
0 99
0 102.5
0 104.5
2 109.5
2 117
0 120
0 120
0 120
0 120
0 120
0 120
0 120
4 120
0 120
1 120
1 120
2 120
0 120
0 120
0 120
0 120
0 120
0 121.5
0 126
0 126
0 126
0 126
0 126
1 126
0 126
0 126
0 126
0 126
2 126
0 126
0 126
0 126
0 126
0 126
0 126
0 126
2 126
0 126
0 126
0 126
0 126
1 126
4 126
0 126
1 126
0 126
0 126
0 126
0 126
0 126
0 126
0 126
0 126
0 126
0 126
0 126
1 126
0 126
0 126
0 127.5
0 129.5
0 129.5
0 129.5
2 131
0 132
0 133
0 133
0 133
0 133
0 135.5
0 143.5
0 148.5
0 149.5
0 150
0 150
0 155.5
0 155.5
0 155.5
0 157.5
0 157.5
2 157.5
0 157.5
0 157.5
0 159.5
0 169
3 170.5
0 172.5
0 173.5
1 173.5
0 180
0 183.5
0 184.5
0 184.5
1 187.5
0 187.5
0 191
0 204.5
1 207
0 210
2 210
1 210
0 210
0 223.5
1 240
0 242.5
0 264
3 264
0 264
0 270
0 270
0 270
0 270
0 270
0 270
0 273.5
0 275.5
0 287.5
4 287.5
0 297
1 301.5
5 324
0 330
1 330
0 330
5 376.5
1 401.5
0 420
1 420
0 420
1 420
0 420
0 420
1 420
0 420
1 420
1 420
0 430.5
0 432
0 441.5
5 450
5 450
1 456
1 464.5
0 469.5
0 492
0 492
0 492
0 498
0 511
0 521
0 540
4 540
0 540
0 540
0 540
0 540
0 540
0 543.5
0 547
3 557
0 557.5
4 564
0 582
7 606
2 609
0 666
1 666
0 692.5
1 708
2 715.5
0 756
0 792
1 792
1 804
0 805.5
0 834
6 863
5 900
2 900
6 924
0 924
0 924
0 943.5
1 969
0 1015.5
5 1053
2 1068
0 1068
2 1088
4 1126.5
5 1195.5
0 1267.5
1 1267.5
0 1267.5
2 1285.5
0 1302
0 1357.5
0 1405.5
1 1416
5 1487.5
0 1487.5
0 1512
1 1717.5
0 1722
1 1722
3 1739.5
3 1969
1 2041.5
4 2099.5
3 2121
0 2244
0 2340
2 2352
3 2519
1 2688
0 2807.5
1 2807.5
1 2807.5
2 2807.5
3 3450
1 4080
0 4680
7 4680
0 4680
1 4680
-------------- next part --------------
> Thanks for all the helpful responses. I include the data file and the > syntax file for reference. Again, if I use the fields function, as is, >It helped to have your example:> data <- read.table("splus-data-black.dat") > dim(data)[1] 606 2> summary(data)V1 V2 Min. :0.0000 Min. : 0.0 1st Qu.:0.0000 1st Qu.: 3.0 Median :0.0000 Median : 36.0 Mean :0.3977 Mean : 247.6 3rd Qu.:0.0000 3rd Qu.: 149.9 Max. :7.0000 Max. :4680.0> table(data[1])0 1 2 3 4 5 6 7 488 64 24 10 7 9 2 2 Tps(x, Y, ...) expects x to be a matrix, in the geostatistical "world view" of the locations of the observations in say 2d or 3d:> Tps(data[,2], data[,1])Error in svd(tempM) : error 159 in dsvdc The distribution of your Y makes it look like counts, both x and Y have smaller numbers of unique values than n=606:> length(unique(data[,1]))[1] 8> length(unique(data[,2]))[1] 161> plot(data[,2], data[,1])shows that the data are rather different from the ones in example(Tps). So this doesn't look like a typical numerical problem, but Tps being used on data for which it was not designed, and failing to trap this state. Roger I get the message:> > Error in svd(tempM) : error 159 in dsvdc > > using traceback, I get: > > > traceback() > 4: stop(paste("error ", z$info, " in dsvdc")) > 3: svd(tempM) > 2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp, scale.type > = scale.type, > outputcall = Tpscall, p = p, ...) > 1: Tps(bvolcap, bdsm) > > > > > if I change the occurrence of svd in the fields package to La.svd, I get > the error message > > > bout <- Tps( bvolcap, bdsm) > Error in "[<-"(*tmp*, (nt + 1):np, (nt + 1):np, value = temp$v) : > number of items to replace is not a multiple of replacement > length > > > > traceback() then gives: > > > traceback() > 2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp, scale.type > = scale.type, > outputcall = Tpscall, p = p, ...) > 1: Tps(bvolcap, bdsm) > > > > Sorry to be so useless, I am a R (S) newbie. Any advice is greatly > appreciated.-- Roger Bivand NHH, Breiviksveien 40, N-5045 Bergen, Norway (travelling but still accessible) -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._