- 博客(12)
- 资源 (5)
- 问答 (2)
- 收藏
- 关注
RSS订阅原创 R语言二次规划之quadprog包
二次规划二次规划是非线性规划中一类特殊的数学规划问题，在如投资组合、约束最小二乘问题中都有应用。二次规划已经成为运筹学，经济数学，管理科学和组合优化科学的基本方法。 二次规划的一般形式： minxq(x)=12xTDx−dTx+cs.t.ATx≥b0min_{x}\quad q(x)=\frac{1}{2}x^TDx-d^Tx+c\\s.t. \quad A^Tx\ge b_0 其中，x为向
2016-12-30 22:04:31
15385
4
原创 python小函数（一）
1. shape函数numpy模块 shape(a) - a: 数组np.shape(np.eye(3)) ## 返回结果 (3L, 3L)np.shape([[1, 2]]) ## 返回结果(1L, 2L)np.shape([0]) ##返回元组## (1L,)shape(0) ##返回空元组## ()a = array([[3,4,5],[6,3,6]
2016-12-28 15:14:31
993
原创 长城评论词云分析
最近《长城》这部电影很火，争议也很大，我根据豆瓣底部的点评来分析这部电影在观众眼中是怎样的。 此为长城短评链接。 下面给出代码：library(XML)library(RCurl)library(stringr)library(Rwordseg)library(tm)library(wordcloud2)library(wordcloud)url <- "https://movie
2016-12-24 22:20:43
670
原创 决策树（一）
1. 特征选择1.1 信息增益熵（entropy）是表示随机变量不确定性的度量，设x是一个取有限个值的离散随机变量，其概率分布为： P(X=xi)=pi,i=1,2,⋯,nP(X=x_i)=p_i,i=1,2,\cdots,n 则随机变量X的熵定义为 H(X)=−∑i=1npilogpiH(X)=-\sum_{i=1}^np_ilog p_i 熵只依赖于X的分布，与X的取值无关，所
2016-12-22 16:54:27
615
原创 K近邻法算法（KNN）及其R实现
1. K近邻算法输入：训练数据集 T={(x1,y1),(x2,y2),⋯,(xN,yN)}T=\{(x_1,y_1),(x_2,y_2),\cdots,(x_N,y_N)\} 其中，xi∈χ⊆Rnx_i \in \chi \subseteq R^n为实例的特征向量，yi∈Y={c1,c2,⋯,cK}y_i\in Y=\{c_1,c_2,\cdots,c_K\}为实例的类别，i=1,2,⋯,N;
2016-12-21 20:50:51
9351
原创 网络数据抓取及其R实现（以链家楼盘为例）
\quad本次数据抓取的链家新房楼盘，由于对网络抓取数据还不熟悉，所以现在只抓取了北京地区的楼盘。我在学习过程中主要参考资料为《基于R语言的自动数据收集》，这本书确实不错，在前八章讲解的基础知识，以本书的配套网站资源为讲解基础,但是在按照书中抓取数据时有些地方和书本讲的不太一样，我认为主要是网络结构的变化引起的，大部分还是可行的。现在我看到了第九章，再往后都是一些网络数据抓取项目实例，所以后面的东西
2016-12-19 21:39:01
4103
原创 朴素贝叶斯法及其R实现
1. 朴素贝叶斯基本方法1.1 贝叶斯分类法基本公式：P(Y=ck|X=x)=P(Y=ck)ΠjP(X(j)=x(j)|Y=ck)∑k(P(Y=ck)ΠjP(Xj=xj|Y=ck)),k=1,2,⋯,KP(Y=c_k|X=x)=\frac{P(Y=c_k)\Pi_jP(X^{(j)}=x^{(j)}|Y=c_k)}{\sum_k(P(Y=c_k)\Pi_jP(X^{j}=x^{j}|Y=c_k))}
2016-12-19 12:44:35
6255
原创 感知机及其R实现
感知机1. 感知机模型定义：假设输入空间（特征空间）是χ⊆Rn\chi\subseteq R^n，输出空间是Y={+1,−1}Y=\{+1,-1\}。输入x∈χx\in \chi表示实例的特征向量，对应于输入空间的点；输出y∈Yy\in Y表示实例的类别。由输入空间到输出空间的如下函数 f(x)=sign(w⋅x+b)f(x)=sign(w\cdot x+b) 称为感知机。其中，ww和bb为感知
2016-12-15 19:14:15
5573
1
原创 R语言 MCMC算法及其实现
Metropolis-Hasting算法：Step 1: 选择一个不可约Markov链转移概率q(i,j),i,j∈\inS.再从S={1，2，… ,m}中选择某个整数. Step 2: 令n=0以及x0=kx_0=k. Step 3: 产生一个随机变量X使得P{X=j}=q(Xn,j)P\{X=j\}=q(X_n,j), 再产生一个随机数U. Step 4: 如果U
2016-12-12 19:09:47
15512
6
原创 R语言 牛顿-拉夫森迭代法求方程组
R语言 牛顿-拉夫森迭代法求方程组牛顿-拉夫森迭代法： xk+1=xk−[f′(x)]−1f(x)x_{k+1}=x_{k}-{[f'(x)]^{-1}}f(x) 其中，f′(x)f'(x)为Jacobi矩阵。 例，设非线性方程组为： x2+y2=1,x^2+y^2=1, y2=2xy^2=2x 求方程组的解。 Jocabi行列式： R代码如下：fun
2016-12-11 23:00:30
15913
3
c语言控制台输出中文乱码问题
2018-03-28
非法类型开始与标识符
2017-08-07
TA创建的收藏夹 TA关注的收藏夹
TA关注的人