On Mon, Feb 13, 2012 at 02:08:52AM -0800, uday wrote:
>
> I have some dataset
>
>
> sci.pre <- 0.30 0.38 0.50 0.65 0.80 1.03 1.33 1.72
> 2.22
> 2.87 3.81 5.06 6.75 9.00 11.97 14.15 16.34 19.04
> 22.27 25.49 29.72 34.67 40.47 47.29 55.29 64.67 75.6
> 88.50 103.50 121.10 141.70 165.80 194.00 227.00 265.00 308.00
> 356.50 411.00 441.60 472.20 506.35 540.50 578.55 616.60 637.75
> 658.90 680.05 701.20 724.65 748.10 771.55 795.00 820.95 846.90
> 872.85 898.80 927.35 955.90 984.45 1013.00 ( length is 60)
>
> sci.avgkernal < - 0.30 0.38 0.50 0.65 0.80 1.03 1.33
> 1.72 2.22
> 2.87 3.81 5.06 6.75 9.00 11.97 14.15 16.34 19.04
> 22.27 25.49 29.72 34.67 40.47 47.29 55.29 64.67 75.65
> 88.50 103.50 121.10 141.70 165.80 194.00 227.00 265.00 308.00
> 356.50 411.00 441.60 472.20 506.35 540.50 578.55 616.60 637.75
> 658.90 680.05 701.20 724.65 748.10 771.55 795.00 820.95 846.90
> 872.85 898.80 927.35 955.90 984.45 1013.00 ( length is 60)
>
> pres.interptimes <- [,1] [,2] [,3] [,4] [,5]
> [,6] [,7] [,8]
> [1,] 1016.1267 1005.9741 989.9127 970.0237 945.6067 880.5082 790.4647
> 675.8315
> [2,] 875.6320 866.8767 853.0258 835.8741 814.8176 758.6784 681.0275
> 582.1712
> [3,] 996.0351 986.0758 970.3201 950.8098 926.8576 862.9984 774.6692
> 662.2184
> [4,] 996.0353 986.0760 970.3203 950.8100 926.8578 862.9987 774.6694
> 662.2187
> [5,] 1008.0222 997.9431 981.9978 962.2527 938.0123 873.3847 783.9926
> 670.1888
> [6,] 999.8343 989.8371 974.0214 954.4367 930.3932 866.2906 777.6247
> 664.7453
> [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16]
> [1,] 544.3248 410.3611 289.4130 237.2794 191.5622 152.3020 119.1842 91.68080
> [2,] 468.7636 353.2371 248.9350 203.9765 164.5513 130.6945 102.1346 78.41645
> [3,] 533.2153 401.8020 283.1566 232.0154 187.1686 148.6559 116.1686 89.18875
> [4,] 533.2155 401.8022 283.1568 232.0157 187.1688 148.6561 116.1688 89.18898
> [5,] 539.6334 406.6388 286.5657 234.8092 189.4227 150.4466 117.5683 90.26387
> [6,] 535.2505 403.3363 284.2386 232.9026 187.8848 149.2253 116.6141 89.53148
> [,17] [,18] [,19] [,20]
> [1,] 69.12170 50.81654 24.30808 0.8657024
> [2,] 58.96213 43.17629 20.31615 0.1001228
> [3,] 67.05913 49.10247 23.09867 0.1025881
> [4,] 67.05936 49.10270 23.09890 0.1028196
> [5,] 67.86797 49.69524 23.37854 0.1057444
> [6,] 67.31751 49.29240 23.18949 0.1057560 ( dim is 6 20)
>
> sci.prediff <- diff(sci.pre)
> sci.prediff <- c(sci.pre[1],sci.prediff)
> sum(sci.avgkernal*sci.prediff )/sum(sci.prediff )
> pres.interptimes <- pres.interptime[,-20]#skip last level
> tm3.avgkernal <-array(NA,c(length(1:nobs),19))
> for (k in 1:nobs){
> for (h in 1:19){
> sel <- sci.pre <= pres.interptimes[k,h] & sci.pre >
> pres.interptimes[k,h+1]
> tm3.avgkernal[k,h] <- sum((sci.avgkernal * sci.prediff)[sel]) /
> sum(sci.prediff[sel])
> }
> }
>
> after running code I get error
>
> Error: subscript out of bounds
>
> How to fix this error ?
I suspect that the error is generated at
pres.interptimes[k,h+1]
since h+1 goes up to 20 and pres.interptimes was restricted to 19
columns in
pres.interptimes <- pres.interptime[,-20]#skip last level
(assuming that "pres.interptime" is a typo and is "pres.interptimes"
in fact).
Hope this helps.
Petr Savicky.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
