- 博客(15)
- 资源 (6)
- 收藏
- 关注
RSS订阅原创 R语言学习笔记——数据转换为正态分布
library(AID)out <- boxcoxnc(data, method = "mle", lambda = seq(-2,2,0.0001), verbose = F, plot = F)out$lambda.hat## [1] -0.0474 library(MASS)out <- boxcox(data~1, lambda = seq(-2,2,0.0001), plotit = F)out$x[which.max(out$y)]## [1] -0.0474 l
2022-05-04 12:01:33
3579
原创 R语言学习笔记——“东北乱炖”
######缺失值#########install.packages("Amelia")#install.packages("mice")#install.packages("VIM")#install.packages("DMwR2")library("Amelia")library("mice")library("VIM")library("DMwR2")#simulated exampleset.seed(123)n.mis<-10x1<-seq(1.51,4,0.
2021-11-25 23:07:12
423
原创 R语言学习笔记——空间自相关：全局Moran’I/局部Moran’I /Geary’c/Moran散点图
data:ChinaRD2<-readRDS("gadm36_CHN_1_sp.rds")# download the file: https://gadm.org/download_country_v3.htmlnames(ChinaRD2)[4]<-c("province_name")popDataSP<-merge(ChinaRD2, data_sp, by=c("province_name"),all.x=TRUE)popDataSP<-popDataSP[!is
2021-03-08 16:52:06
5219
3
原创 R语言学习笔记——估计概率密度函数
估计概率密度函数——R语言学习笔记参数法正态分布方法1:fit.normal <- fitdist(data_tmp$VE,"norm")summary(fit.normal)方法2:fit.normal_b<-Bnormal(data_tmp$VE, priors=list(muMean=13.651040, muSD=3.998584))summary(fit.normal_b)伽马分布方法1:fit.gamma <- fitdist(data_tmp$VA,
2021-03-08 16:00:25
8271
原创 R语言学习笔记——项目里的一些小函数
统计缺失个数na_count<-function(data,x){sum(data[,x]==""|is.na(data[,x]))}eg:na_count(Data_Tmp,“hypertension_before”)插补centralImputation<-function (data,xx){for (i in c(xx)) {if (any(idx <- which(data[, i]==""|is.na(data[, i]))))data[idx,i]<
2021-01-03 15:47:26
611
原创 R语言学习笔记——中国省/市热力图
require(maps)require(ggmap)library(mapdata)library(maps)library(ggplot2)library(mapdata) library(autoimage) library(RColorBrewer)library(ggplot2)library(ggthemes)# devtools::install_github("sjmgarnier/viridis") #这个好像没用到library(viridis) #devtool
2020-06-28 06:48:36
5283
3
原创 R语言学习笔记——数据拆分
有一种分隔符data5[,1]<-as.character(data5[,1])month<-t(as.data.frame(strsplit(data5[,1],split="/")))[,2]有两种分隔符a<-strsplit(data5[,1],split="/|-")#数据中有/和-两种分隔符huoguo<-array()xiangguo<-arr...
2020-04-05 11:51:21
3980
原创 R语言学习笔记——去除重复的行/计算重复次数
去除重复的行data <- data.frame(id1 = c(1, 1, 1, 2, 2), id2 = c(2, 2, 3, 2, 2))library(dplyr)data %>% group_by(id1, id2) %>% filter(row_number() == 1) %>% ungroup()...
2020-04-05 11:27:37
7522
1
原创 R语言学习笔记——GBK和UTF-8转换
背景介绍使用Mac打开用Windows系统写的script 和数据时会出错，主要是中文显示为乱码1、R script打开 r scriptfile- reopen with encoding- GB18030重新打开后中文显示正常2、数据读取csv文件data1<-read.csv(“data_normal.csv”,fileEncoding = ‘GBK’)其中data_...
2020-04-05 11:19:13
6341
原创 R语言学习笔记——添加ID
patient<-data.frame(data$时间,data$医院名称,data$患者姓名,data$孕前体重.kg.,data$身高.cm.,data$年龄.岁.)###根据数据中的检查时间、医院名称、患者姓名、体重、身高来判断是否是同一个人colnames(patient)<-c("time","hospital","name","weight","height","age"......
2019-11-20 04:33:04
2540
原创 R语言学习笔记——合并两个公式
form1<-y~x1+x2form2<-y~x3+x4merge.formula <- function(form1, form2, ...){ # get character strings of the names for the responses # (i.e. left hand sides, lhs) lhs1 <- deparse(form1[...
2019-10-15 06:34:49
890
原创 R语言学习笔记——数据处理
读入数据：data1<-read.csv("city_hospital.csv")data2<-read.csv("data_province_old.csv")#data1和data2共同的列名：医院名称#给data2的医院添加医院等级，医院等级存在data1中##left模式匹配data3<-merge(data2,data1,all.x=TRUE,sort=TR...
2019-09-30 23:17:21
495
原创 R语言学习笔记——在图形下面添加表格
CV_1 <- 0.2CV_2 <- 0.3Mean <- 65sigma_1 <- sqrt(log(1 + CV_1^2))mu_1 <- log(Mean) - sigma_1^2 / 2sigma_2 <- sqrt(log(1 + CV_2^2))mu_2 <- log(Mean) - sigma_2^2 / 2q <- c(...
2019-09-24 00:47:56
2893
原创 R学习笔记-可视化-带有置信区间的条形图
带有置信区间的条形图library(gplots)attach(airquality)heights <- tapply(Temp,Month,mean)lower <- tapply(Temp,Month,function(v) t.test(v)$conf.int[1])upper <- tapply(Temp,Month,function(v) t.test(v)...
2019-09-13 23:16:47
8182
3
原创 缺失值处理方法——基于R语言
缺失值处理方法——基于R语言缺失值处理方法有很多，除了删除、利用平均值或中位数等插补外，本文简要介绍一些基于模型的插补方法。判别如果因变量可以由若干个非缺失变量解释的话，选择非缺失变量为自变量，含缺失值的变量为因变量，建立合适的判别模型进行判别。热卡基于GIBBS抽样KNN（未完待续）如何插入一段漂亮的代码片去博客设置页面，选择一款你喜欢的代码片高亮样式，下面展示同样高亮的 代码...
2019-09-12 09:54:38
2535
Intellectual property rights in innovation management research-A review.pdf
2020-04-06
Correction for population stratification in random forest analysis
2019-10-10
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人