Блог

Главная - Вопросы по программированию - Есть ли способ получить значения градиента из моего графика регрессии в ggplot2?

Есть ли способ получить значения градиента из моего графика регрессии в ggplot2?

#r #ggplot2 #regression

#r #ggplot2 #регрессия

Вопрос:

Я строю свои результаты регрессии, используя ggplot2 в R. У меня есть данные, разделенные на мужчин и женщин, есть ли способ получить коэффициенты корреляции линий тренда для этих двух групп отдельно, т.Е. Один для женщин и один для мужчин?

Я использую текущий сценарий ниже:

 model <- lm(BMD ~ Sex   Age   PercentDepth   
              Sex*Age   Sex*PercentDepth   Age*PercentDepth, 
            data = SD_MFC_DATA_LENGTHv6)
summary(model)

library(ggplot2)
ggplot(SD_MFC_DATA_LENGTHv6, 
       aes(x = Age,
           y = BMD,
           group = Sex,
           color = Sex,
           shape = Sex))   
  geom_point()  
  geom_smooth(method = lm, aes(fill = Sex)))

Пример данных из SD_MFC_DATA_LENGTHv6:

Пример данных pt 1Пример данных, часть 2

Комментарии:

1. Привет, Сэмми! В stackoverflow вам нужно поделиться воспроизводимым примером. Пожалуйста, добавьте вывод dput(SD_MFC_DATA_LENGTHv5) . Кроме того, ваша модель и ваш график не показывают одно и то же. В lm модели ваша цель BMD , а на графике Tb.Th . Как линейная модель связана с графиком?

2. Привет, Сэмми. У вас есть три разные переменные-предсказателя плюс взаимодействия, поэтому ваша модель фактически описывает две разные плоскости в трехмерном пространстве, а не две линии тренда в 2D-пространстве. Это потому, что для каждого пола будет своя «линия тренда» в зависимости от возраста и МПК для каждого бесконечно малого значения PercentDepth. Если вы оставите PercentDepth из своей модели, вы можете отобразить то, что просите. Если вы оставите это, то построение 2 линий тренда просто не имеет смысла.

3. Привет, Аллан, извини, я должен был изменить его, чтобы они оба были связаны с BMD. Ах, хорошо, трехмерное пространство имеет смысл. К сожалению, мне нужно включить процентную глубину, поскольку это очень важно для результатов. Я могу попытаться создать воспроизводимую версию своих данных, данных много, поэтому я не смог включить их сюда

4. @AllanCameron Я добавил в некоторые примеры данных, к сожалению, я не смог отформатировать их в вопросе, поэтому они добавлены как ссылка