Compare commits

..

2 Commits

+15 -1
View File
@@ -16,7 +16,21 @@ ggheatmap<-function(df, x=NULL, y=NULL, value=NULL, grouping="mean", exclude_gro
if (grouping == "median"){
df<-df %>% summarise(Value=median(Value)) %>% ungroup
}
order<-clustsort(df %>% spread(VarY,Value) %>% select(!all_of(exclude_group)) %>% as.data.frame)
if (length(unique(df$VarX)) > 1 & length(unique(df$VarY)) > 1){
order<-clustsort(df %>% spread(VarY,Value) %>% select(!all_of(exclude_group)) %>% as.data.frame)
}else{
order<-list()
if (length(unique(df$VarX)) > 1){
xhclust<-clustsort(df %>% spread(VarY,Value) %>% select(!all_of(exclude_group)) %>% as.data.frame)
order[["x"]]<-pull(xhclust, 1)[hclust(dist(xhclust %>% select(-1)))$order]
}else{order[["x"]]<-df %>% pull(VarX) %>% unique}
if (length(unique(df$VarY)) > 1){
yhclust<-clustsort(df %>% spread(VarY,Value) %>% select(!all_of(exclude_group)) %>% as.data.frame)
order[["y"]]<-colnames(yhclust)[2:ncol(yhclust)][hclust(dist(t(yhclust %>% select(-1))))$order]
}else{order[["y"]]<-df %>% pull(VarY) %>% unique}
}
if (scale != "none"){
if (scale == "rows"){