统计编程的框架与R语言统计分析基础——摘(1)

   2023-02-09 学习力0
核心提示:清屏命令ctrl+L一、基础1、产生数据结构a、直接输入b、冒号,1:10c、seq函数d、rep函数 1:10 [1]123456789 10 10:2[1] 1098765432 seq(1,10) [1]123456789 10 seq(1,10,2)[1] 1 3 5 7 9 rep(2,4)[1] 2 2 2 2 rep(1:3,times=3)[1] 1 2 3 1 2 3 1 2 3 rep(x=1:

清屏命令ctrl+L

一、基础

1、产生数据结构

a、直接输入

b、冒号,1:10

c、seq函数

d、rep函数

> 1:10
 [1]  1  2  3  4  5  6  7  8  9 10
> 10:2
[1] 10  9  8  7  6  5  4  3  2
> seq(1,10)
 [1]  1  2  3  4  5  6  7  8  9 10
> seq(1,10,2)
[1] 1 3 5 7 9
> rep(2,4)
[1] 2 2 2 2
> rep(1:3,times=3)
[1] 1 2 3 1 2 3 1 2 3
> rep(x=1:3, each=3)
[1] 1 1 1 2 2 2 3 3 3
> rep(1:3,1:3)
[1] 1 2 2 3 3 3

 

2、数据产生

a、c()函数产生向量

b、matrix()函数产生矩阵

c、data.frame()函数产生数据框

d、factor()函数产生因子

e、list()函数产生列表

f、ts()函数产生时间序列

> x = c(3,4,6)
> x
[1] 3 4 6
> matrix(1:10,2)
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    3    5    7    9
[2,]    2    4    6    8   10
> matrix(1:10,nrow=2,ncol=5,byrow=T)
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    2    3    4    5
[2,]    6    7    8    9   10
> x = data.frame(1:5, 6:10)
> x
  X1.5 X6.10
1    1     6
2    2     7
3    3     8
4    4     9
5    5    10
> x = cbind(x, c('a','b','c','d','e')) #绑上一列字符
> x
  X1.5 X6.10 c("a", "b", "c", "d", "e")
1    1     6                          a
2    2     7                          b
3    3     8                          c
4    4     9                          d
5    5    10                          e
> dimnames(x)                    #查看x的行列名
[[1]]
[1] "1" "2" "3" "4" "5"

[[2]]
[1] "X1.5"                                 "X6.10"                               
[3] "c(\"a\", \"b\", \"c\", \"d\", \"e\")"

> colnames(x)                    #只看列名
[1] "X1.5"                                 "X6.10"                               
[3] "c(\"a\", \"b\", \"c\", \"d\", \"e\")"
> colnames(x) = c('X1','X2','X3')        #改列名
> x
  X1 X2 X3
1  1  6  a
2  2  7  b
3  3  8  c
4  4  9  d
5  5 10  e

 

3、运算

%%  余数

%/%  整数商

^    乘方

& | !  逻辑运算

 

4、下标的使用(获取元素)  用中括号[]

> x = 1:10
> x
 [1]  1  2  3  4  5  6  7  8  9 10
> x[x>6]
[1]  7  8  9 10
> x[x>6 & x<9]
[1] 7 8
> x = matrix(1:20, 4)
> x
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    5    9   13   17
[2,]    2    6   10   14   18
[3,]    3    7   11   15   19
[4,]    4    8   12   16   20
> x[x>2 & x<14]
 [1]  3  4  5  6  7  8  9 10 11 12 13
> x[x>2 & x<14] = NA
> x
     [,1] [,2] [,3] [,4] [,5]
[1,]    1   NA   NA   NA   17
[2,]    2   NA   NA   14   18
[3,]   NA   NA   NA   15   19
[4,]   NA   NA   NA   16   20

 

二、一些数字和统计函数

1、

max()     min()    mean()    

标准差sd()    方差var()    相关系数cor()

求和sum()    求积prod()

中位数median()   分位数quantile()

对数log()    指数exp()

排列factorial()  组合choose()

四舍五入round()  向下取整floor()  向上取整ceiling()

总结summary()

2、

累加cumsum()  秩rank()  排序sort()  倒序rev()

矩阵转置t()  逆矩阵solve()  特征根eigen()

 

三、流程控制

> x = 6
> while(x>1 & winDialog("yesno",message="you are sb")!='YES')
+ x = x-1

 

 
反对 0举报 0 评论 0
 

免责声明:本文仅代表作者个人观点,与乐学笔记(本网)无关。其原创性以及文中陈述文字和内容未经本站证实,对本文以及其中全部或者部分内容、文字的真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
    本网站有部分内容均转载自其它媒体,转载目的在于传递更多信息,并不代表本网赞同其观点和对其真实性负责,若因作品内容、知识产权、版权和其他问题,请及时提供相关证明等材料并与我们留言联系,本网站将在规定时间内给予删除等相关处理.

  • 拓端tecdat|R语言VAR模型的不同类型的脉冲响应
    原文链接:http://tecdat.cn/?p=9384目录模型与数据估算值预测误差脉冲响应识别问题正交脉冲响应结构脉冲反应广义脉冲响应参考文献脉冲响应分析是采用向量自回归模型的计量经济学分析中的重要一步。它们的主要目的是描述模型变量对一个或多个变量的冲击的演化
    03-16
  • Visual Studio 编辑R语言环境搭建
    Visual Studio 编辑R语言环境搭建关于Visual Studio 编辑R语言环境搭建具体的可以看下面三个网址里的内容,我这里就讲两个问题,关于r包管理和换本地的r的服务。1.r包管理:Ctrl+72.R本地服务管理:Ctrl+9Visual Studio R官方帮助文档(中文): https://docs
    03-16
  • 拓端tecdat|R语言代写实现向量自回归VAR模型
    原文链接:http://tecdat.cn/?p=8478 澳大利亚在2008 - 2009年全球金融危机期间发生了这种情况。澳大利亚政府发布了一揽子刺激计划,其中包括2008年12月的现金支付,恰逢圣诞节支出。因此,零售商报告销售强劲,经济受到刺激。因此,收入增加了。VAR面临的批
    03-16
  • [译]用R语言做挖掘数据《五》 r语言数据挖掘简
    一、实验说明1. 环境登录无需密码自动登录,系统用户名shiyanlou,密码shiyanlou2. 环境介绍本实验环境采用带桌面的Ubuntu Linux环境,实验中会用到程序:1. LX终端(LXTerminal): Linux命令行终端,打开后会进入Bash环境,可以使用Linux命令2. GVim:非常好
    03-08
  • 拓端tecdat|Mac系统R语言升级后无法加载包报错 package or namespace load failed in dyn.load(file, DLLpath = DLLpath, ..
    拓端tecdat|Mac系统R语言升级后无法加载包报错
    问题重现:我需要安装R软件包stochvol,该软件包 仅适用于3.6.0版的R。因此,我安装了R(3.6.0 版本),并使用打开它 RStudio。但是现在  ,即使我成功 使用来 安装软件包,也无法加载任何库 。具体来说,我需要加载的库是stochvol  ,Rcpp和 caret
    03-08
  • 拓端数据tecdat|R语言k-means聚类、层次聚类、主成分(PCA)降维及可视化分析鸢尾花iris数据集
    拓端数据tecdat|R语言k-means聚类、层次聚类、
    原文链接:http://tecdat.cn/?p=22838 原文出处:拓端数据部落公众号问题:使用R中的鸢尾花数据集(a)部分:k-means聚类使用k-means聚类法将数据集聚成2组。 画一个图来显示聚类的情况使用k-means聚类法将数据集聚成3组。画一个图来显示聚类的情况(b)部分:
    03-08
  • 《R语言数据挖掘》读书笔记:七、离群点(异常值)检测
    《R语言数据挖掘》读书笔记:七、离群点(异常值
    第七章、异常值检测(离群点挖掘)概述:        一般来说,异常值出现有各种原因,比如数据集因为数据来自不同的类、数据测量系统误差而收到损害。根据异常值的检测,异常值与原始数据集中的常规数据显著不同。开发了多种解决方案来检测他们,其中包括
    03-08
  • 拓端数据tecdat|R语言中实现广义相加模型GAM和普通最小二乘(OLS)回归
    拓端数据tecdat|R语言中实现广义相加模型GAM和
    原文链接:http://tecdat.cn/?p=20882  1导言这篇文章探讨了为什么使用广义相加模型 是一个不错的选择。为此,我们首先需要看一下线性回归,看看为什么在某些情况下它可能不是最佳选择。 2回归模型假设我们有一些带有两个属性Y和X的数据。如果它们是线性
    03-08
  • 拓端数据tecdat|R语言时间序列平稳性几种单位根检验(ADF,KPSS,PP)及比较分析
    拓端数据tecdat|R语言时间序列平稳性几种单位根
    原文链接:http://tecdat.cn/?p=21757 时间序列模型根据研究对象是否随机分为确定性模型和随机性模型两大类。随机时间序列模型即是指仅用它的过去值及随机扰动项所建立起来的模型,建立具体的模型,需解决如下三个问题模型的具体形式、时序变量的滞后期以及随
    03-08
  • 拓端tecdat|R语言风险价值VaR(Value at Risk)和损失期望值ES(Expected shortfall)的估计
    拓端tecdat|R语言风险价值VaR(Value at Risk)
    原文链接: http://tecdat.cn/?p=15929 风险价值VaR和损失期望值ES是常见的风险度量。首先明确:时间范围-我们展望多少天?概率水平-我们怎么看尾部分布?在给定时间范围内的盈亏预测分布,示例如图1所示。  图1:预测的损益分布 给定概率水平的预测的分
    03-08
点击排行