Диссертация (1138200), страница 27

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 27)

число исходящих реберfor (i in 1:len){just$fr[i] <- a$lengths[which(a$values == just$id[i])]# вершины, к которым идут ребра из вершиныjust$id[i]temp <- lcut[lcut$from == just$id[i], "to"]if(length(temp) == 0) { just$justfr[i] <- NA} # если нет исходящих реберjust$justfr[i] <- just$justfr[i] + sum(just[just$id %in% temp, "just"])rm(temp)}rm(i,len)just$justfr_a <- (just$justfr / just$fr) # средняя лояльность окружения к коррупции# write.table(just, file = "just_nn.csv", sep = ";", col.names = T, qmethod = "double",row.names=FALSE)library(xlsx)write.xlsx(just, file = "just_nn.xlsx", row.names=FALSE)# корреляцииrequire(Hmisc)# лояльности к коррупции и лояльности к коррупции по окружениюrcorr(just$just, just$justfr_a, type="spearman")# cor.test(just$just, just$justfr_a, method = "spearman")# лояльности к корруцпии и числа мощности узлаrcorr(just$just, just$degree, type="spearman")167# отклонение от среднего по соседям:just$dev <- just$just-just$justfr_arequire(Publish)ci.mean(just$dev, na.rm=T, normal=FALSE)ci.mean(just[just$just > .5, "dev"], na.rm=T, normal=FALSE) #для лояльных к коррупцииci.mean(just[just$just <= .5, "dev"], na.rm=T, normal=FALSE) #для нелояльных к коррупции# матрица переходовmat <- matrix(as_adjacency_matrix(nn), nrow = v, ncol = v)colnames(mat) <- c(nodescut$id)rownames(mat) <- c(nodescut$id)# для justnotna <- c(just$id)theta <- just$justA <- subset(mat, colnames(mat) %in% notna, rownames(mat) %in% notna)Aout <- A %*% thetaAin <- t(A) %*% theta# минимизация нормы эпсилон:norm_eps <- function(alpha) sqrt(sum((theta-alpha*Aout-(1-alpha)*Ain)^2))# t <- nlm(norm_eps, c(0:1))solution <- optimize(norm_eps, lower = 0, upper = 1, maximum = FALSE)alpha <- solution$minimumeps <- theta-alpha*Aout-(1-alpha)*Ainsave.image(file = "december11.Rdata")# neigh$q13fr <- sum(neigh[which(neigh$from)==neigh$from[]]), q13to])#nodes$q13fr <- as.integer(NA)len <- length(nodes$id)for (i in 1:len){for (j in 1:len){if (are_adjacent(nn, i, j)){ #nodes[is.na(nodes$q13fr[i])] <- 0nodes$q13fr[i] <- nodes$q13fr[i]+nodes$q13[j]}}}.

Характеристики

Список файлов диссертации