are these the same?

Dear members,

I was trying to simulate W which is iid and depends on X and Y. 

Here are 2 methods:

Method 1:
x<-rnorm(100)
y<-rnorm(100)
w<-rnorm(100, 2x+y,1)

Method 2:
x<-rnorm(100)

y<-rnorm(100)

w<-2x+y+rnorm(100,0,1)

Are these methods comparable? 

Since x and y are vectors, the term 2x+y would return a vector. It seems that,

would come from a different population with means depending on the value of x

Method 1 equally acceptable as Method 2 in creating W?

I'd really appreciate clues to this puzzle! Thanks!

set.seed(1001); w1 <- rnorm(100,2*x+y)
set.seed(1001); w2 <- 2*x+y+rnorm(100)
var(cbind(x,y,w1))

var(cbind(x,y,w2))