R语言_找波峰波谷求波长_果蝇心管数据处理.docx
R语言一找波峰波谷求波长一果蝇心管数据处理ZW2024-03-04加载xlsx包&导入Excel数据# 载入R包Iibrary(Xlsx)library(ggplot2)library(ggalt)ggalt,:# #RegisteredS3methodsoverwrittenby# #methodfrom# #grid.draw.absoluteGrobggplot2# #grobHeight.absoluteGrobggplot2# #grobWidth.absoluteGrobggplot2# #grobX.absoluteGrobggplot2# #grobV.absoluteGrobggplot2# 导入心管界面内径的测量数据sheets<-getSheets(IoadWorkbook("E:Desktopheart_tube.×lsx")# ExceL表格名称names(sheets)# #1,Sheetl""Sheet2""Sheet3""Sheet4",Sheet5,'# 每个表格生成变量for(iinl:length(sheets)assign(paste('Dm_heart_dt'i,ijsep=''),read.xlsx("E:WDesktopWheart_tube.×lsx",sheetindex="header=F)# 所有表格合并为一个ListDm_heart_dt<-list(sheetl=Dm_heart_dt_l$X2>sheet2=Dm_heart_dt_2$X2,sheet3=Dm_heart_dt_3$X2>sheet4=Dm_heart_dt_4$X2>sheet5=Dm_heart_dt_5$X2)# 作折线图plot(Dm_heart_dt_l$X2>type=',b")# 美化折线图ggplot(Dm-heart-dt-ljaes(Xl,X2)+geom_xspline(spline_shape=0.5)+labs(×='Time,jy="Length')+theme_test(base_size=20)+theme_classic(base_size=18,base_family=,serif')# #Warning:Usingthe'size'aestheticinthisgeomwasdeprecatedinggplot23.4.0.# #iPleaseuse'linewidth'inthe'default_aes'fieldandelsewhereinstead.# #Thiswarningisdisplayedonceevery8hours.# #Call'lifecycle:last_lifecycle_warnings()'toseewherethiswarningwas# #generated.# 保存图片ggsave("Dm_heart_dt_l.png9width=8,height=5)波峰波谷变量x,波峰=(x-l)<=x&x>=(x+I),波谷=(x-l)>×S×<(x+l)# 波峰peak<-function(x)x_maxs<-c()x_max<-c()x_i<-c()for(iin2:(length(x)-l)if(xi-l<=xi&xi>=×i+l)x_max<-xix_maxs<-c(x_maxs,x_max)x_i<-C(X_i,i)y<-data.frame(peak_x=x_i,peak_y=x_maxs)return(y)peak_test<-peak(Dm_heart_dt_l$X2)# 施客valley<-function(x)x_mins<-c()x_min<-c()×.i<-c()for(iin2:(length(x)-l)if(xi-l>×i&xi<xi+l)x_min<-xix_mins<-c(x_mins,x_min)x_i<-c(x_i,i)y<-data.fname(valley_x=x_i,valley_y=x_mins)return(y)# 找波峰波谷简化：# 波峰peakl<-function(x)x_maxs<-c()x_max<-c()for(iin2:(length(x)-1)if(xi-l<=xi&xi>=xi+l)x_max<-×ix_maxs<-C(X_maxs,x_max)return(x_maxs)peakl_data<-Iapply(Dm_heart_dt,peakl)# 波容valleyl<-function(x)x_mins<-c()x_min<-c()×.i<-c()for(iin2:(length(x)-l)if(xi-l>xi&xi<xi+l)x_min<-xix_mins<-c(x_mins,x_min)x_i<-c(x_i,i)y<-data.fname(valley_x=x_i,valley_y=x_mins)return(y)# 找波峰波谷求波峰持续时间# 分类汇总:求波峰个数,适合波峰值不一致的数据peak_count_a<-aggregate(×=peak_test$peak_y,by=list(peak_test$peak_y),FUN=length)# rLe函数计算向量中连续重复元素的长度:求波峰个数peak_count_r<-rle(peak_test$peak_y)$lengthpeak_counts<-lapply(peakl_data,rle)# 茄诙持续时间计算peak_length<-list()for(iin1:length(peak_counts)peak_lengthi<-peak_countsi$lengths*0.4)# 改list内的名称names(peak_length)<-names(sheets)capture.output(peak_length,file="peak_length.txt")求邻近波峰间距# 求邻近波峰间距wavelength<-function(x)wl<-c()for(iin2:(length(x)-2)if(xi-l<=xi&xi>xi+l)xiapo<-ifor(minxiapo:(length(x)-2)if(xm-l<xm&xm>=xm+l)shangpo<-mbreakwl<-c(wljshangpo-xiapo)return(wl)# 示例wavelength(Dm_heart_dt_3$X2)# #1565656566# 批量wave_length<-list()for(iinl:length(sheets)Dm_h_d<-get(paste("Dm_heart_dt_n,i,sep=",)wave_lengthi<-wavelength(Dm_h_d$X2)*0.4# 改List内的名称names(wave_length)<-names(sheets)capture.output(wave_length,file="wave_length.txt")
