1

Quiero obtener z score agrupando por groups. Este es mi código:

data<-read.csv("C:/Users/usuario/Desktop/Tandas/TANDA JAPON/metabolomics/data/fwdreandres/R_PCA_tapR_sini.csv",sep=";",row.names = 1)

head(data) 

La foto adjunta es una representacion de mis datos pues los mismos son demasiado largos como para postear

introducir la descripción de la imagen aquí

data<-scale(data,center = TRUE, scale = TRUE)

Este código calcula z score pero sin tener en cuenta los grupos. Me podeis ayudar?

1
  • 1
    Hola. Este foro sólo usa español; por favor, traduce tu pregunta.
    – Alfabravo
    el 19 nov. 2019 a las 15:00

1 respuesta 1

0

Una forma bien simple en este caso es apoyarse en las rutinas de dplyr, algo así:

library(dplyr)

data %>% 
    group_by(group) %>% 
    mutate(var1 = scale(var1),
           var2 = scale(var2),
           var3 = scale(var3),
           var4 = scale(var4))

O bien

data %>% 
    group_by(group) %>% 
    mutate_at(vars(matches('var')), list(scaled = scale))

Básicamente:

  • Agrupamos por la variable group mediante el verbo group_by
  • Por cada grupo, calculamos una nueva variable con el mismo nombre de la original aplicando scale() al grupo
  • La versión con mutate_at() te aplica scale a todas las columnas con el nombre var* y crea una nueva versión con el nombre *_scaled
  • Recuerda que puedes instalar dplyr con install.packages("dplyr") o install.packages("tidyverse") que agrega varios paquetes muy útiles
1
  • Muchas gracias por la ayuda!!
    – Andrew27
    el 2 mar. 2020 a las 11:44

Tu Respuesta

Al pulsar en “Publica tu respuesta”, muestras tu consentimiento a nuestros términos de servicio, política de privacidad y política de cookies

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.