# Reading Melissa data to understand what she did # Revised to read "Arena_LaurenJoe.txt" # Revised to read a list of Green Crab .txt files rm(list=ls()) require(rgl) List=dir(patt='[.]txt') for (i in 1:length(List)){ XX=read.csv(List[i], header=FALSE) XX=as.matrix(XX) XX[XX[,1]==999.99,]=NA XXo=XX # print(attributes(XXo)) lst=(1:length(XXo[,1]))[is.na(XXo[,1])] cat(List[i],"NAs at:",lst,'\n') ntps=length(lst) open3d("windowRect"=c(50,50,550,550)) par3d("windowRect"=c(50,50,550,550)) legend3d("topright", legend = paste('Type', c('A', 'B', 'C')), pch = 16, col = rainbow(3), cex=1, inset=c(0.02)) bgplot3d({ plot.new() title(main = paste(List[i],'recto'), line = 3) }) rgl.points(XXo[lst[1]:lst[3],], color=rainbow(lst[3]), size=10) rgl.linestrips(XXo[lst[1]:lst[2],], color='black', lwd=2) rgl.linestrips(XXo[lst[2]:lst[3],], color='red', lwd=2) text3d(XXo[lst[1]:lst[2],],texts=lst[1]:lst[2],adj=1.5, color='black',cex=0.75) text3d(XXo[lst[2]:lst[3],],texts=lst[2]:lst[3],adj=-0.75,color='red',cex=0.75) browser("Type 'c' to exit!", skipCalls=0) rm(XX, XXo) rgl.close() } # end of List