本文首发于“生信补给站”:https://mp.weixin.qq.com/s/lpkWwrLNtkLH8QA75X5STw 生存分析作为分析疾病/癌症预后的出镜频率超高的分析手段,而其结果展示的KM曲线也必须拥有姓名和颜值! 生存分析相关推文: 生存分析和KM曲线:R|生存分析(1) 分析结果一键输出:R|生存分析-结果整理 时间依赖生存分析:R|timeROC-分析 一 数据和R包 为方便,使用内置lung数据集 #载入所需的R包library("survival")library

R语言作为BI中ETL的工具,增删改 R语言提供了强大的R_package与各种数据库进行数据交互. 外加其强大数据变换清洗函数,为ETL提供一条方便快捷的道路. RODBC ROracal RMysql Rmongodb http://mirrors.ustc.edu.cn/CRAN/web/packages/rmongodb/vignettes/rmongodb_cheat_sheet.pdf step1 新建连接con,并查看其信息 library(RODBC) con<-odbcConn

概率图模型 基于R语言 这本书中的第一个R语言程序 prior <- c(working =0.99,broken =0.01) likelihood <- rbind(working = c(good=0.99,bad=0.01),broken =c(good=0.6,bad=0.4)) data <- c("bad","bad","bad","bad") bayes <- function(prio

目的: 1. 计算自定义模序在所有蛋白质的匹配位点和次数 2. 输出超过阈值的蛋白质序列到Hit_sequences.fasta 3. Hit_sequences.fasta中序列用小写字母,匹配用大写字母 4. 返回一个数据框,内容包存储ID.注释行(anno)括--.长度(len).匹配位置(Position),匹配次数(Hits),相应序列(tag) 一.问题思考: 1. 如何快速计算匹配位点 2. 输出文件如何构建 >序列ID(ACCESSION) 序列内容

目的:读入序列文件(fasta格式),返回一个数据框,内容包括--存储ID.注释行(anno).长度(len).序列内容(content) 一.问题思考: 1. 如何识别注释行和序列内容行 2. 如何快速定位序列内容所在位置 二.你可能需要的知识--基本的R语言基础 1. R语言基本数据类型 2. 会使用帮助(help,?)及网络资源 3. 其他的部分可能需要你针对自己看到的问题自己想办法解决或者留言 ##--构建函数--## seq_import <- function( file ){ se

问题: 使用forecast.Arima对带xreg的arima模型进行预测,报xreg Error pre.m4x <- forecast.Arima(m4x, h = 20, xreg = seq(429,448,by = 1)) 报错: Error in if (ncol(xreg) != ncol(object$call$xreg)) stop("Number of regressors does not match fitted model") : argument i

1. 在matlab中将数据保存到*.mat 文件夹 save("data.mat","data","label")#将data和label两个变量保存到data.mat文件夹中 2.在R语言中安装R.matlab包 install.packages('R.matlab') 3.读取*.mat 文件中的数据 library(R.matlab) ob<-readMat("data.mat")# 返回的是一个列表,通 $进行

方法一: 按照R导论中的方法,使用RODBC包, library(RODBC) channel<-odbcConnectExcel("file.xlsx") da2<-sqlFetch(channel,"Sheet1") odbcClose(channel) 结果: 未能成功,原因是 odbcConnectExcel is only usable with 32-bit Windows 即只能在32位的电脑中使用,64位的则不可以. 方法二: 使用xls

镜像下载.域名解析.时间同步请点击阿里云开源镜像站 方法一: 打开RStudio,输入options()$repos查看默认镜像源情况 options()$repos 打开tools工具栏,找到Global Options 找到Packages,点击右边的add 出现添加镜像源的对话框,name处填入aliyun,url填入:https://mirrors.aliyun.com/CRAN/, 点击ok保存 在命令栏输入options()$repos查看镜像源情况 方法二: 找到C盘,progra

1.用0替代数据框中的缺失值NA 生成数据框: > m <- matrix(sample(c(NA, :), , replace = TRUE), ) > d <- as.data.frame(m) V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 1 4 3 NA 3 7 6 6 10 6 5 2 9 8 9 5 10 NA 2 1 7 2 3 1 1 6 3 6 NA 1 4 1 6 4 NA 4 NA 7 10 2 NA 4 1 8 5 1 2 4 NA 2 6

判断ori_data[,1]中是否存在元素a: a %in% ori_data[,1] 如果存在返回 true,否则返回 false    

使用attach(file)时,一定要配合使用detach(file),否则再此运行程序时极易出现问题,The following objects are masked ... 此外工作空间中不能有与file里的变量同名的全局变量存在,如果有会显示masked attach()与with的使用方法与区别可参见http://www.cnblogs.com/Make-It-Fun/p/5265396.html

在R语言编译器中,设置当前工作文件夹可以用setwd()函数. > setwd("e://桌面//")> setwd("e:\桌面\")> setwd("e:/桌面/") 这三种结构都是可以编译通过的, 但是在VS C#中却不行,只有一种能运行成功. (PS:R语言在VS中运行要先配置环境,还没配置的童鞋先要配置好,才可运行,如有问题可看我前面的随笔.) 就是这种结构,engine.Evaluate("setwd('e

1.生存曲线 基础包survival+扩展包survminer. survival包内置肺癌数据集lung. library(survival) library(survminer) str(lung) #拟合模型 fit <- survfit(Surv(time,status)~sex,data=lung) #绘制生存曲线 ggsurvplot(fit, pval = TRUE, #添加log rank检验的p值 conf.int = TRUE, #添加置信区间 risk.table = TR

本文目录 公共参数列表 par 颜色相关 字体相关 字体大小相关 线条相关 符号相关 线条和符号大小相关 结束 本文首发:program-dog.blogspot.com 注1:本文也曾在csdn发布,不过无法忍受csdn超长时间的审核,迁移到博客圆了. 本作品采用知识共享署名-非商业性使用-相同方式共享 4.0 国际许可协议进行许可. 这一篇介绍par参数比较基础的几个参数用法,涉及颜色,字体,线条和符号,坐标轴,添加图例,组合做图留到下一篇文章. 上一篇文章已经详细的介绍了R语言可视化技术的

code{white-space: pre;} pre:not([class]) { background-color: white; }if (window.hljs && document.readyState && document.readyState === "complete") { window.setTimeout(function() { hljs.initHighlighting(); }, 0);}.main-container {

ddply和aggregate是两个用来整合数据的功能强大的函数. aggregate(x, ...) 关于aggregate()函数的使用在<R语言实战>中P105有简单描述,这里重新说一下.此函数主要有一下几种用法: ## Default S3 method: aggregate(x, ...) ## S3 method for class 'data.frame' aggregate(x, by, FUN, ..., simplify = TRUE, drop = TRUE) ## S3

Iris的R语言命令工具箱(1) 最近在做数据分析,使用了R语言做了些数据处理和可视化,在此记下遇到过的问题.应用过的命令.处理方式以及工具包- *版权声明:本文为博主原创文章,转载请注明本文地址.http://www.cnblogs.com/o0Iris0o/p/6365927.html* 1.CSV格式文件数据导入导出 --------------- ** 1.1 导入CSV数据(导入CSV数据之后储存为"Data Frame数据框"数据结构)** ## 导入csv数据 ## 第一

时间序列: (或称动态数列)是指将同一统计指标的数值按其发生的时间先后顺序排列而成的数列.时间序列分析的主要目的是根据已有的历史数据对未来进行预测.(百度百科) 主要考虑的因素: 1.长期趋势(Long-term trend) : 时间序列可能相当稳定或随时间呈现某种趋势. 时间序列趋势一般为线性的(linear),二次方程式的 (quadratic)或指数函数(exponential function). 2.季节性变动(Seasonal variation) 按时间变动,呈现重复性行为的序列

一.scale函数 R语言base库中自带数据标准化接口scale函数,函数介绍如下 Usage scale(x, center = TRUE, scale = TRUE) Arguments x: a numeric matrix(like object). center: either a logical value or a numeric vector of length equal to the number of columns of x. scale: either a logic

本文为带大家了解R语言以及分段式的步骤教程! 人们学习R语言时普遍存在缺乏系统学习方法的问题.学习者不知道从哪开始,如何进行,选择什么学习资源.虽然网络上有许多不错的免费学习资源,然而它们多过了头,反而会让人挑花了眼. 为了构建R语言学习方法,我们在Vidhya和DataCamp中选一组综合资源,帮您从头学习R语言.这套学习方法对于数据科学或R语言的初学者会很有用;如果读者是R语言的老用户,则会由本文了解这门语言的部分最新成果. R语言学习方法会帮助您快速.高效学习R语言. 前言 在开始学习之前