ковариационная матрица и корреляционная матрица

Question

Здравствуйте, я использую данные дистрофии из пакета ipred.Я использовал подмножество для отделения от носителей и нормальное значение:

carrier = subset(dystrophy,dystrophy$Class == "carrier")
normal = subset(dystrophy,dystrophy$Class == "normal")

, и я сократил эти данные, выбирая только пациентов с 1 посещением в больнице:

carrier = subset(carrier,carrier$OBS == "1")
normal = subset(normal,normal$OBS == "1")

ТакТеперь я хотел бы попрактиковаться в расчете вектора средних, ковариационной матрицы и корреляционной матрицы белков, но по отдельным группам (фактор класса).

Я пробовал с помощью cor и cov, но я думаю, что я делаючто-то не так.Любая помощь будет оценена.спасибо !!

Answer 1

Это может помочь вам начать.Используя ваши переменные, вы можете получить средства для каждого из белков, используя:

sapply(carrier[,6:9], mean, na.rm=T)
sapply(normal[,6:9], mean, na.rm=T)

Для корреляции и ковариации вы можете использовать:

cor(carrier[,6:9], use="pairwise.complete.obs")
cor(normal[,6:9], use="pairwise.complete.obs")

cov(carrier[,6:9], use="pairwise.complete.obs")
cov(normal[,6:9], use="pairwise.complete.obs")

Часть 6:9 естьограничить вычисление белками и не включать другие функции, такие как возраст.Часть use="pairwise.complete.obs" предназначена для обработки пропущенных значений.