는

Question

내가이 예를 들어, 데이터 프레임

> head(df, 10) 
    individual events 
1   1 72.0 
2   1 75.5 
3   1 87.5 
4   2 3.0 
5   2 14.5 
6   2 16.5 
7   2 32.0 
8   2 45.5 
9   2 50.0 
10   2 70.5 

은 내가 하나를 얻을 수 있도록 이벤트의 누적 단계 그래프를 플롯 싶습니다 제공

set.seed(12345) 
n1 <- 3 
n2 <- 10 
n3 <- 60 

times <- seq(0, 100, 0.5) 

individual <- c(rep(1, n1), 
       rep(2, n2), 
       rep(3, n3)) 

events <- c(sort(sample(times, n1)), 
      sort(sample(times, n2)), 
      sort(sample(times, n3))) 

df <- data.frame(individual = individual, events = events) 

을 말해봐 개인 당 라인은 이벤트가 "발생"될 때마다 1 씩 올라갑니다.

예를 들어, 개인 1은 72.0으로 0이 될 것이고 그래프의 끝까지 87.5로 3이 될 때까지 75.5가 될 때까지 최대 1이 될 것입니다.

가장 쉬운 방법은 무엇입니까?

Answer 1

df$step <- 1 

library(plyr) 
df <- ddply(df,.(individual),transform,step=cumsum(step)) 

plot(step~events,data=df[df$individual==1,],type="s",xlim=c(0,max(df$events)),ylim=c(0,max(df$step)),xlab="time",ylab="step") 
lines(step~events,data=df[df$individual==2,],type="s",col=2) 
lines(step~events,data=df[df$individual==3,],type="s",col=3) 

Answer 2

사용 ggplot2는 :

library(ggplot2) 

# Add step height information with sequence and rle 
df$step <- sequence(rle(df$individual)$lengths) 

# plot 
df$individual <- factor(df$individual) 
ggplot(df, aes(x=events, group=individual, colour=individual, y=step)) + 
    geom_step() 

Answer 3

또한 통계 패키지의 stepfun 기능이있다. 그 사용하면 해당 개체 클래스의 plot 방법을 사용할 수 있습니다

는

sdf <- split(df, individual) 

plot(1, 1, type = "n", xlim = c(0, max(events)), ylim = c(0, max(table(individual))), 
    ylab = "step", xlab = "time") 

sfun <- lapply(sdf, function(x){ 
    sf <- stepfun(sort(x$events), seq_len(nrow(x) + 1) - 1) 
    plot(sf, add = TRUE, col = unique(x$individual), do.points = FALSE) 
})