[R-br] Pre-tratamento de dados
Benilton Carvalho
beniltoncarvalho em gmail.com
Sexta Junho 22 11:16:50 BRT 2012
E partindo da resposta do Tura, mais simplificacoes podem vir.
Exemplos abaixo, nao testados, a partir da matriz B1 proposta pelo
Tura:
set.seed(1)
B1 <- matrix(rpois(1209,10),nrow=39,ncol=31)
library(matrixStats)
N1 <- scale(B1, colQuantiles(B1, 0), rowDiffs(colQuantiles(B1, c(0, 1))))
N2 <- scale(B1, colMeans(B1), colQuantiles(B1, 1)-colMeans(B1))
N3 <- scale(B1, colMeans(B1), rowDiffs(colQuantiles(B1, c(0, 1))))
N4 <- scale(B1, colMeans(B1), colVars(B1))
b
> Tales sua resposta é a seguinte:
>
> set.seed(1)
> B1 <- matrix(rpois(1209,10),nrow=39,ncol=31)
> colunas <- ncol(B1)
> linhas <-nrow(B1)
>
> Mmin <- matrix(rep(apply(B1,2,min),linhas),ncol=colunas,byrow=T)
> Mmax <- matrix(rep(apply(B1,2,max),linhas),ncol=colunas,byrow=T)
> Mmean <- matrix(rep(apply(B1,2,mean),linhas),ncol=colunas,byrow=T)
> Mvar <- matrix(rep(apply(B1,2,var),linhas),ncol=colunas,byrow=T)
>
> N1<-(B1-Mmin)/(Mmax-Mmin)
> N2<-(B1-Mmean)/(Mmax-Mmean)
> N3<-(B1-Mmean)/(Mmax-Mmin)
> N4<-(B1-Mmean)/Mvar
> N5<-sqrt(B1)
> N6<-log(B1)
> N7<-log(B1,10)
> N8<-1/sqrt(B1)
> N9<-1/B1
> N10<-log((B1-Mmin)/(Mmax-B1))
>
> Observe que:
> 1- as duas primeiras linhas foram feitas para eu conseguir um exemplo
> minimamente reprodutível da matriz B1, que você deveria ter fornecido.
>
> 2- se tivesse descrito as matrizes - algo como Mmin= menor valor de cada
> coluna de B1 - me teria facilitado entender o código
>
> 3- A familia de comandos *pply resolve muita coisa vale a pena conhecê-los
>
> []s
> Tura
> P.S concordo com os outros que você poderia explicar melhor o problema...
>
>
>
>
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
> mínimo reproduzível.
Mais detalhes sobre a lista de discussão R-br