[R-es] Error lectura

[Manuel Máquez]

Los datos originales están Tdist.csv que adjunto, el código es:
mc <- matrix(0,39,48,byrow = T)           # Matriz Cta
mp <- matrix(0,39,48,byrow = T)           # Matriz Pr
ms <- matrix(0,39,48,byrow = T)           # Matriz Sumas
mtzbin <- matrix(0,39,48,byrow = T)       # Matriz Binom
vec <- 0                                  # Valor temporal
dt <- read.csv('Tdist.csv', header = F)   # Datos
lg <- 0                                   # Inicial
lg <- length(dt$V1)                       # Valor longitud
for(i in 1:lg) {                          # Registro
    for(j in 1:39) {                      # Elemento
        da <- dt[i ,j]                    # Dato
        if(da == 0) {next
            }else{
            mc[j, da] <- mc[j, da] + 1
            ms[j, da] <- ms[j, da] + da
            it <- da/mc[j, da]
            mp[j, da] <- it
        }
    }
}
sma <- rowSums(ms, 1)
sr <- max(sma)
write.csv(ms, file = 'dato.csv')
datIn <- read.table('dato.csv', sep = ',', as.is = T)
for(i in 1:39) {
    for(j in 1:48) {
        vec <- datIn[i, j]
    #    mtzbin[i, j] <- ((vec/sma[j])^vec)*((1-
vec/sma[j])^(sma[j]-vec))*choose((sr-1),(vec-1))
        mtzbin[i, j] <-
((vec/sma)^vec)*((1-vec/sma)^(sma-vec))*choose((sr-1),(vec-1))
    }
}
Se trata de leer los datos por medio de la primera parte del código, para
después construir una tabla de probabilidad binomial negativa, vía la
segunda parte del mismo.
Esa primera parte da por resultado 'dato.csv' que también adjunto.
Carlos, nuevamente muchas gracias.
Atentamente;
*MANOLO MÁRQUEZ P.*
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20170908/98915176/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Tdist.csv
Type: text/csv
Size: 60123 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20170908/98915176/attachment-0002.csv>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: dato.csv
Type: text/csv
Size: 4975 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20170908/98915176/attachment-0003.csv>