Re: [R] loop to subtract arrays / error

Irucka Embry Mon, 19 Nov 2012 05:59:00 -0800

Hi Rui, thank-you. 

That was simple and worked great.


Irucka 

<-----Original Message-----> 
>From: Rui Barradas [ruipbarra...@sapo.pt]
>Sent: 11/19/2012 4:13:24 AM
>To: iruc...@mail2world.com
>Cc: r-help@r-project.org
>Subject: Re: [R] loop to subtract arrays / error
>
>Hello,
>
>Or simpler, since Vobsr only has one column:
>
>Xjj <- as.vector(Vobsr) - Vsimr
>
>Hope this helps,
>
>Rui Barradas
>Em 19-11-2012 10:05, Rui Barradas escreveu:
>> Hello,
>>
>> Try the following.
>>
>> Xjj <- matrix(nrow = 101, ncol = 1000)
>> for (i in 1:dim(Vsimr)[2]) {
>> Xjj[, i] <- Vobsr - Vsimr[, i]
>> }
>>
>> Hope this helps,
>>
>> Rui Barradas
>> Em 19-11-2012 01:41, iembry escreveu:
>>> Hi everyone, I am having trouble with creating a loop to subtract 
>>> arrays.
>>>
>>> In R, this is what I have done:
>>>> Vobsr <- read.csv("Observed_Flow.csv", header = TRUE, sep =",") # 
>>>> see data
>>>> below
>>>> Vsimr <- read.csv("1000Samples_Vsim.csv", header = TRUE, sep
>=",") # 
>>>> see
>>>> data below
>>>> Vobsr <- as.matrix(Vobsr[,-1]) # remove column 1 from analysis thus
>>>> Vobsr
>>>> is 101x1 double matrix (column 1 is date information)
>>>> Vsimr <- as.matrix(Vsimr[,-1]) # remove column 1 from analysis thus
>>>> Vsimr
>>>> is 101x1000 double matrix (column 1 is date information)
>>>> Vobsr - Vsimr
>>> Error in Vobsr - Vsimr : non-conformable arrays
>>>
>>> Thus I attempted to create the loop below to perform the subtraction
>>> operation for each of the 1000 columns.
>>>
>>>> dim(Vsimr)[2]
>>> [1] 1000
>>>
>>>> for (i in 1:dim(Vsimr)[2]) {
>>> Xjj <- Vobsr - Vsimr[,i]
>>> }
>>> Xjj is a 101x1 double matrix rather than a 101X1000 double matrix
>>>
>>> How can I subtract each column of Vsimr from the single column of 
>>> Vobsr over
>>> the 1000 columns present?
>>>
>>> I would like to thank each of you in advance for your assistance.
>>>
>>> How can I subtract each column of Vsimr from the single column of 
>>> Vobsr over
>>> the 1000 columns present?
>>>
>>> I would like to thank each of you in advance for your assistance.
>>>
>>>
>>> I am including some of the data from the files that I am operating
on 
>>> below:
>>> 1 column of Observed_Flow.csv
>>> 81.071
>>> 73.187
>>> 66.991
>>> 62.482
>>> 59.662
>>> 58.529
>>> 59.085
>>> 61.328
>>> 65.259
>>> 70.878
>>> 78.184
>>> 87.179
>>> 97.862
>>> 110.23
>>> 124.29
>>> 140.08
>>> 157.57
>>> 176.76
>>> 197.63
>>> 220.18
>>> 244.4
>>> 270.31
>>> 297.88
>>> 327.14
>>> 358.09
>>> 390.71
>>> 425.03
>>> 461.03
>>> 498.72
>>> 538.09
>>> 579.16
>>> 621.91
>>> 666.35
>>> 712.48
>>> 760.29
>>> 809.8
>>> 860.99
>>> 913.87
>>> 968.44
>>> 1024.7
>>> 1082.6
>>> 1142.3
>>> 1203.6
>>> 1266.6
>>> 1331.3
>>> 1397.7
>>> 1465.7
>>> 1535.5
>>> 1606.9
>>> 1680.1
>>> 1754.9
>>> 1831.4
>>> 1907.1
>>> 1981.9
>>> 2055.9
>>> 2129
>>> 2201.3
>>> 2272.7
>>> 2343.3
>>> 2413.1
>>> 2482
>>> 2550.1
>>> 2617.3
>>> 2683.7
>>> 2749.2
>>> 2813.9
>>> 2877.8
>>> 2940.8
>>> 3003
>>> 3064.3
>>> 3124.8
>>> 3184.4
>>> 3243.2
>>> 3301.1
>>> 3358.2
>>> 3414.5
>>> 3469.9
>>> 3524.4
>>> 3578.2
>>> 3631
>>> 3683.1
>>> 3734.3
>>> 3784.6
>>> 3834.1
>>> 3882.8
>>> 3930.6
>>> 3977.6
>>> 4023.7
>>> 4069
>>> 4113.4
>>> 4157
>>> 4199.8
>>> 4241.7
>>> 4282.7
>>> 4323
>>> 4362.3
>>> 4400.9
>>> 4438.6
>>> 4475.4
>>> 4511.4
>>> 4546.6
>>>
>>>
>>> 2 columns of 1000 columns of 1000Samples_Vsim.csv
>>> 81.07 81.07
>>> 73.19 73.19
>>> 65.81 67.16
>>> 58.93 63
>>> 52.55 60.7
>>> 46.68 60.25
>>> 41.31 61.67
>>> 36.44 64.95
>>> 32.08 70.08
>>> 28.22 77.08
>>> 24.86 85.94
>>> 22.01 96.65
>>> 19.65 109.23
>>> 17.8 123.67
>>> 16.46 139.96
>>> 15.61 158.12
>>> 15.27 178.14
>>> 15.43 200.02
>>> 16.1 223.75
>>> 17.27 249.35
>>> 18.94 276.81
>>> 21.11 306.13
>>> 23.79 337.31
>>> 26.97 370.34
>>> 30.65 405.24
>>> 34.84 442
>>> 39.52 480.62
>>> 44.71 521.1
>>> 50.41 563.44
>>> 56.61 607.64
>>> 63.31 653.7
>>> 70.51 701.62
>>> 78.21 751.4
>>> 86.42 803.04
>>> 95.13 856.53
>>> 104.35 911.89
>>> 114.06 969.11
>>> 124.28 1028.2
>>> 135.01 1089.1
>>> 146.23 1151.9
>>> 157.96 1216.6
>>> 170.19 1283.1
>>> 182.93 1351.5
>>> 196.16 1421.7
>>> 209.9 1493.8
>>> 224.15 1567.8
>>> 238.89 1643.6
>>> 254.14 1721.3
>>> 269.89 1800.8
>>> 286.15 1882.2
>>> 302.91 1965.5
>>> 320.17 2050.6
>>> 337.18 2134.8
>>> 353.93 2218.1
>>> 370.44 2300.4
>>> 386.69 2381.8
>>> 402.7 2462.3
>>> 418.45 2541.8
>>> 433.95 2620.4
>>> 449.2 2698.1
>>> 464.2 2774.9
>>> 478.94 2850.7
>>> 493.44 2925.6
>>> 507.68 2999.5
>>> 521.67 3072.6
>>> 535.41 3144.7
>>> 548.9 3215.8
>>> 562.14 3286.1
>>> 575.12 3355.4
>>> 587.86 3423.8
>>> 600.34 3491.2
>>> 612.57 3557.7
>>> 624.55 3623.3
>>> 636.28 3688
>>> 647.76 3751.7
>>> 658.98 3814.5
>>> 669.96 3876.4
>>> 680.68 3937.3
>>> 691.15 3997.3
>>> 701.37 4056.4
>>> 711.34 4114.6
>>> 721.06 4171.8
>>> 730.52 4228.1
>>> 739.74 4283.4
>>> 748.7 4337.9
>>> 757.41 4391.4
>>> 765.87 4443.9
>>> 774.08 4495.6
>>> 782.04 4546.3
>>> 789.74 4596.1
>>> 797.2 4644.9
>>> 804.4 4692.8
>>> 811.35 4739.8
>>> 818.05 4785.9
>>> 824.5 4831
>>> 830.7 4875.2
>>> 836.64 4918.5
>>> 842.33 4960.8
>>> 847.78 5002.2
>>> 852.97 5042.7
>>> 857.91 5082.3
>>>
>>>
>>>
>>> -- 
>>> View this message in context: 
>>>
http://r.789695.n4.nabble.com/loop-to-subtract-arrays-error-tp4650001.ht
ml
>>> Sent from the R help mailing list archive at Nabble.com.
>>>
>>> ______________________________________________
>>> 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.
>>
>> ______________________________________________
>> 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.
>
>.
> 


<span id=m2wTl><p><font face="Arial, Helvetica, sans-serif" size="2" 
style="font-size:13.5px">_______________________________________________________________<BR>Get
 the Free email that has everyone talking at <a href=http://www.mail2world.com 
target=new>http://www.mail2world.com</a><br>  <font color=#999999>Unlimited 
Email Storage &#150; POP3 &#150; Calendar &#150; SMS &#150; Translator &#150; 
Much More!</font></font></span>
        [[alternative HTML version deleted]]

______________________________________________
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.

Re: [R] loop to subtract arrays / error

Reply via email to