Použití randomizací k vyhodnocení asortativnosti

V tomto cvičení zjistíš, jak pravděpodobná je pozorovaná asortativnost v síti přátelství s ohledem na pohlaví vrcholů – a to pomocí randomizační procedury. Pohlaví vrcholů v síti náhodně permutuj 1000krát a pro každou náhodnou síť znovu vypočítej asortativnost.

Toto cvičení je součástí kurzu

Analýza sítí v R

Zobrazit kurz

Pokyny k cvičení

Pomocí assortativity() vypočítej asortativnost objektu grafu g1 na základě pohlaví s využitím objektu values vypočítaného v předchozím cvičení a výsledek přiřaď do objektu observed.assortativity.
Uvnitř smyčky for vypočítej asortativnost sítě g1 pomocí assortativity() a pokaždé náhodně permutuj objekt values pomocí sample().
Zobraz distribuci hodnot asortativnosti z této permutační procedury pomocí hist() a přidej červenou svislou čáru odpovídající původní pozorované hodnotě asortativnosti sítě g1 uložené v observed.assortativity.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# Calculate the observed assortativity
observed.assortativity <- ___(g1, values)

# Calculate the assortativity of the network randomizing the gender attribute 1000 times
results <- vector('list', 1000)
for(i in 1:1000){
  results[[i]] <- ___(g1, sample(___))
}

# Plot the distribution of assortativity values and add a red vertical line at the original observed value
___(unlist(results))
abline(v = ___, col = "red", lty = 3, lwd=2)

Upravit a spustit kód