zaza0_0的博客

原创 【python】可变变量与不可变变量

但是，不可变类型的关键特点是它们的值不能被改变。当您修改这个列表（例如 lens[0] = 0），更改会反映在这个列表对象上，因此 counts 也会看到这个更改，因为它们指向的是同一个对象。可变对象（如列表）的变量名是指向对象的引用。多个变量可以指向同一个可变对象，因此一个变量对对象的修改可以通过其他指向同一对象的变量看到。当您创建一个不可变对象（如整数、字符串）时，Python实际上是将这个值存储在内存中的某个位置，并将变量名指向这个值。不可变对象（如整数、字符串）的变量名可以被认为是指向其值的引用。

原创 【python】数据转换basis function并计算转换之后的维数

【代码】【python】数据转换basis function并计算转换之后的维数。

原创 【Github】将本地仓库同步到github上

许久没有用GitHub了，怎么传仓库都忘记了。在这里记录一下。

原创 【python】理解什么时候函数内部的变量会被改变

原创 【python】使用spyder运行文件

【代码】【python】使用spyder运行文件。

原创 【pytorch】安装只需要cpu的pytorch

【代码】【pytorch】安装只需要cpu的pytorch。

原创 【slurm】常用操作

sinfo: 查看nodes的使用情况scontrol show node xxx: 查看某个node分配了多少cpu,内存

原创 【python】面向对象的编程，除内部reassignment/multiprocessing之外传递给函数的可更改对象在函数内部会被修改

【代码】【python】面向对象的编程，除内部reassignment之外传递给函数的可更改对象在函数内部会被修改。

原创 【Git】两行代码恢复删除的stash文件

【代码】【Git】恢复删除的stash文件。

原创 【R语言】facet并不包含全部level的处理, 不同facet调节ylim，ggarrange合并两个ggplot保留配色，并且使用全部水平legend

ggplot合并两个facet, facet每个图调刻度

原创 【R语言】ggplot pie plot +facet_wrap + txt

ggplot pie chart

原创 【Slurm-GCP】暂停正在运行或排队的任务

slurm-gcp暂停正在运行或排队的任务

原创 【slurm-gcp】使用putty/winscp进行文件传输

slurm-gcp ssh远程连接

原创 【GCP-slurm】远程文件传输WinSCP

slurm-gcp远程文件传输

原创 【python】cannot import name ‘asarray‘ from ‘numpy.core.multiarray‘

cannot import name 'asarray' from 'numpy.core.multiarray'

原创 【Eigen库】检查是否某一值在向量或矩阵中

# 检查是否至少有一个，注意为逻辑array，但凡有一个以上sum都是1(I.array() != i).sum(); # 得到有多少个(I.array() != i).cast <int> ().sum()

原创 【python】numpy多维数组的理解

数组引用一篇比较好的文章 https://www.cnblogs.com/yqxg/p/10582215.html。引用首先是最外层的括号，然后向内部。比如3D数组，引用的顺序是层-行-列。print(E4)[[[1] [1]] [[3] [2]]]print(E4.shape)(2, 2, 1)比如E4是一个3D的矩阵。从最外层看，有2层；然后进入一下个括号，有两行；最后进入最后一个括号，有一列；所以维度是(2,2,1)或者也可以从内向外数维度，最内层括号1个元素，第二层括

原创 【R语言】对模拟结果的3维array进行latex表格制作

由于有的表格比较复杂，直接用xtable不好输出，可以先生成csv文件，再用latextablegenerator之类的在线latex表格生成器生成。n_list = c(30, 90)method_list = c(1, 3, 7, 2, 4, 8, 5, 6)metric_list = c(1, 3, 4:(4+p.x-1), 10)weight=5m = 20gen.table <- function(n_list, method_list, metric_list){ tab

原创 【R语言】ggplot对3维array模拟结果数据绘图

library(ggplot2)n_list = c(30, 60, 90)method_list = c(1, 3, 7, 2, 4, 8, 5, 6)metric_list = c(1, 3, 4:(4+2*p.x))weight=5m = 20for(n in n_list){ file_name = paste("results/p.x", p.x, "p.z",p.z, "n",n,"m",m, "sigma", sigma,"nrep",nrep, "weight", weig

原创 【latex】subfigure分多页

\begin{figure}[htbp]\ContinuedFloat*\centering\subfigure[$1a$]{\label{fig:1a}\includegraphics[width=0.45\textwidth]{fig/t_1a.jpg}}\subfigure[$1b$]{\label{fig:1b}\includegraphics[width=0.45\textwidth]{fig/t2_1b.jpg}}\subfigure[$1c$]{\label{fig:1c

原创 【R语言】处理3维数组(array)进行折线图绘制的函数

lineplot.process=function(n_list, method_list, metric_list){ plot_meandat = NULL plot_sddat =NULL for(n in n_list) { file_name = paste("results/p.x", p.x, "p.z",p.z, "n",n,"m",m, "sigma", sigma,"nrep",nrep, "weight", weight,"one",onelevel, sep

原创 【R语言】Error in eval(.FLXgetGroupingVar(formula), data) : 数值‘envir‘参数的长度不是一；Error in eval(.FLXgetGr

把矩阵变为Data.frame形式就行了。

原创 【R语言】正则表达式提取特定字符串

> tmp = strsplit(as.character(colname[i, ]),'[0-9]') # 先拆分出想要的子字符串> tmp[[1]][1] "" "" [3] " - " "" [5] "

原创 【github】Failed to connect to github.com port 443: Timed out

实在无解，无奈以下方法貌似每次都能成功…先设置$ git config --global http.proxy 127.0.0.1:1080再取消…$ git config --global --unset http.proxy

原创 【github】fatal: unable to access ‘https://github.com/.......‘: OpenSSL SSL_read: Connection w

git config --global http.sslVerify “false”之后再次git push即可

原创 【ggplot】合并两个ggplot图片

library(ggpubr)p=ggarrange(p1, p2, labels = c("A", "B"), ncol = 1, common.legend = TRUE, legend = "none")p

原创 【R语言】创建设计阵model.matrix

x = model.matrix(~.,as.data.frame(x))[,-1]

原创 【R语言】查找重复的长度

group_size(group_by(df,a,b))

原创 【R语言】代码美化，等号=变为箭头＜-, 对齐空格

formatR::tidy_app()

原创 【R语言】产生R包的pdf说明文档

原创 【R语言】一键读取合并同一目录下不同子目录中所有的Excel文件

参考自：https://stackoverflow.com/questions/32888757/how-can-i-read-multiple-excel-files-into-r数据在一个大的目录下有很多子目录，每个子目录里有很多excel文件。R语言读取setwd('C:/Users/test/Desktop/20200331 data')file.list <- list.files(pattern='*.xls', recursive = T)# recursive=T读取所有

原创 【R语言】Error in newdata[[nameZ[i]]] : subscript out of bounds

newdata需要是data.frame

原创 【R语言】设置显示语言

Sys.setenv(LANG = “en”)

原创 【ggplot2】改坐标轴breaks标签，加数学表达式，数学表达式legend对齐(align the text of legend)

figure_display_batch <- function(res) { res_tmp <- res size <- dim(res_tmp)[3] res_list <- list() for (i in 1:size) { res_list[[i]] <- res_tmp[, , i] } res_list <- do.call("rbind", res_list) col_name <- c("Prediction

原创 【ggplot2】箱线图geom_boxplot+facet_wrap时每个facet保留离群值，不去除离群值；或者去掉离群值使每个facet的ylim范围不包含离群值

1. 保留离群值将coef设置为5，如果没有fact，直接在 geom_boxplot(coef = 5)。 calc_stat <- function(x) { coef <- 5 n <- sum(!is.na(x)) # calculate quantiles stats <- quantile(x, probs = c(0.1, 0.25, 0.5, 0.75, 0.9)) names(stats) <- c("ymin", "lowe

原创 【R语言】gather函数宽数据变长数据的理解

假设我们有如下数据：希望用education_level作为横坐标，每一行作为一个分组，绘制折线图。这就需要将数据变为长数据格式，产生三列，一列是数据，一列是分组的标签（“Never”, “Once”,…），一列是教育水平，这是最终效果：实现方法：plot_dat = gather(dat, claim_amount, counts, -education_level)claim_amount：原来除去要用来做横坐标的education_level之外的列们统一叫这个 名字；counts：那些

原创 【Rmarkdown】中文，PDF的设置

pdf输出title: “HW 5”author: “Lian”date: “2020/9/3”output:pdf_document:latex_engine: xelatex中文输出title: “Wilcoxon rank sum 与 t test的功效比较”author: “Lian”date: “2020/9/22”output:pdf_document:includes:header-includes: \usepackage{xeCJK}keep_tex: y

原创 【Eigen库】矩阵取列赋值崩掉.leftCols()

加一个取值：.leftCols().eval()

原创 c++11: error: ‘random_device’ is not a member of ‘std’ error: ‘mt19937’ is not a member of ‘std’

Rcpp: error: ‘random_device’ is not a member of ‘std’ error: ‘mt19937’ is not a member of ‘std’需要支持c++11需要支持c++11需要在r里面执行下面这句：Sys.setenv(“PKG_CXXFLAGS”="-std=c++11")

原创 RcppEigen比较vector和matrix

前一篇介绍了查找某元素是否在vector,matrix中https://blog.csdn.net/zaza0_0/article/details/104581400这一篇写逐元素比较。//[[Rcpp::depends(RcppEigen)]]#include<RcppEigen.h>#include<Rcpp.h>using namespace std;us...