- 博客(85)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 `n()` must only be used inside dplyr verbs.
data <- data.frame(value = 1:9, # Create example data group = letters[1:3])data # Print example datainstall.packages("dplyr") # Install & load dplyr packagelibrary("dplyr")install.pa
2021-08-14 11:47:35
2988
原创 【R】 heatmap.2()参数
gplots (version 3.1.1)heatmap.2: Enhanced Heat MapThis heatmap provides a number of extensions to the standard R heatmap function.heatmap.2 (x, # dendrogram control Rowv = TRUE, Colv=if(symm)"Rowv" else TRUE,
2021-07-15 15:20:49
3548
原创 R package
library(devtools)library(roxygen2)devtools::build()devtools::document()devtools::load_all()
2021-06-26 17:06:06
397
原创 Linux安装R之后设置环境变量并写入bashrc中
输入vim .bashrc进入编辑模式，按i插入这一句export PATH=/home/…/R/bin:$PATH(…是R的安装# Source global definitionsif [ -f /etc/bashrc ]; then. /etc/bashrcfiexport PATH=/home/…/R/bin:$PATH
2021-06-18 13:05:02
2697
原创 指数分布族的自然参数空间 $A(\eta)$（与充分统计量的关系）
自然参数空间 A(η)A(\eta)A(η)1. 利用A(η)A(\eta)A(η)可以计算矩例如对于二项分布η(p)=logp1−p,A(η)=nlog(1+eη)\eta(p) = \log\frac{p}{1-p},A(\eta)=n\log(1+e^{\eta})η(p)=log1−pp,A(η)=nlog(1+eη)首先有p=eη1+eηp = \frac{e^{\eta}}{1+e^{\eta}}p=1+eηeηA′(η(p))=neη(p)1+eη(p)=npA′′(η(p)
2021-04-13 11:56:42
1673
原创 计算指数分布族的矩（统计推断）
计算指数分布族的矩（统计推断）假定指数分布族的形式如下：f(x∣θ)=h(x)c(θ)exp(∑i=1kwi(θ)ti(x))f(x|\pmb\theta) = h(x)c(\pmb\theta)\exp\left(\sum^k_{i=1}w_i(\pmb\theta)t_i(x) \right)f(x∣θθθ)=h(x)c(θθθ)exp(i=1∑kwi(θθθ)ti(x))要保证h(x)h(x)h(x)和c(θ)c(\pmb\theta)c(θθθ)恒大于等于0则，均值和方差计算方式如
2021-04-09 16:37:45
1520
原创 R语言 solve()函数 求解线性方程组和矩阵的逆 C++
# solve(A,b) 求解线性方程组# solve(A) 求矩阵A的逆V <- matrix(...)B <- c(...)bett <- solve(V, B)J <- solve(V)C++#include <Eigen/Dense>求解线性方程组V.colPivHouseholderQr().solve(B)求逆V.inverse()QR分解速度较慢，可选择cholesky分解V.llt().solve(B)...
2021-02-22 19:37:27
8896
原创 C++ 求正态分布密度
#include <iostream>#include <cmath> //exp() log()#define _USE_MATH_DEFINES //pidouble dnorm(double value){ double fx = (1.0/(sqrt(2*M_PI)))*exp(-pow(value,2)/2); return(fx);}// dnorm(0) = 0.3989423 in RR语言 dnorm()..
2021-02-21 16:14:06
807
原创 C++ 计算正态分布积分
//https://stackoverflow.com/questions/2328258/cumulative-normal-distribution-function-in-c-c#include <cmath> //exp() log()#include <algorithm> //nth_element()#include <iomanip> //erfcdouble normalCFD(double value){
2021-02-20 22:59:34
1413
原创 C++ 求分位数
double quantile(vector<double>&x, double q){ assert(q >= 0.0 && q <= 1.0); const int n = x.size(); double id = (n-1)*q; int lo = floor(id); int hi = ceil(id); double qs = x[lo]; double h = (id-lo); r
2021-02-18 23:16:15
3133
1
原创 Warning in install.packages : package ‘limma’ is not available (for R version 4.0.2)
问题：> install.packages("limma")Warning in install.packages : package ‘limma’ is not available (for R version 4.0.2)解决：To install this package, start R (version “4.0”) and enter:if (!requireNamespace("BiocManager", quietly = TRUE)) install.pac
2021-01-13 13:08:34
4032
原创 【转】python图片切分（一张图片切分为多张小图片）
转自其他人的博客or网页，找不到原出处了qaq# -*- coding: utf-8 -*-"""Spyder EditorThis is a temporary script file."""import osfrom PIL import Imagedef splitimage(src, rownum, colnum, dstpath): img = Image.open(src) w, h = img.size if rownum <= h and
2021-01-11 19:00:22
1981
1
原创 R语言 【rexp】
R文档生成指数分布随机数的函数rexp(n, rate = 1)n表示生成的随机数个数x <- rexp(100,1/10) #生成100个均值为10的指数分布随机数
2020-12-23 13:54:35
11724
原创 R语言 【integrate】函数
R文档Integration Of One-Dimensional Functions一维函数的积分在有限或无限区间上一元函数的自适应求积integrate(f, lower, upper, …, subdivisions = 100L, rel.tol = .Machine$double.eps^0.25, abs.tol = rel.tol, stop.on.error = TRUE, keep.xy = FALSE, aux = NULL)Augme
2020-12-23 10:56:50
10092
4
原创 Soft Thresholding 求解优化问题例子
软阈值(SoftThresholding)可以求解如下优化问题：arg minx∣∣x−B∣∣22+λ∣∣x∣∣1\argmin_x ||x-B||^2_2 + \lambda ||x||_1xargmin∣∣x−B∣∣22+λ∣∣x∣∣1的解为：soft(B,λ/2)={B+λ/2,B<−λ/20,∣B∣<λ/2B−λ/2,B>λ/2soft(B,\lambda / 2)=\left\{\begin{aligned}B+\lambda /2 &,&
2020-12-23 09:19:43
592
原创 未知错误，可能由于拨号连接未创建成功
关于校园网客户端登陆失败未知错误，可能由于拨号连接未创建成功不知道该说是win10垃圾还是浙大校园网客户端垃圾，或者是我太垃圾QAQ重点是：cc98er都是人才存到博客里，以防我的win10又出同样的毛病感恩的心Windows Registry Editor Version 5.00[HKEY_LOCAL_MACHINE\SYSTEM\CurrentControlSet\Services\EventLog\Parameters]“ServiceDLL”=hex(2):25,00,53,
2020-12-16 18:07:31
7109
4
原创 电脑自带的相机软件打开黑屏
对着网上好多办法试了，都不管用最后，发现键盘F8上面有一个禁用相机的图标哦 原来是我不知道什么时候把它禁用了按一下就好了被自己笨死
2020-11-18 19:14:21
3589
6
原创 关于修改Jupyter notebook默认路径
对着网上各种方法试了一遍，都改不过来最后发现是config里路径格式写的不对原来用的是F:\\jupyter notebook改为F:/jupyter notebook成功我也不知道为什么系统是win10
2020-11-07 15:33:47
341
原创 R包安装失败 package ‘mculst’ is not available (for R version 4.0.0)
错误代码：WARNING: Rtools is required to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding:https://cran.rstudio.com/bin/windows/Rtools/Warning in install.packages : package ‘mc
2020-06-20 12:08:48
3611
1
原创 1. R语言入门
文章目录1. R语言入门1.1 一个简单例子1.2 简单会话在这里插入图片描述1.3 函数入门1.4 数据结构1. R语言入门R可以在两种模式下运行：交互模式和批处理模式（类似Python吧）1.1 一个简单例子考虑一个标准正态分布 XXX ，生成一个新的随机变量 Y=∣X∣Y=|X|Y=∣X∣ ，下面来计算 YYY 的均值,我们的方法基于模拟 N(0,1)N(0,1)N(0,1) 分布随机变量的取值：> means(abs(rnorm(100)))# [1] 0.731308这行代
2020-05-18 00:29:21
1237
原创 TF-IDF+sklearn 读取文件夹中文件并给出每个文件主题
# 第一步：import ospath = 'export_output_clean' #文件夹目录fileLists= os.listdir(path) #得到文件夹下的所有文件名称print(fileLists)s = []for file in fileLists: #遍历文件夹 if not os.path.isdir(file): #判断是否是文件夹，不是...
2020-04-22 18:36:54
385
原创 ArcGIS 10.6 教程 第03章【基础篇】坐标系
文章目录坐标系地理坐标系椭球体大地基准面中国常用坐标系投影坐标系高斯克吕格投影GK墨卡托投影UTMGIS操作Gauss_Kruger 变为GCS_Xian_1980GCS_Xian_1980 变为 Gauss_Kruger 20WGS1984 变为 UTM 50N坐标系地理坐标系 Geographical Coordinate System投影坐标系 Projected Coordi...
2020-04-15 14:37:14
2362
1
原创 ArcGIS 10.6 教程 第02章【基础篇】基础操作 2
文章目录三、编辑几何数据（shapefile）四、编辑属性表（字段操作、属性表编辑）五、数据选择（属性选择、空间选择）六、符号化显示七、专题图制作三、编辑几何数据（shapefile）导入点、线、面 .shp文件，外加一份谷歌地图，用来划分区域插入影像图时如果不显示，可以右键选择缩放至图层创建要素修改属性打开属性表可以看到不同属性，不同颜色属性模板...
2020-04-13 20:02:12
1007
1
原创 在Anaconda下安装tensorflow 2.0 （超级简单，不用pip，不用额外的安装cuda、cudnn！）
我的安装历程：安装了tensorflow 2.0 CPUCPU运行了一个CNN之后感觉要炸了新建环境安装 tensorflow 2.0 GPU版本运行了一个程序之后发现我的显卡驱动和cuda版本不匹配查看 cuda是10.0版本，而我的 930m 仅支持9.0更换 cuda 为 9.0 同时它把我的 tensorflow 2.0 也更改为了 1.14， 好吧运行了一个 CNN ，...
2020-04-06 09:17:54
1302
2
原创 import cv2 , ImportError: DLL load failed
Tensorflow 2import cv2error， so I install opencvthen，DLL load failedso I install py-opencvlibopencv 是在 install py-opencv 时自动 install 的成功解决不知道为什么就这样解决了 error，py-opencv 我是随便下的 QAQ另外，需要 i...
2020-04-05 16:13:42
104
原创 module ‘tensorflow._api.v2.train‘ has no attribute ‘AdamOptimizer‘ （不必更改tensorflow 2）
Tensorflow 2.0model.compile(optimizer = tf.train.AdamOptimizer(), loss = 'sparse_categorical_crossentropy', metrics=['accuracy'])出错改为：model.compile(optimizer = tf.opt...
2020-04-04 22:08:54
27061
14
原创 fashion_mnist 数据集加载 附网盘链接
Tensorflow学习链接：https://pan.baidu.com/s/1c4IZITtl3dUbR8p4qMg7lQ提取码：8tn0或（嗯…通常network error，我科学上网了也还是加载不出来）fashion_mnist = keras.datasets.fashion_mnist(train_images, train_labels),(test_images, tes...
2020-04-04 20:41:41
1407
3
原创 Berkeley CS 61A Lecture10-15
http://inst.eecs.berkeley.edu/~cs61a/sp18/文章目录CS 61A LECTURE 10 - CODEData AbstractionPairsReducing to Lowest TermsAbstraction BarriersViolating Abstraction BarriersData RepresentationsRational Data ...
2020-03-19 17:46:00
304
原创 reshape(-1, 1)
cs231n assignment 1这里-1是指未设定行数，程序自动计算，所以这里-1表示任一正整数examplereshape(-1, 1) 表示（任意行，1列），4行4列变为16行1列reshape(1, -1) 表示（1行，任意列），4行4列变为1行16列reshape(-1, 2) 表示（任意行，2列），4行4列变为8行2列reshape(-1, 1) 等多用于不知道具体行数...
2020-03-17 12:49:39
1373
原创 plt.subplot()
cs 231n assignment 1plt.subplot()函数用于直接指定划分方式和位置进行绘图# subplot(numRows, numCols, plotNum)# or subplot(numRows numCols plotNum)图表的整个绘图区域被分成 numRows 行和 numCols 列然后按照从左到右，从上到下的顺序对每个子区域进行编号，左上的子区域的编号...
2020-03-17 10:38:23
4732
原创 np.random.choice()
cs 231n assignment 1# np.random.choice(a, size=None, replace=True, p=None)idxs = np.random.choice(idxs, samples_per_class, replace=False)从 a （ndarray 、 int 、list 或 tuple等均可，需要是一维的） 中，以概率p随机选取大小为si...
2020-03-16 20:45:57
300
原创 np.flatnonzero()
cs 231n assigment 1函数输入一个矩阵，返回扁平化后矩阵中非零元素的位置（index）x = np.arange(-2, 3) # array([-2, -1, 0, 1, 2])np.flatnonzero(x) # Non zero : array([0, 1, 3, 4])import numpy as npd = np.array([1,2,3,4,...
2020-03-16 20:22:13
1426
原创 Berkeley CS 61A Lecture 9
http://inst.eecs.berkeley.edu/~cs61a/sp18/CS 61A LECTURE 9Function AbstractionsChoosing namesWhich values deserve a nameTestingthe greatest common divisor of two integers m and n:#testpyth...
2020-03-12 20:42:10
142
原创 Berkeley CS 61A Lecture 8
http://inst.eecs.berkeley.edu/~cs61a/sp18/CS 61A LECTURE 8The Cascade Functiondef cascade(n): if n < 10: print(n) else: print(n) cascade(n // 10) print(n)>>> cascade(5)5&g...
2020-03-12 14:11:32
195
原创 Berkeley CS 61A Lecture 7
http://inst.eecs.berkeley.edu/~cs61a/sp18/CS 61A LECTURE 7Recursive Function递归函数Definition: A function is called recursive if the body of that function calls itself, either directly or indirectly...
2020-03-11 23:15:00
178
原创 Berkeley CS 61A Lecture 6 （Mario！！！）
http://inst.eecs.berkeley.edu/~cs61a/sp18/CS 61A LECTURE 6Iteration Example: Fibonacci NumbersLecture 4讲过Fibonacci Numbersdef fib(n): """Compute the nth Fibonacci number. >>> fib(0)...
2020-03-11 20:57:50
246
原创 Berkeley CS 61A Lecture 5
http://inst.eecs.berkeley.edu/~cs61a/sp18/CS 61A LECTURE 5Environments Enable Higher-Order FunctionsHigher-order function: A function that takes a function as an argument value or returns a functio...
2020-03-11 17:38:37
160
原创 Python处理csv
Python处理csv与专有的二进制文件格式不同，CSV通常用于在电子表格软件和纯文本之间交互数据。CSV文件内容仅仅是一些用逗号分隔的原始字符串值。专门用于解析和生成CSV的库。将数据转成csv格式写出并再次读取：'''导入csv模块以及distutils.log.warn()后者作为print语句或函数的代理'''import csvfrom distutils.log i...
2020-03-11 15:22:13
149
Andrew Ng tensorflow P1 - gpu.ipynb
2020-04-06
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人