porco_bosso-优快云博客

原创 RMysql

安装时注意mysql为32位，rjava也要32位。#连接到服务器ucscDb<-dbConnect(MySQL(),user="genome",host="genome-mysql.cse.ucsc.edu") #获得数据库 result<-dbGetQuery(ucscDb,"show databases;") #连接到指定数据库 hg19<-dbConnect(MySQL(),user="

2014-04-22 20:14:36 730

原创 Rmysql

安装时注意mysql为32位，#连接到服务器ucscDb<-dbConnect(MySQL(),user="genome",host="genome-mysql.cse.ucsc.edu") #获得数据库 result<-dbGetQuery(ucscDb,"show databases;") #连接到指定数据库 hg19<-dbConnect(MySQL(),user="genome",db=

2014-04-22 20:14:09 687

source("http://bioconductor.org/biocLite.R") biocLite("rhdf5") library(rhdf5) create = h5createFile("example.h5") getwd() create=h5createGroup("example.h5","foo") create=h5createGroup("example.h5","ba

2014-04-22 20:11:22 1237

原创 week 2

GIT 下载git:http://git-scm.com/downloads 安装git:默认下一步使用: 1.设置用户名，邮箱： git config --global user.name "" git config --global user.email "" 检查: git config --list 2.退出 exit 新建本地git mkdir

2014-04-17 20:13:38 685

原创 week1

目录操作 getwd(); setwd('../') if(!file.exists(directoryname)){ dir.create(directoryname) } 下载

2014-04-13 10:10:48 878

原创 Week 1

目录 getwd(); setwd('c:/');

2014-04-08 22:15:58 570

原创杂

随机数生成 %产生1到m之间的随机数生成器 rand_indices = randperm(m); %生成100个，返回一个向量 sel = X(rand_indices(1:100),:);

2014-04-07 16:15:56 703

原创 linear regression几个代码

J = 0; grad = zeros(size(theta)); %对于y=1，cost为 -log[h(x)] % y=0，cost为 -log[1-h(x)] %J可以表示为 - log[| if（y==0，1，0） -h(x) |] if(判断，true的输出值，false的输出值) z = X*theta; neg=find(y==0); tmp

2014-03-30 11:10:04 1014

原创牛顿法求解柯西分布参数的最大似然估计（P055）

1.牛顿法求方程根假设f(x)在x=a处取到0，即f(a)=0。已知c在a的附近，有a=c+Δx。根据泰勒公式展开，省略高次项，得 2.简单代码实现 newt_cauchy.R cauchy_mle=function(x,theta){ sum = 0; for(i in 1:length(x)){ sum = sum+(x[i]-theta)

2014-03-27 22:06:33 10043 2