2016-10-03 4 views
0

Я хотел бы создать график строки x-y для 4 пациентов (PatientID), которые были с течением времени (Age_0 - Age_10) для уровня CRP (CRP_0 - CRP_10).Графический график X-Y с разными значениями X для каждого значения Y

PatientID Age_0 Age_1 Age_2 Age_3 Age_4 Age_5 Age_6 Age_7 Age_8 Age_9 Age_10 CRP_0 CRP_1 CRP_2 CRP_3 CRP_4 CRP_5 CRP_6 CRP_7 CRP_8 CRP_9 CRP_10 
1 22 24 24 28 30 31 0 0 0 0 0 3 9 2 1 0 1 0 0 0 0 0 
3 27 28 29 32 33 35 0 0 0 0 0 2 10 2 1.2 0.1 0 0 0 0 0 0 
4 37 38 39 40 42 43 44 45 0 0 0 8 0 7 7 0 0 7 2 0 0 0 
5 33 35 36 38 39 40 41 0 0 0 0 2 2 5 0.2 0 0 0 0 0 0 0 

structure(list(PatientID = c(1L, 3L, 4L, 5L), Age_0 = c(22L, 
27L, 37L, 33L), Age_1 = c(24L, 28L, 38L, 35L), Age_2 = c(24L, 
29L, 39L, 36L), Age_3 = c(28L, 32L, 40L, 38L), Age_4 = c(30L, 
33L, 42L, 39L), Age_5 = c(31L, 35L, 43L, 40L), Age_6 = c(0L, 
0L, 44L, 41L), Age_7 = c(0L, 0L, 45L, 0L), Age_8 = c(0L, 0L, 
0L, 0L), Age_9 = c(0L, 0L, 0L, 0L), Age_10 = c(0L, 0L, 0L, 0L 
), CRP_0 = c(3L, 2L, 8L, 2L), CRP_1 = c(9L, 10L, 0L, 2L), CRP_2 = c(2L, 
2L, 7L, 5L), CRP_3 = c(1, 1.2, 7, 0.2), CRP_4 = c(0, 0.1, 0, 
0), CRP_5 = c(1L, 0L, 0L, 0L), CRP_6 = c(0L, 0L, 7L, 0L), CRP_7 = c(0L, 
0L, 2L, 0L), CRP_8 = c(0L, 0L, 0L, 0L), CRP_9 = c(0L, 0L, 0L, 
0L), CRP_10 = c(0L, 0L, 0L, 0L)), .Names = c("PatientID", "Age_0", 
"Age_1", "Age_2", "Age_3", "Age_4", "Age_5", "Age_6", "Age_7", 
"Age_8", "Age_9", "Age_10", "CRP_0", "CRP_1", "CRP_2", "CRP_3", 
"CRP_4", "CRP_5", "CRP_6", "CRP_7", "CRP_8", "CRP_9", "CRP_10" 
), class = "data.frame", row.names = c(NA, -4L)) 

Итак, я хотел бы построить график, в котором X есть возраст, а Y - это CRP. Например, Age_0 vs CRP_0 - это одна точка данных для каждого идентификатора пациента. Затем, вторая точка данных Age_1 против CRP_1 и т.д. В конце концов, я хотел бы иметь сюжет что-то вроде этого:

enter image description here

Я был бы рад, если вы можете помочь мне в этом , Большое спасибо.

+1

вообще вопрос лучше воспринимается, если вы объясните, что вы пробовали и где проблема, что вам не могу решить ;-) – agenis

+0

Я забыл упомянуть, что значения «0» фактически отсутствуют или значения NA, которые мы можем вызывать. – Lothlorian

+1

Я просто построил первый график значений (данные $ CRP_0, data $ Age_0, type = "l") ', а затем перечеркнул другие значения с помощью строк:' lines (data $ CRP_1, data $ Age_1) ' – Choubi

ответ

2

Вы можете сделать так (хотя я чувствую, что это, вероятно, не самый эффективный способ):

library(ggplot2); library(reshape2): library(dplyr) 
d1=df %>% select(matches("Age|PatientID")) %>% melt(id.vars=1, value.name="Age") %>% select(-variable) 
d2=df %>% select(matches("CRP|PatientID")) %>% melt(id.vars=1, value.name="CRP") %>% select(-variable, -PatientID) 
toplot = cbind(d1, d2) %>% arrange(PatientID, Age) %>% mutate(PatientID=factor(PatientID)) 
ggplot(toplot[toplot$Age!=0,]) + aes(x=Age, y=CRP, color=PatientID) + geom_point(pch=15, size=4) + geom_line(size=1.5) + coord_cartesian(xlim=c(0, 50)) + xlab("Time") 

Здесь я предположил, что точки с Age=0 отсутствуют значения, но точки с ненулевым возраста и CRP=0 действительны.

enter image description here

+0

Спасибо большое. Именно то, что я искал. Я проверю это сейчас. Приветствия. – Lothlorian

+0

@ Бахти хорошо! Если какой-либо из ответов решает ваш pb, вы можете отметить его как принятый, чтобы другие люди видели, что он разрешен. Если не сообщите нам – agenis

+0

Я оптимизировал ваши команды и Sotos 'один ниже и работал красиво. Большое спасибо за Вашу помощь. Я действительно оценил это. – Lothlorian

2

Я не могу увидеть изображение вы предоставляете, но вот идея,

df1 <- data.frame(PatientID = df$PatientID, age = stack(df[,grep('Age', names(df))])[,1], 
        CRP = stack(df[,grep('CRP', names(df))])[,1], stringsAsFactors = FALSE) 

head(df1) 
# PatientID age CRP 
#1   1 22 3 
#2   3 27 2 
#3   4 37 8 
#4   5 33 2 
#5   1 24 9 
#6   3 28 10 

#converting PatientID to factor, 
df1$PatientID <- as.factor(df1$PatientID) 

#Plot, 

library(ggplot2) 
ggplot(df1, aes(x = age, y = CRP, color = PatientID))+geom_line() 
+0

Это кажется хорошей альтернативой для примера agenis. Я проверю оба, и вы снова узнаете. Большое спасибо Сотос. Я ценю вашу помощь. Лучшие, Бахти – Lothlorian

Смежные вопросы