[R] Error in Fields TPS function {svd ...} again
Jason Bond
jbond at arg.org
Wed Oct 30 21:34:57 CET 2002
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 --------------
More information about the R-help
mailing list