Dataset_Doron_Shin_et_al_PRhL1_Science_2020/Text_exports/00_exportcode.r

newindex =data.frame()
index <- read.csv("~/github_repo/07_larkumlabrepos/Doron_Shin_PRhL1_paper_data/Text_exports/00_matfiles_index.csv")

g= index$matfile_provenance

for (j in c(1:length(g))){
  data <- readMat(g[j])
  nd=names(data)
  for (i in c(1: length(nd))){
    varname=nd[i]
    filename = paste0(j,"_",i,"_",varname,".tsv")
    if ("matrix" %in% class(data[[i]]) ){
      write.table(data[[i]], paste0("Text_exports/",filename), row.names = FALSE, col.names = FALSE, sep="\t")
      newindex = rbind(newindex, c(filename,g[j],varname, ""))
    } else {
      write.table(as.character(list2DF(data[[i]])), paste0("Text_exports/",filename), row.names = FALSE, col.names = FALSE, sep="\t")
      newindex = rbind(newindex, c(filename,g[j],varname, "no direct export possible"))
    }
    
  }

}
    

names(newindex) = names (index)
write.table(newindex, paste0("Text_exports/","01_file_index.tsv"), row.names = FALSE, col.names = FALSE, sep="\t")

data <- readMat("Figure 3/Fig3G.mat")
nd=names(data)
for (i in c(1: length(nd))){
  varname=nd[i]
  filename = paste0(j,"_",i,"_",varname,".tsv")
  data2= list2DF(data[[i]])
  write.table(t(data2), paste0("Text_exports/",filename), row.names = FALSE, col.names = FALSE, sep="\t")
}


# ## code used to join exports from matlab, files deleted
# a=list.files("Text_exports/fig3_s3_data")
# 
# setwd("Text_exports/fig3_s3_data")
# dataexp= data.frame()
# for (i in c(1: length(a))){
# #  i=1
#   datam=as.data.frame(t(read.csv(a[i], header = FALSE)))
#   colnames(datam)= as.character(c(1:150))
#   print(c(ncol(datam), nrow(datam)))
#   name = stringr::str_split(a[i],"_")
#   datam$mouse=as.character(name[[1]][1])
#   datam$timingcat= as.character(stringr::str_split(name[[1]][2],".csv")[[1]][1])
#   dataexp=rbind(dataexp,datam)
# }
# 
# write.csv(dataexp, "Fig3_s3_alldata.csv")