2

在 NetLogo 中建立了一个不同年龄代理的社交网络,类似如下,从而形成一个由链接连接的代理圈。这些链接的总体目的是代表这些链接之间的联系。该模型模拟了感染通过网络的传播。代理开始时是易感的,如果他们与传染性链接邻居接触,就有可能被感染。例如,我想模拟受感染个体的隔离或隔离。即他们与其他人的链接将被完全停用,或者至少他们的大部分链接将被停用。理想情况下,我会按下观察者界面上的按钮来停用受感染代理的链接。我' d 还希望能够模拟关闭学校,例如,大多数幼儿和儿童以及他们之间的联系将被停用,从而阻止儿童之间的感染能力。如果它们数量较少,它们与成年人的联系可能应该保持开放,尽管现在一次只关注一个品种之间的联系可能会更好。

breed [ toddlers toddler ]
breed [ children child ]
breed [ adults adult ]
breed [ over45s over45 ]

globals
[
  num-nodes
  num-infected
  num-susceptible  
  prob-infection-toddler
  prob-infection-child
  prob-infection-adult
  prob-infection-over45
  force-of-infection 
]

turtles-own
[
  temp-infected?    ;; temporary infected state
  infected?         ;; true if agent has been infected
  susceptible?
  num-infected-neighbors
]

links-own
[
 closed? 
]


to setup
  clear-all
  create-toddlers 20
  create-children 20
  create-adults 20
  create-over45s 20
  create-network
  reset-ticks
  layout-circle (sort turtles) max-pxcor - 8
  ask turtles
  [
     facexy 0 0
     if who mod 2 = 0 [fd 4]
  ]
  ask turtles
    [set susceptible? true]

  ask one-of turtles
  [
    set infected? true
    set susceptible? false
    set color red
  ]
  display
end

to create-network

  let connexions (list
    (list toddlers toddlers 5)
    (list toddlers children 2)
    (list toddlers adults 2)
    (list toddlers over45s 1)
    (list children toddlers 3)
    (list children children 8)
    (list children adults 5)
    (list children over45s 1)
    (list adults toddlers 1)
    (list adults children 3)
    (list adults adults 6)
    (list adults over45s 3)
    (list over45s toddlers 1)
    (list over45s children 1)
    (list over45s adults 5)
    (list over45s over45s 5)
  )

  foreach connexions [
    let source-breed item 0 ?
    let target-breed item 1 ?
    let num-connexions item 2 ?
    let reverse-num-connexions item 2 first filter [
      item 0 ? = target-breed and item 1 ? = source-breed
    ] connexions
    ask source-breed [
      repeat num-connexions [
        let possible-targets other target-breed with [
          (not member? myself link-neighbors) and
          (count link-neighbors with [ breed = source-breed ] < reverse-num-connexions)
        ]
        let target one-of possible-targets
        if target != nobody [ create-link-with target ]
      ]
    ]
  ]

  ask links [set closed? false]
end

to spread
  ;;; there is one of these for each age group as they have different probabilities for infection
  ask toddlers with [ susceptible? = true ]
    [

   ;;tried changing to something like this but this is the line that throws up an error  
   ask my-links with [closed? = false][
       set num-infected-neighbors count (link-neighbors with [infected? = true])
   ]
   ;;should only include active links
   set force-of-infection (prob-infection-toddler) ^ num-infected-neighbors ;;paraphrased equation but num-infected-neigbours is important

   if ( random-float 1 <= force-of-infection)  ;; infect with probability p
     [

            set temp-infected? true
            set susceptible? false

     ]
 ]
   ask turtles with [temp-infected? = true]
    [
      set infected? true
      set temp-infected? false
    ]
end

to isolate-infected
   ;;for all infected individuals deactivate their links
end


to close-schools
   ask toddlers 
   [
     ask my-links [set closed? true]
   ]

end

to close-offices
   ;; e.g cut the number of links between adults and so on
end

正如您所看到的,这里有两个问题,一个是根据终端节点所处的状态或品种来停用链接(理想情况下,一旦学校停课结束或一旦受感染的个体恢复,它们可以再次打开)他们是。

第二个问题是计算受感染链接邻居的数量,而不包括停用的链接。隐藏链接并不能阻止它的存在,它只是变得不可见,因此不会停止联系。我还考虑过简单地将停用链接的颜色更改为例如黑色。有没有一种方法可以重写受感染邻居的计数以仅计算非黑色或隐藏的链接,例如set num-infected-neighbors count (link-neighbors with [infected? = true]......并且链接隐藏?= false .... 或链接颜色“!=”黑色?)

我觉得第二个问题可能是两者中较容易的一个,但我可能弄错了。在这个阶段我遇到了太多的砖墙,我的头被这件事炒了。任何帮助将不胜感激。感谢您花时间阅读本文,我意识到这有点咆哮 :) 再次感谢 Nicolas Payette 之前的帮助

编辑:为关闭添加了链接自己的布尔值?

试图更改计算打开链接上受感染邻居数量的部分,但出现错误 this code can't be run by a link error while link 14 15 running SET called by procedure SPREAD called by Button 'Spread'

还在“关闭学校”中添加了一个功能,该功能通过设置关闭来提供所有幼儿链接的非常基本的关闭?为真

编辑:这会是一个可能的解决方案吗

set num-infected-neighbors 0
   ask my-links with [closed? = false][
      ask other-end [if (infected?) [set num-infected-neighbors num-infected-neighbors + 1 ]]
   ]
   set force-of-infection 1 - (1 - prob-infection-adult) ^ num-infected-neighbors

编辑:据我所知,num-infected-neighbors 应该是一个turtles-own变量,但是当我将手表放在乌龟身上并运行模拟时,乌龟的 num-infected-neighbors 似乎一直超过实际数量乌龟拥有的链接邻居。这也是错误的。虽然我不明白为什么....

编辑:

let infectnum 0
  set num-infected-neighbors 0
   ask my-links with [closed? = false][
      ask other-end [if (infected?) [set infectnum infectnum + 1 ]]
       ;;set num-infected-neighbors count (link-neighbors with [infected? = true])
   ]
   set num-infected-neighbors infectnum

似乎也不能正常工作......

编辑:供将来参考问题已解决 -

set num-infected-neighbors length filter [[infected?] of ?]
                        [other-end] of my-links with [not closed?]

首先,由未关闭链接上的链接邻居组成一个列表,然后对其进行过滤以仅提供受感染的链接。num-infected-neighbors 是该列表的长度

4

1 回答 1

2

首先是一些一般性建议:最好给出一个最小工作示例(MWE),即一个被剥离到最低限度以显示问题的示例。这段代码有太多细节,所以人们不太可能回答。此外,有时一旦你把它拆开,你就会找出答案。此代码也不起作用:缺少全局变量的定义。

您可能想查看 NetLogo 附带的模型库中的“网络上的病毒”。它做了类似的事情,但使用了更“NetLogoey”的编码风格。

您可以测试链接的颜色。即,您可以使用链接是“停用”颜色的事实来测试它是否应该被计算或遵循。一个更好的主意是使用 为链接提供一个变量links-own,并对其进行测试。但这里有一些示例代码(根据“网络上的病毒”中的过程修改)为此目的使用颜色:

to spread-virus
  ask turtles with [infected?]
    [ ask my-links with [color != red]
      [ask other-end
        [ if (not resistant?) and (random-float 100 < virus-spread-chance)
            [ become-infected ] ] ] ]
end

诀窍是,我不是直接看 ,而是link-neighbors向每只乌龟询问它的my-links,然后测试它们的颜色。然后对于通过测试的每个链接,我向它询问它的other-end,即不是通过 找到链接的海龟的海龟my-links。然后病毒被传播给那只乌龟——但前提是链接没有被禁用。(我认为这段代码工作正常——似乎。但如果你做这样的事情,你应该彻底测试它。)使用计数禁用或非禁用链接my-links应该更容易。因为isolate-infected你可以做一些类似ask turtles with [infected?] [ask my-links [set颜色/停用/等的事情。]]. 没有测试过。

对于响应我的回答的编辑版本中的新错误,我认为问题可能出在以下代码:

ask my-links with [closed? = false][
  set num-infected-neighbors count (link-neighbors with [infected? = true])
]

这要求一些链接运行程序link-neighbors,但该程序仅设计为由海龟运行。(在 NetLogo 中,某些程序仅被设计为“由代理”运行,例如在 . 之后的大括号内ask ...。然而,其中一些程序仅适用于海龟代理,而其他程序仅适用于链接代理。)您可能应该使用other-end,这是一个设计为由链接运行的过程。(不过,这是一个有趣的过程,因为它仅被设计为在ask请求链接运行的命令other-endask乌龟调用时由链接运行。所以other-end只有当你有类似的东西时才有效ask turtle 1 [ask mylinks [ ... other-end ...]]。)

(可以在第一个代理编辑的另一个代理中引用代理ask。对于其中一个级别,您可以使用myself。对于asks 的其他级别,您可以在外部级别定义变量,然后在内部级别中引用它们。 )

于 2013-08-22T04:35:56.607 回答