Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- object1 <- rnorm(500,mean=1,sd=0.1)
- object2 <- rnorm(500,mean=2,sd=0.1)
- object3 <- rnorm(500,mean=3,sd=0.1)
- estimations <- data.frame(object1,object2,object3)
- colnames(estimations) <- 1:3
- SEframe <- data.frame()
- SEframe <- rbind(SEframe,c(1,1,mean(object1),0.1))
- SEframe <- rbind(SEframe,c(2,2,mean(object2),0.15))
- SEframe <- rbind(SEframe,c(3,3,mean(object3),0.25))
- colnames(SEframe) <- c("ID","True.value","Estimated.value","SE")
- estMelted <- melt(estimations)
- estMelted$variable <- as.factor(estMelted$variable)
- p <- ggplot(estMelted, aes(x=variable, y=value)) +
- geom_violin()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement