Corrigindo... Tinha um erro no cálculo do nlink, que retornava sempre 4... Segue a modificação. # <BEGIN> #install.packages('repmis') library(repmis) FinURL <- paste0(" https://dl.dropboxusercontent.com/u/61883020/2013-08-25-r-br.csv") data <- repmis::source_data(FinURL, sep = ",", header = TRUE) dim(data); head(data) data$new = 0; data$prop = 0; data$medlnk= 0; data$nlink=0 str(data) for (i in 1:length(data[,1])){ # para cada linha col = data[i,2]+2 # coluna do grupo dado dado = data[i,col] # dado para a coluna calculada nlink = length(which(data[i,3:6]>0)) # núm de grupos com valor (ligações) data[i,8] = dado # aloja o dado data[i,9] = dado/(data[i,7]-dado) # total - dado do grupo data[i,10] = data[i,7]/nlink # total/núm ligações data[i,11] = nlink } tail(data) #<END>

tail(data) pessoa pessoa.grupo g.1 g.2 g.3 g.4 tc new prop medlnk nlink 15 q 3 8 49 157 17 231 157 2.1216216 57.75000 4 16 r 2 27 138 33 27 225 138 1.5862069 56.25000 4 17 s 2 32 82 25 77 216 82 0.6119403 54.00000 4 18 t 3 12 16 183 3 214 183 5.9032258 53.50000 4 19 u 4 0 54 33 90 212 90 0.7377049 70.66667 3 20 v 3 0 0 135 11 211 135 1.7763158 105.50000 2

-- Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]