[Rsiena-help] Error message selecting actors

[Maarten Van Zalk]

Dear ruth,

An additional question on selecting specific actors:

As I run the following syntax, I get the error:

 > library(RSiena)
 > library(snow)
 >  library(rlecuyer)
 > setwd('/Users/maartenselfhout/Desktop/CrowdStudy/justin')
 > net.w1 <- as.matrix(read.table("W3SIENA.dat"))
 >  net.w2 <- as.matrix(read.table("W4SIENA.dat"))
 >  net.w3 <- as.matrix(read.table("W5SIENA.dat"))
 >  rad <- as.matrix(read.table("radical.dat"))
 > depr <- as.matrix(read.table("depress.dat"))
 >  self <- as.matrix(read.table("selfesteem.dat"))
 > shy <- as.matrix(read.table("shy.dat"))
 > net.w1[net.w1 %in% c(6,9)] <- NA
 >  net.w2[net.w2 %in% c(6,9)] <- NA
 >  net.w3[net.w3 %in% c(6,9)] <- NA
 > rad[rad %in% c(99)] <- NA
 > depr[depr %in% c(99)] <- NA
 >  self[self %in% c(99)] <- NA
 >  shy[shy %in% c(99)] <- NA
 > use <- rad[, 1] == 1
 > friendship <- sienaNet(array(c(net.w1[use, use], net.w2[use, use],  
net.w3[use, use]),dim=c(796, 796, 3)))
 > shy <- sienaNet(shy[use, ], type='behavior')
 > depr <- varCovar(val = depr[use, ])
 >  self <- varCovar(val = self2[use, ])
Error in varCovar(val = self2[use, ]) : object 'self2' not found
 >  self <- varCovar(val = self[use, ])
 > mymodel <- model.create(useStdInits = TRUE, projname = 'crowd')
 > mydata <- sienaDataCreate(friendship,depr,self,shy,rad)
Error: invalid object in sienaDataCreate matrix
 > mydata <- sienaDataCreate(friendship,depr,self,shy)
 > myeff <- getEffects(mydata)
Error in storage.mode(test) <- "logical" : subscript out of bounds
 >

Do you have any idea why this error arises?

It does NOT arise when I run the syntax without the 'use' command

Thanks!

Kind regards,

Maarten Selfhout-Van Zalk


On Feb 15, 2010, at 12:59 , Ruth M. Ripley wrote:

> Dear Maarten,
>
> There are no special facilities: you would need to extract the  
> subsets of
> the matrices yourself before creating the Siena network object.  
> Using the
> command line, after reading in some data:
>
> friend.data.w1 <- as.matrix(read.table("s50-network1.dat"))
> friend.data.w2 <- as.matrix(read.table("s50-network2.dat"))
> drink <- as.matrix(read.table("s50-alcohol.dat"))
> smoke <- as.matrix(read.table("s50-smoke.dat"))
>
> To select the actors with the first column of smoke equal to 1, say,
>
> use <- smoke[, 1] == 1
>
> friendship <- sienaNet(array(c(friend.data.w1[use, use],
>       friend.data.w2[use, use]), dim=c(38, 38, 2)))
>
> alcohol <- sienaNet(drink[use, ], type="behavior")
>
> To do this via the gui, you would need to do the subsetting outside R.
>
> Regards,
>
> Ruth
>
> PS I have copied this message to the rsiena-help mailing list. I see  
> two
> benefits to using the lists: others might find the information  
> useful too,
> and if I am busy someone else might be able to help!
> --
> Ruth M. Ripley,                         Email:ruth at stats.ox.ac.uk
> Dept. of Statistics,                    http://www.stats.ox.ac.uk/~ruth/
> University of Oxford,                   Tel:   01865 282851
> 1 South Parks Road, Oxford OX1 3TG, UK  Fax:   01865 272595
>
>
> On Mon, 15 Feb 2010, Maarten Van Zalk wrote:
>
>> Dear Ruth,
>>
>> Is it possible in RSiena to select certain actors for the analysis?  
>> So
>> if you have a variable X you want to use to select individuals (for
>> example, all individuals with a score of 1), Can you select only  
>> these
>> individuals?
>>
>> Thanks for your help!
>>
>> Kind regards,
>>
>> Maarten Van Zalk
>>
>