<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;
        mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:#0563C1;
        text-decoration:underline;}
span.EmailStyle19
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi all,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I was wondering if anyone had any insights into this?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Cheers,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">David<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="mso-fareast-language:EN-GB">From:</span></b><span lang="EN-US" style="mso-fareast-language:EN-GB"> Fisher, David <david.fisher@abdn.ac.uk>
<br>
<b>Sent:</b> 02 November 2020 10:23<br>
<b>To:</b> statnet_help@u.washington.edu<br>
<b>Subject:</b> [statnet_help] Modelling number of events and event size in bipartite network<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hi folks, <br>
<br>
I posed this query in July, but thought I would ask again with code to perhaps demonstrate my problem a bit better.
<o:p></o:p></p>
<p class="MsoNormal"><br>
I've recently started playing around with modelling actor by event bipartite networks using "ergm". My first mode are actors with a single attribute, which is a 2-level factor (lets call it "class", with reds and blues), while my 2nd mode are events which have
 no attributes. I am interested in modelling 2 processes in a two-mode network that can lead to a higher degree in the projected one-mode network. Actors can get a higher one-mode degree by: 1) attending more events & 2) attending events attended by more others.
 I am interested in specifically modelling these two processes in the two-mode network so I don’t want to project to the one-mode then model a difference in degree.<br>
<br>
I am reasonably happy I get at 1) with ‘b1factor("class")’ <br>
<br>
But right now I cannot work out how to model 2). I am certain that there is a term for it and I am just not understanding the documentation. I don't think it is "b2factor("class")" as events do not vary by class. I don't think it is any of the nodematch functions
 (e.g. "b1nodematch("class")") as I am not interested (for now) in whether reds attend events with other reds. I think I want to compare the number of two paths between the two groups, but I can't see a function to do that. Alternatively, perhaps you could
 frame it as a b2 k-star question, and does the number of different k-stars where the events are the central node depend on the class of those nodes that they are attached too. I’ve tried “b2starmix(2,"class", diff=T)” and this gives 2 different estimates for
 the different classes, which may allow one to say whether the 2 classes are similar in terms of the size of the event, but I am not sure. I have also wondered about adding the event degree as an attribute, and then seeing if reds or blues differ in their preference
 to connect with events differing in that attribute.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Code recreating where I am at using the Davis Southern Women dataset here (I assume you can access the Davis dataset via statnet instead of tnet, but I haven’t been able to track down how; “data(davis)” doesn’t work for me:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">library(tnet)<o:p></o:p></p>
<p class="MsoNormal">library(statnet)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"># Load Davis' (1940) Southern Women Dataset<o:p></o:p></p>
<p class="MsoNormal">data(tnet)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">sw_ebi = as.matrix( table(Davis.Southern.women.2mode) ) #18 women (rows), 14 social events (columns)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">sw_bip = network(sw_ebi, bipartite = T)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">#randomly assign ‘class’ as “red” or “blue”. Code events as “black”<o:p></o:p></p>
<p class="MsoNormal">sw_bip %v% "class" = c(sample(c("red","blue"), nrow(sw_ebi), replace=T), rep("black", ncol(sw_ebi)))
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">plot(sw_bip, vertex.col = "class")<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">#Fit Bipartite ERGM, model differences in number of connections between classes<o:p></o:p></p>
<p class="MsoNormal">ergm_bip = ergm(sw_bip ~ edges + <o:p></o:p></p>
<p class="MsoNormal">                   b1factor("class") , <o:p></o:p></p>
<p class="MsoNormal">                 reference=~Bernoulli)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">summary(ergm_bip) #so this is good, but I also want to model whether blues or reds attend larger events.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><br>
Any help would be greatly appreciated. <br>
<br>
Best regards, <br>
<br>
David<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">David N. Fisher<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">he/him<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">Research Fellow<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">The School of Biological Sciences<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">University of Aberdeen<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB"><a href="https://evoetholab.com/">Web</a> |
<a href="https://scholar.google.com/citations?user=1TuwKzQAAAAJ&hl=en">GS</a> | <a href="https://twitter.com/DFofFreedom">
Tw</a> | <a href="https://www.researchgate.net/profile/David_Fisher30">RG</a> | <a href="https://orcid.org/0000-0002-4444-4450">
Or</a><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB"><br>
<br>
The University of Aberdeen is a charity registered in Scotland, No SC013683.<br>
Tha Oilthigh Obar Dheathain na charthannas clÓraichte ann an Alba, └ir. SC013683.
<o:p></o:p></span></p>
</div>
<br>
<br>
The University of Aberdeen is a charity registered in Scotland, No SC013683.<br>
Tha Oilthigh Obar Dheathain na charthannas clÓraichte ann an Alba, └ir. SC013683.
</body>
</html>