- 找不到某些GEO数据的表达矩阵肿么办
小洁忘了怎么分身
举个栗子,GSE99410和GSE78200。library(GEOquery)gse="GSE78200"eSet<-getGEO(gse,destdir='.',getGPL=F)exp<-exprs(eSet[[1]])#Errorinexp[1:4,1:4]:subscriptoutofboundsexp[1:4,1:4]boxplot(exp)#Errorinnames(groups)<
- R语言学习-GEO数据挖掘
米粥的味道
根据教程准备试着处理一下GSE的其他数据,发现对于多组数据很难简单的复现,待后续探索完善再补更;先复习LINUX;(1)下载数据及读取rm(list=ls())options(stringsAsFactors=F)library(GEOquery)eSet<-getGEO("GSE94648",destdir='.',getGPL=F)exp<-exprs(eSet[[1]])exp[1:4,1:
- 7+非肿瘤+WGCNA+机器学习+诊断模型,构思巧妙且操作简单
生信风暴
论文阅读
今天给同学们分享一篇生信文章“Platelets-relatedsignaturebaseddiagnosticmodelinrheumatoidarthritisusingWGCNAandmachinelearning”,这篇文章发表在FrontImmunol期刊上,影响因子为7.3。结果解读:DEGs和血小板相关基因的鉴定作者通过R的“GEOquery”和“limma”包,在数据预处理后获得了
- GEOquery 下载 GEO 数据
名本无名
前言NCBIGeneExpressionOmnibus(基因表达综合数据库,GEO)公开了很多高通量基因组学数据集。尽管名字很好听,但这个数据库并不是专门针对基因表达数据的。NCBIGEONCBIGEO是以样本的形式组织起来的,这些样本被归入series。对于较大的实验,有SubSeries和SuperSeries两种一个SuperSeries就是一篇论文的所有实验,一个SuperSeries可以
- 2018-10-11
一路向前_莫问前程_前程似锦
使用GEOquery包if(T){source("http://bioconductor.org/biocLite.R")options(BioC_mirror="http://mirrors.ustc.edu.cn/bioc/")biocLite("GEOquery")library(GEOquery)eSeteSet$GSE62832_series_matrix.txt.gzExpressio
- R 报错:ERROR: failed to lock directory XXX 00LOCK
11的雾
Error:install.packagesfailedtolockdirectory在使用bioconductor安装一些包的时候,比如我安装GEOquery时,会有ErrorinunpackPkgZip(foundpkgs[okp,2L],foundpkgs[okp,1L],lib,libs_only,:ERROR:failedtolockdirectory‘C:\ProgramFiles\R
- GEO通用
pedagogoose
数据下载:可以通过GEOquery这个R包进行下载,或者直接在网页上下载---下载的内容包括:==表达矩阵==和==注释文件==rm(list=ls())#一键清空环境中的对象options(stringsAsFactors=F)#默认的设置library(GEOquery)gse="GSE107943"#需要找到需要下载的数据的GSE号eSet<-getGEO(gse,destdir='.',g
- aglient芯片原始数据处理
小洁忘了怎么分身
本文讲的是aglient芯片原始数据的处理,参考资料是limma的userguide文档。GEO数据库下载的表达矩阵不符合预期,比如是空的,或者是有负值的,那我们就处理一下它的原始数据。aglient的芯片应用也很广泛,举个OSCC的栗子:GSE23558,跟着学习学习。1.下载和读取数据1.1获取临床信息数据从前,提到GEO数据下载,我们只有GEOquery,神功盖世,但是死于网速。后来就有了中
- GEO数据挖掘--下载数据的3种方式
Seurat_
rm(list=ls())##魔幻操作,一键清空~当前环境中对象全部删除options(stringsAsFactors=F)#在调用as.data.frame的时,将stringsAsFactors设置为FALSE可以避免character类型自动转化为factor类型数据下载下载方式1GEOquery#library(GEOquery)gse_number="GSE42872"#修改号码#eS
- GEO生信数据挖掘(五)提取临床信息构建分组,分组数据可视化(绘制层次聚类图,绘制PCA图)
人工智能学术前沿(真)
数据挖掘人工智能生信分析GEO数据挖掘基因数据分析
检索到目标数据集后,开始数据挖掘,本文以阿尔兹海默症数据集GSE1297为例上节做了很多的基因数据清洗(离群值处理、低表达基因、归一化、log2处理)操作,本节介绍构建临床分组信息。我们已经学习了提取表达矩阵的临床信息#安装并加载GEOquery包library(GEOquery)#指定GEO数据集的IDgse_id与是一致的。需要的分组信息已经提取完毕。分组后箱线图可视化上一节保存了数据清洗后的
- GEO生信数据挖掘(一)数据集下载和初步观察
人工智能学术前沿(真)
生信分析GEO数据挖掘基因数据分析
检索到目标数据集后,开始数据挖掘,本文以阿尔兹海默症数据集GSE1297为例目录GEOquery简介安装并加载GEOquery包getGEO函数获取数据(联网下载)更换下载数据源对数据集进行初步观察处理GEOquery简介GEOquery是一个在生物信息学中常用的R语言包,用于从NCBIGeneExpressionOmnibus(GEO)数据库中获取和分析基因表达数据。以下是GEOquery包的简
- 数据挖掘学习笔记——GEO数据库:芯片数据分析
福旺旺
生物信息学数据挖掘学习数据库
数据挖掘数据挖掘学习笔记——GEO数据库:芯片数据分析文章目录数据挖掘一、芯片基础知识1.1、背景二、GEO数据库概述2.1、基础简介2.2、检索页面展示三、GSE项目的三种下载方式3.1、主页下载原始数据3.2、主页下载表达矩阵3.3、GEOquery包下载表达量四、基因名与探针ID的转换技巧4.1、获取对照关系4.1.1、利用Bioconductor中汇总的R包4.1.2、利用平台的数据4.1
- 190304 用GEOquery从GEO数据库下载数据
森尼啊
①下载得到对象gds858<-getGEO(‘GDS858’,destdir=“.”)##根据GDS号来下载数据,下载soft文件gpl96<-getGEO(‘GPL96’,destdir=“.”)##根据GPL号下载的是芯片设计的信息!gse1009<-getGEO(‘GSE1009’,destdir=“.”)##根据GSE号下载数据,下载_series_matrix.txt.gz下载的文件都会
- 从GEO数据库下载数据并处理
刘大帅1
rm(list=ls())Versioninfo:R3.2.3,Biobase2.30.0,GEOquery2.40.0,limma3.26.8RscriptsgeneratedSunAug1123:08:50EDT2019################################################################rm(list=ls())#使用limma包对数
- 从GEO数据库下载数据
YTXu
从GEO数据库下载数据有很多种方法,主要用的是jimmy在生信技能树里给出的和果子老师现成的代码。Jimmy的更变通,果子老师的就只教了一种傻瓜式方法。下面是jimmy的,Ref:http://www.bio-info-trainee.com/bioconductor_China/software/GEOquery.html只需要记住三个函数,以及每个函数返回的对象该如何处理即可getGEO/ge
- R GEOquery: Timeout of 60 seconds was reached 解决办法
生信摆渡
`感谢各位简友一直的关注。请关注:感谢陪伴``RgetOption("timeout")options(timeout=1000)默认下载时限为60秒,修改为1000秒即可
- GEO数据挖掘
看远方的星
复现此篇文章:https://mp.weixin.qq.com/s/5Bdx7YyAk8InYs2MFIUUVQ下载相应的R包options(BioC_mirror="http://mirrors.ustc.edu.cn/bioc/")#设置镜像BiocManager::install('stringi')BiocManager::install('GEOquery')BiocManager::i
- GEOquery包、affy包、DESeq2包的加载
研平方
GEOquery包的加载library(Biobase)library(BiocGenerics)library(parallel)library(GEOquery)affy包的加载library(BiocGenerics)library(parallel)library(Biobase)library(affy)DESeq2包的加载library(S4Vectors)library(stats4
- 2019-12-15小白新手笔记 GEO下载数据及初步处理可视化,差异表达
海阔天空周
今天介绍另外一种方式具体哪种好用,我还在检测中rm(list=ls())options(stringsAsFactors=F)library(GEOquery)setwd("C:\\Users\\zhouwenqing789\\Desktop\\DFEF")eset1<-getGEO("GSE83521",destdir=".",getGPL=F)#下载了平台注释文件以及不是矩阵的形式,list。
- 2020-12-28 GEO下载代码
海阔天空周
setwd("D:\\GSE140358衰弱")library(GEOquery)getGEO(GEO="GSE140358",)gset140358<-getGEO('GSE140358',GSEMatrix=TRUE,destdir=".",AnnotGPL=T,getGPL=T)save(gset140358,file="GSE140358.Rdata")#取列表中的元素a=gset1403
- 2019-07-06学习总结
桃浪桃浪
GEOquery包使用说明书https://www.jianshu.com/p/0d32fd410bcf作业6>rm(list=ls())##魔幻操作,一键清空~>options(stringsAsFactors=F)>f#https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE17215>library(GEOquery)>#这个包需要注意两个
- linux学习100篇8:ERROR: configuration failed for package 'curl'
Seurat_Satija
ERROR:dependencies'httr','xml2'arenotavailableforpackage'GEOquery'*removing'/usr/local/lib/R/site-library/GEOquery'Thedownloadedsourcepackagesarein'/tmp/RtmpkVS6KY/downloaded_packages'Warningmessages:
- 【install.packages(“openssl“)不成功,提示tools/version.c:1:10: fatal error: openssl/opensslv.h: 不存在】
weixin_42001274
ubuntulinux运维
问题描述:安装GEOquery这个R包过程中,显示依赖的openssl的R包安装不成功.环境中已经有安装openssl.错误信息如下:BiocManager::install("GEOquery")tryingURL'https://bioconductor.org/packages/3.10/bioc/src/contrib/GEOquery_2.54.1.tar.gz'Contenttype'
- 2020-02-10
Clouddong
GSE125042+WGCNA数据下载+预处理rm(list=ls())options(stringsAsFactors=F)library(GEOmirror)library(AnnoProbe)library(WGCNA)library(GEOquery)gset1)idxhead(g_l)#gene_idlength#11000096001842#21000096092538#3100009
- Ubuntu20.04 安装 DeSeq2报错解决
Angeladaddy
问题描述ubuntu20.04,安装好了R(4.2.0)和RstudioServer,按照如下脚本安装包时报错:if(!requireNamespace("BiocManager",quietly=TRUE))install.packages("BiocManager")BiocManager::install("DESeq2")BiocManager::install("GEOquery")Bi
- 通过3个数据集取差异基因的交集复现一篇文中的韦恩图
小梦游仙境
老大让复现的是一张韦恩图,来自3个数据集,分别找每个数据集的差异基因,然后取UP和DOWN的交集。原文链接:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6054764/#.原文中的图片如下image-20191121212758565想先用循环下载含有表达矩阵和临床信息的一个list,然后试着用循环library(GEOquery)ls2]=2n[n0
- 2019-12-13 小白新手笔记 GEO数据库下载及ID转换
海阔天空周
GEO数据库下载library(Biobase)library(BiocGenerics)library(parallel)library(GEOquery)gest<-getGEO("GSE9579",destdir=".",GSEMatrix=T,AnnotGPL=FALSE,getGPL=F)#"."可填保存的路径从GPL平台直接下载平台基因信息gpl570=read.delim("C:\\
- limma | 分层样本的差异分析这样搞!(二)
生信漫卷
后端
1写在前面上期介绍了用limma包做配对样本的差异分析。本期介绍一下Multi-level如何处理吧。应用场景:Control和Diseased的T细胞和B细胞分层对比。2用到的包rm(list=ls())library(tidyverse)library(limma)library(GEOquery)3示例数据这里我们还是利用上期介绍的GEO数据库上的dataset。在3个样本中对T细胞和B细胞
- 生信技能树——GEO芯片数据的合并
一指流年,一纸沙
生信技能树r语言
GSE83521和GSE89143数据合并1.下载数据rm(list=ls())library(GEOquery)library(stringr)gse="GSE83521"eSet1<-getGEO("GSE83521",destdir='.',getGPL=F)eSet2<-getGEO("GSE89143",destdir='.',getGPL=F)#(1)提取表达矩阵expexp1<-ex
- 5. 生信技能树——GEO转录组RNA_seq_GSE162550
一指流年,一纸沙
生信技能树R语言GEOr语言
和生信技能树GEO转录组“GSE150392“分析类似,唯一区别就是在数据处理和ID转换这一环节略微有区别1.数据下载最方便的是xena。可以网页下载,也可以用代码下载。proj="DHA"2.生存信息与临床信息这里仅仅是查看一下,到生存信息部分再整理。library(GEOquery)eSet=getGEO("GSE162550",destdir=".",getGPL=F)eSet=eSet[[
- 继之前的线程循环加到窗口中运行
3213213333332132
javathreadJFrameJPanel
之前写了有关java线程的循环执行和结束,因为想制作成exe文件,想把执行的效果加到窗口上,所以就结合了JFrame和JPanel写了这个程序,这里直接贴出代码,在窗口上运行的效果下面有附图。
package thread;
import java.awt.Graphics;
import java.text.SimpleDateFormat;
import java.util
- linux 常用命令
BlueSkator
linux命令
1.grep
相信这个命令可以说是大家最常用的命令之一了。尤其是查询生产环境的日志,这个命令绝对是必不可少的。
但之前总是习惯于使用 (grep -n 关键字 文件名 )查出关键字以及该关键字所在的行数,然后再用 (sed -n '100,200p' 文件名),去查出该关键字之后的日志内容。
但其实还有更简便的办法,就是用(grep -B n、-A n、-C n 关键
- php heredoc原文档和nowdoc语法
dcj3sjt126com
PHPheredocnowdoc
<!doctype html>
<html lang="en">
<head>
<meta charset="utf-8">
<title>Current To-Do List</title>
</head>
<body>
<?
- overflow的属性
周华华
JavaScript
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml&q
- 《我所了解的Java》——总体目录
g21121
java
准备用一年左右时间写一个系列的文章《我所了解的Java》,目录及内容会不断完善及调整。
在编写相关内容时难免出现笔误、代码无法执行、名词理解错误等,请大家及时指出,我会第一时间更正。
&n
- [简单]docx4j常用方法小结
53873039oycg
docx
本代码基于docx4j-3.2.0,在office word 2007上测试通过。代码如下:
import java.io.File;
import java.io.FileInputStream;
import ja
- Spring配置学习
云端月影
spring配置
首先来看一个标准的Spring配置文件 applicationContext.xml
<?xml version="1.0" encoding="UTF-8"?>
<beans xmlns="http://www.springframework.org/schema/beans"
xmlns:xsi=&q
- Java新手入门的30个基本概念三
aijuans
java新手java 入门
17.Java中的每一个类都是从Object类扩展而来的。 18.object类中的equal和toString方法。 equal用于测试一个对象是否同另一个对象相等。 toString返回一个代表该对象的字符串,几乎每一个类都会重载该方法,以便返回当前状态的正确表示.(toString 方法是一个很重要的方法) 19.通用编程:任何类类型的所有值都可以同object类性的变量来代替。
- 《2008 IBM Rational 软件开发高峰论坛会议》小记
antonyup_2006
软件测试敏捷开发项目管理IBM活动
我一直想写些总结,用于交流和备忘,然都没提笔,今以一篇参加活动的感受小记开个头,呵呵!
其实参加《2008 IBM Rational 软件开发高峰论坛会议》是9月4号,那天刚好调休.但接着项目颇为忙,所以今天在中秋佳节的假期里整理了下.
参加这次活动是一个朋友给的一个邀请书,才知道有这样的一个活动,虽然现在项目暂时没用到IBM的解决方案,但觉的参与这样一个活动可以拓宽下视野和相关知识.
- PL/SQL的过程编程,异常,声明变量,PL/SQL块
百合不是茶
PL/SQL的过程编程异常PL/SQL块声明变量
PL/SQL;
过程;
符号;
变量;
PL/SQL块;
输出;
异常;
PL/SQL 是过程语言(Procedural Language)与结构化查询语言(SQL)结合而成的编程语言PL/SQL 是对 SQL 的扩展,sql的执行时每次都要写操作
- Mockito(三)--完整功能介绍
bijian1013
持续集成mockito单元测试
mockito官网:http://code.google.com/p/mockito/,打开documentation可以看到官方最新的文档资料。
一.使用mockito验证行为
//首先要import Mockito
import static org.mockito.Mockito.*;
//mo
- 精通Oracle10编程SQL(8)使用复合数据类型
bijian1013
oracle数据库plsql
/*
*使用复合数据类型
*/
--PL/SQL记录
--定义PL/SQL记录
--自定义PL/SQL记录
DECLARE
TYPE emp_record_type IS RECORD(
name emp.ename%TYPE,
salary emp.sal%TYPE,
dno emp.deptno%TYPE
);
emp_
- 【Linux常用命令一】grep命令
bit1129
Linux常用命令
grep命令格式
grep [option] pattern [file-list]
grep命令用于在指定的文件(一个或者多个,file-list)中查找包含模式串(pattern)的行,[option]用于控制grep命令的查找方式。
pattern可以是普通字符串,也可以是正则表达式,当查找的字符串包含正则表达式字符或者特
- mybatis3入门学习笔记
白糖_
sqlibatisqqjdbc配置管理
MyBatis 的前身就是iBatis,是一个数据持久层(ORM)框架。 MyBatis 是支持普通 SQL 查询,存储过程和高级映射的优秀持久层框架。MyBatis对JDBC进行了一次很浅的封装。
以前也学过iBatis,因为MyBatis是iBatis的升级版本,最初以为改动应该不大,实际结果是MyBatis对配置文件进行了一些大的改动,使整个框架更加方便人性化。
- Linux 命令神器:lsof 入门
ronin47
lsof
lsof是系统管理/安全的尤伯工具。我大多数时候用它来从系统获得与网络连接相关的信息,但那只是这个强大而又鲜为人知的应用的第一步。将这个工具称之为lsof真实名副其实,因为它是指“列出打开文件(lists openfiles)”。而有一点要切记,在Unix中一切(包括网络套接口)都是文件。
有趣的是,lsof也是有着最多
- java实现两个大数相加,可能存在溢出。
bylijinnan
java实现
import java.math.BigInteger;
import java.util.regex.Matcher;
import java.util.regex.Pattern;
public class BigIntegerAddition {
/**
* 题目:java实现两个大数相加,可能存在溢出。
* 如123456789 + 987654321
- Kettle学习资料分享,附大神用Kettle的一套流程完成对整个数据库迁移方法
Kai_Ge
Kettle
Kettle学习资料分享
Kettle 3.2 使用说明书
目录
概述..........................................................................................................................................7
1.Kettle 资源库管
- [货币与金融]钢之炼金术士
comsci
金融
自古以来,都有一些人在从事炼金术的工作.........但是很少有成功的
那么随着人类在理论物理和工程物理上面取得的一些突破性进展......
炼金术这个古老
- Toast原来也可以多样化
dai_lm
androidtoast
Style 1: 默认
Toast def = Toast.makeText(this, "default", Toast.LENGTH_SHORT);
def.show();
Style 2: 顶部显示
Toast top = Toast.makeText(this, "top", Toast.LENGTH_SHORT);
t
- java数据计算的几种解决方法3
datamachine
javahadoopibatisr-languer
4、iBatis
简单敏捷因此强大的数据计算层。和Hibernate不同,它鼓励写SQL,所以学习成本最低。同时它用最小的代价实现了计算脚本和JAVA代码的解耦,只用20%的代价就实现了hibernate 80%的功能,没实现的20%是计算脚本和数据库的解耦。
复杂计算环境是它的弱项,比如:分布式计算、复杂计算、非数据
- 向网页中插入透明Flash的方法和技巧
dcj3sjt126com
htmlWebFlash
将
Flash 作品插入网页的时候,我们有时候会需要将它设为透明,有时候我们需要在Flash的背面插入一些漂亮的图片,搭配出漂亮的效果……下面我们介绍一些将Flash插入网页中的一些透明的设置技巧。
一、Swf透明、无坐标控制 首先教大家最简单的插入Flash的代码,透明,无坐标控制: 注意wmode="transparent"是控制Flash是否透明
- ios UICollectionView的使用
dcj3sjt126com
UICollectionView的使用有两种方法,一种是继承UICollectionViewController,这个Controller会自带一个UICollectionView;另外一种是作为一个视图放在普通的UIViewController里面。
个人更喜欢第二种。下面采用第二种方式简单介绍一下UICollectionView的使用。
1.UIViewController实现委托,代码如
- Eos平台java公共逻辑
蕃薯耀
Eos平台java公共逻辑Eos平台java公共逻辑
Eos平台java公共逻辑
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
蕃薯耀 2015年6月1日 17:20:4
- SpringMVC4零配置--Web上下文配置【MvcConfig】
hanqunfeng
springmvc4
与SpringSecurity的配置类似,spring同样为我们提供了一个实现类WebMvcConfigurationSupport和一个注解@EnableWebMvc以帮助我们减少bean的声明。
applicationContext-MvcConfig.xml
<!-- 启用注解,并定义组件查找规则 ,mvc层只负责扫描@Controller -->
<
- 解决ie和其他浏览器poi下载excel文件名乱码
jackyrong
Excel
使用poi,做传统的excel导出,然后想在浏览器中,让用户选择另存为,保存用户下载的xls文件,这个时候,可能的是在ie下出现乱码(ie,9,10,11),但在firefox,chrome下没乱码,
因此必须综合判断,编写一个工具类:
/**
*
* @Title: pro
- 挥洒泪水的青春
lampcy
编程生活程序员
2015年2月28日,我辞职了,离开了相处一年的触控,转过身--挥洒掉泪水,毅然来到了兄弟连,背负着许多的不解、质疑——”你一个零基础、脑子又不聪明的人,还敢跨行业,选择Unity3D?“,”真是不自量力••••••“,”真是初生牛犊不怕虎•••••“,••••••我只是淡淡一笑,拎着行李----坐上了通向挥洒泪水的青春之地——兄弟连!
这就是我青春的分割线,不后悔,只会去用泪水浇灌——已经来到
- 稳增长之中国股市两点意见-----严控做空,建立涨跌停版停牌重组机制
nannan408
对于股市,我们国家的监管还是有点拼的,但始终拼不过飞流直下的恐慌,为什么呢?
笔者首先支持股市的监管。对于股市越管越荡的现象,笔者认为首先是做空力量超过了股市自身的升力,并且对于跌停停牌重组的快速反应还没建立好,上市公司对于股价下跌没有很好的利好支撑。
我们来看美国和香港是怎么应对股灾的。美国是靠禁止重要股票做空,在
- 动态设置iframe高度(iframe高度自适应)
Rainbow702
JavaScriptiframecontentDocument高度自适应局部刷新
如果需要对画面中的部分区域作局部刷新,大家可能都会想到使用ajax。
但有些情况下,须使用在页面中嵌入一个iframe来作局部刷新。
对于使用iframe的情况,发现有一个问题,就是iframe中的页面的高度可能会很高,但是外面页面并不会被iframe内部页面给撑开,如下面的结构:
<div id="content">
<div id=&quo
- 用Rapael做图表
tntxia
rap
function drawReport(paper,attr,data){
var width = attr.width;
var height = attr.height;
var max = 0;
&nbs
- HTML5 bootstrap2网页兼容(支持IE10以下)
xiaoluode
html5bootstrap
<!DOCTYPE html>
<html>
<head lang="zh-CN">
<meta charset="UTF-8">
<meta http-equiv="X-UA-Compatible" content="IE=edge">