Entering edit mode
20 months ago
Confused_human
▴
30
Hello everyone
I am trying to make a cluster dendrogram for my rna seq gene expression data.
This is the script am using
set.seed(786)
file_loc <- 'seeds.txt'
seeds_df <- read.table('data.txt', header =T)
str(seeds_df)
summary(seeds_df)
any(is.na(seeds_df))
seeds_label <- seeds_df$type.of.seed
seeds_df$type.of.seed <- NULL
str(seeds_df)
seeds_df_sc <- as.data.frame(scale(seeds_df))
summary(seeds_df_sc)
dist_mat <- dist(seeds_df_sc, method = 'euclidean')
hc <- hclust(dist_mat)
plot(hc)
but every time I am having this error:
R aborted or encountered a fatal error
Please help on this troubleshoot
Have you run the script line by line and found which line is causing the error?
you may be producing a super large file and you computer can't handle it