qq_38308359-优快云博客

原创 R语言做的坐标轴归一化程序

我这里设计了一个基于LOESS拟合的小程序，可以重新定义坐标轴，这样就可以把不同设备采集的曲线数据用相同坐标轴展示，从而可以进行例如光透过物质后光谱的模拟计算。light

2024-08-14 21:03:38 412

advance3Dlensmatrix<-matrix(data=NA, nrow = Nrow, ncol = Ncol, byrow = FALSE, dimnames = NULL)#创造空矩阵。其实，以这些图片的m行n列像素做平均，在图片没有锐利变化的情况下不会损失图片细节，所获较小矩阵也能获得有关图片的准确信息，下面介绍一个m行n列像素合成一个像素的算法。lengthfactor<-4#这里行列均一致，用4*4=16个原像素完成一个新像素。

2023-11-02 16:31:17 152

原创有趣的R语言方阵均值抽样函数

M<-9#set substract intervalMinterval<-M-1N<-226#FINAL MATRIX LENGTHfor(i in 1:N){ for(j in 1:Minterval){ a<--1*(i+j) z1<-z1[,a] z1<-z1[a,] }}这个函数实现的功能是将方阵中的间隔行列移除，以此代码行向量为例，第一次移除2-9行后，第10行自动称为第2行，之后第二次移除就是3-10行...

2023-11-02 16:22:52 137

原创 R语言persp3D参数效果分析

border = "black"，此语句省略后无网格线theta = 110，管水平旋转phi = 20管垂直旋转clim=c(0,3000)确认颜色统一坐标轴colvar=I，I为数据分布矩阵，z为三维曲面值；

2023-11-02 16:22:26 345

原创用for语句实现的R语言调试好的段落的连续运行

对于调试好的代码，用for(i in 1){}将其囊括，可以节省单步运行的时间

2022-03-01 20:59:01 373 1

原创 R语言reshape包加载出现环境参数错误的解决方式

程序中突然出现 library(reshape)错误: package or namespace load failed for ‘reshape’ in library.dynam(lib, package, package.lib):没有这个DLL ‘plyr’：是不是没有为此架构安装？此外: Warning message:程辑包‘reshape’是用R版本4.1.2 来建造的。重新安装plyr包后，出现：library(reshape)错误: package or namesp.

2022-02-25 16:24:32 5259 1

原创基于R语言的低通滤波算法

这是一个非常经典的巴特沃斯滤波器代码，就我所知，目前使用频域分析和滤波主要是用MATLAB来完成，从百度上查询或是csdn中还未有用R语言设计巴特沃斯滤波器的详细样例，本例用大量的解释和调试工作明确了r语言巴特沃斯滤波器的设计

2022-01-17 11:17:17 1570

原创峰峦图样例及注释

#代码参考《R语言数据可视化之美-专业图表绘制指南》例4.2.1，主要代码链接参见#https://github.com/EasyChart/Beautiful-Visualization-with-R,#在添加了详细的注释基础上，讨论了mydata数据构建的另一种代码

2022-01-07 09:07:10 533

原创常用英文论文佳句

Ultimately, these developments should lead to improvements in characterization of a and b, with concomitant clinical impact.Herein, we derive a a law for measurement of x based on the c technique in turbid tissues, and we validate the approach.Importan

2021-12-31 14:45:33 319 1

原创 LED阵列PCB图设计流程

LED阵列设计有大量重复操作，因此其科学设计具有特殊性，在此总结经验如下1.通过实际实验，确定LED间隔2.通过原理图，确定串并联情况（红线为走底板线）3.建立不连线原理图SCH，4.建立PCB，确认板尺寸和坐标体系，设定单位为mm，设定走线规则，线径等5 导入SCH6 按整体坐标体系布灯，其中包括：截取一排总长度需要的相应数量的灯，M O整体旋转元件到合适角度，E,G,B,水平对齐排列图标，放置在板上；头灯，尾灯输入坐标确定起始和结束位置，EGB水平对齐，EGD水平均匀排列；7

2021-12-25 12:08:43 2795

原创 altium designer2021版本使用探索

PCB部分：PCB上等距布线的最好办法是设置一段长度为所要设置距离的辅助线，然后连辅助线一起复制，粘贴时定位在辅助线末端，并在完成粘贴后ctrl+x清除辅助线；手动布线被挪到了Route-》interactive routing 中；查询元件位置快捷键E J C可以全选准备删除的线段，统一ctrlx删除电路板外形确认在design board shape中侧边的properties中修改公英制，暂时更新这么多，再有新的，在回复中更新...

2021-12-22 13:36:45 692 6

原创基于fanplot包绘制的一个非常好的fanchart样例及详细注释

这是一个非常经典的fanchart图绘制，就我所知，从百度上查询或是csdn中还未有用R语言绘制fanchart图的详细样例，本例用大量的解释和调试工作明确了r语言fanchar图的绘制

2021-12-11 18:12:53 654

原创三个不同供电范围的LED并联共地接三路电流/电压源的BUG问题

三个不同供电范围的LED并联共地接三路电流/电压源的BUG问题

2021-12-10 17:37:42 631

原创 R语言曲面拟合代码详细分析（1）

#代码参考《R语言数据可视化之美-专业图表绘制指南》例4.2.1，主要代码链接参见#https://github.com/EasyChart/Beautiful-Visualization-with-R,#在添加了详细的注释基础上，对一处BUG代码进行了勘误

2021-11-21 14:35:48 1183

原创将各列变量进行分类的代码

mydataln1<-data.frame(mydata$x,mydata$lnmodel1)mydataln1$class<-"lnmodel1"names(mydataln1)<-c("x","y","class")mydataln2<-data.frame(mydata$x,mydata$lnmodel2)mydataln2$class<-"lnmodel2"names(mydataln2)<-c("x","y","class")mydatasum&

2021-11-20 11:55:36 509

原创小提琴图等需要注意的分类排序问题

当画好小提琴图后，需要用scale_x_discrete进行排序，如：p+scale_x_discrete(limits = c( "emodel2", "lmodel2","limodel2","emodel1","lmodel1","limodel1"))才能实现有序排布

2021-11-19 15:22:56 790

原创 COOK距离（influencePlot图的标签设定问题）

influencePlot图设定标签要注意配合row.name命令设定行名，才能和标签一一对应，否则就仅是对应点图mydata3<-mydatamydata3rowname<-round(mydata3$x,3)row.names(mydata3)<-mydata3rownamepar(mfrow = c(2,3))#画图设置influencePlot(lm(y1~0+x,data=mydata3),main="linemodel1")influencePlot(lm(y1~

2021-11-19 14:47:48 1280

原创《机器学习理论和实践》ACF图勘误（P189）

189页代码为plot(acf(linemodel$residuals))实际上acf(linemodel$residuals)本身就可以绘制ACF图，如果使用plot(acf(linemodel$residuals))，在组图显示的时候，会同时绘制2张一样的ACF图。

2021-11-18 14:49:46 506

原创 ggplot与geom中设定aes,图例与facet_zoom

p<-ggplot(mydata)+ geom_line(aes(x=d,y=log10Φ1,color="model1"),size=0.9)+ geom_line(aes(x=d,y=log10Φ2,color="model2"),size=0.9)+ geom_line(aes(x=d1,y=log10Φ1real,color="model1"),size=3,alpha=0.5)+ geom_line(aes(x=d1,y=log10Φ2real,color="model2")...

2021-11-17 19:37:45 644

原创修改文本上下标指令

labs(y=bquote(lnΦ[modelx](d)))+ labs(x=bquote(lnΦ^'modelx'[(d)]))

2021-11-17 10:35:55 148

原创有用的R语言初始设置

rm(list = ls())#清空工作区，可在括号内加装变量等定向清除某个变量等file_path<-file.choose()#在此标注加载文件目录，如E:\工作library(readxl)mydata <- read_excel(file_path)par(mfrow = c(1,1))#画图设置path<-readline()#括号不写内容，默认读取下一行E:\工作\#设定工作目录，输出图片等需要到的目录Mypath<- gsub("\\\\", "/",p

2021-11-16 10:33:26 868

原创 R语言和excel联用时的对数函数问题

R语言中log（x）函数求的是自然对数ln；而excel则是ln()函数，excel的log()函数求的是以10为底的对数，在两种软件联用时易出现计算错误。R语言中用log(x,n)来求以n为底的对数

2021-11-16 10:24:41 1761

原创线性、对数指数拟合模型的建立和验证

#这段程序集成了模型的建立，预测和R2值的计算，发现了一个有趣的问题，R语言输出的参数在EXCEL里进行计算会出现相当大的偏差，因此无法将建立的模型通过excel重新计算输出，全部都应该在R语言内实现rm(list = ls())library(ggplot2)file_path<-file.choose()library(readxl)mydata <- read_excel(file_path)#E:\训练集linemodel<-lm(y1~x,data=mydata)

2021-11-14 10:57:41 1022

原创 r语言对模型输出预测值与实际值计算R2，MS残，RMSE和NRMSD，实现模型评估

#此函数可由模型输出的预测值和实际值计算r2或均方根误差RMSE或归一化后的均方根偏差NRMSD#公式参考吴今朝译《R语言机器学习》323-324页r2.test<-function(y_actual,y_predicted){ avr_y_actual <- mean(y_actual) ss_total <- sum((y_actual - avr_y_actual)^2) ss_regression <- sum((y_predicted - avr_y_act...

2021-11-14 09:32:36 7122 4

原创加载see包后画云雨图特别容易了

#ggplot(mydata, aes(x = num, #y = Age, fill="red")) + #geom_violindot() + #theme_modern()+ #coord_flip()

2021-10-29 21:48:02 157

原创 R语言坐标轴刻度向内指令

theme(axis.ticks.length=unit(-0.25, "cm"), axis.text.x = element_text(margin=unit(c(0.5,0.5,0.5,0.5), "cm")), axis.text.y = element_text(margin=unit(c(0.5,0.5,0.5,0.5), "cm")) )参考ggplot2拟合曲线图（从0开始，刻度向内） - 简书...

2021-10-24 15:13:30 3111

原创 R语言绘图心得2021年10月22日：线性拟合、指数拟合、对数拟合、乘幂拟合及有用的小程序块

rm(list = ls())library(ggplot2)file_path<-file.choose()mydata<-read.csv(file_path,header = T,stringsAsFactors=FALSE)#mydata<-read.csv("figure2b1.csv",stringsAsFactors=FALSE)这个命令是手动拾取文件，能够机动灵活，将文件放在自己指定的文件夹内，下面通过注释可以指出文件位置。linemodel<-..

2021-10-22 16:10:34 2681 3

qq_38308359的博客