- 博客(22)
- 收藏
- 关注
RSS订阅原创 library(reshape) #变量名修改
library(reshape) #变量名修改newdata <- rename(newdata,c(ku2 = "ifinternet"))table(newdata$ifinternet)
2022-05-29 11:31:23
117
原创 stata读取csv数据，变量名不对，数据无法读取
//help psmatch2//帮助文件//ssc install psmatch2//安装clearimport excel D:\2、R+Stata\statawork\all_sports_2018.xlsx, firstrow ***——————————————————————————————————————————**** clearinsheet using D:\2、R+Stata\statawork\all_sports_2018.csv,clear namedesc.
2022-05-04 09:41:13
2088
原创 漫画app分析
//help psmatch2//帮助文件/*——————————————————————————————————————————————————————————————————————————————————————————————————————————————* clearinsheet using C:\Users\Administrator\Desktop\本科生论文\2、漫画可持续使用\manhua_screen.csvsave C:\Users\Administrator\D.
2022-04-21 19:19:21
1347
原创 如何将r语言结果 输出结果到WORD
##输出结果到WORD!library(stargazer) ##生成结果用word打开!stargazer( fit2, fit3, fit4, fit5, title = "results", align = F, type = "text", no.space = TRUE, out = "fit.html")
2022-03-06 12:26:00
3337
原创 视频分析代码douyin
## clear data & close graphs--------rm(list=ls())graphics.off()#install.packages("dplyr")#install.packages("pastecs")#install.packages("haven")library("dplyr")library("pastecs")#设置当前工作空间-----------getwd()#setwd("D:/ⅡR+Stata/R语言代码2021/R语言包 C.
2022-03-04 08:35:12
345
原创 计算kappa值
data0<-read.csv("C:\\Users\\Administrator\\Desktop\\xl1.csv")dim(data0)data0[,1]bianma1<-data0[1:30,2:12]bianma2<-data0[31:60,2:12]bianma1colnames(bianma1)vec.bianma1 <- c(as.matrix(bianma1))length(vec.bianma1)vec.bianma2 <- c(as..
2022-02-26 18:00:11
566
原创 频率分类llll
##2.游戏频率编码)table(data$D2b_6)data$notuse <- as.numeric(data$D2b_6==0) data$rarely <- as.numeric(data$D2b_6==1) #很少data$occasional <- as.numeric(data$D2b_6>=2&data$D2b_6<=3) #经常data$usually <- as.numeric(data$D2b_6>=4&dat.
2022-01-21 20:49:58
1850
原创 核心价值观 分组
##核心价值观（12组）name.hxjzg<-paste0('g7_',c(2,3,5,6,8,9,11,12,14,15,17,18)) #12个选项在data中的名称data.hxjzg<-as.matrix(data[,name.hxjzg]) #12个选项单独存为一个数据，便于后面的NA处理dim(data.hxjzg) #查看维度data.omit.NA<-gtools::na.replace(data.hxjzg,0) #矩阵中的缺失值用0代替sum.hxjzg.
2022-01-21 14:42:12
237
原创 stata 均值和标准差
*（描述统计1，(与R语言一致)self_hly ifinternet gender age huji edu lnincome_fam a36 v458 class east west tabstat self_hly ,stats(mean n q max min sd var cv) *（描述统计2 ，summarize命令可以得出所有处理过的变量的平均值summarize...
2022-01-19 08:58:42
15784
原创 cfps数据库
##(9)创建区域(东中西)data<-data[data$provcd20>=0,] #剔除选项为负数的(1831),注:ear201a不行(负数太多)dim(data)table(data$provcd20)data$east=rep(NA,nrow(data)) #东部地区data$west=rep(NA,nrow(data)) #西部地区data$central=rep(NA,nrow(data)) #中部地区for(i in 1:nrow(data)){ if(d..
2022-01-04 23:16:56
3696
原创 stata 各种回归模型
ologit a36 whxf sex age hunyin a15 a17 edu a35 lnincome_fam east westlogit a36 whxf sex age hunyin a15 a17 edu a35 lnincome_fam east westprobit a36 whxf sex age hunyin a15 a17 edu a35 lnincome_fam east westoprobit a36 whx...
2021-12-30 10:54:35
3128
原创 广义线性回归 logit和probit
## 广义线性回归___Logit模型(y为体育锻炼强度)。#install.packages("ResourceSelection") library(ResourceSelection) #Hosmer-Lemeshow拟合优度检fit4<- glm(fitness2~hlw,data = newdata, binomial(link="logit"))summary(fit4) hl1 <- hoslem.test(fit4$y,fitted(fit4),g=10) ##.
2021-12-29 13:00:10
566
原创 R语言ATT
#五、倾向得分匹配&ATT##1-1平衡性检验#install.packages("MatchIt")library("MatchIt")for(i in 1:nrow(newdata)){ if(newdata$a285[i]==1){newdata$ifinternet[i]=0} else {newdata$ifinternet[i]=1}}table(newdata$ifinternet)#1-2基本处理##以互联网为例，对新媒介进行匹配；#nearest加了r.
2021-12-29 12:55:31
914
原创 stata稳健性、内生性、logit模型
//help psmatch2//帮助文件//ssc install psmatch2//安装clearinsheet using D:\ⅡR+Stata\statawork\mergeelderly.csvsave D:\ⅡR+Stata\statawork\mergeelderly.dta use D:\ⅡR+Stata\statawork\mergeelderly.dta ** 将数据随机排序set seed 10101gen ranorder = runiform().
2021-12-29 12:50:33
13599
1
原创 corrplot画图
##(六)描述统计(画图)library("corrplot")library("dplyr")##(1)画图：混合方法之上三角为圆形，下三角为数字coroutput <-cor(newdata[,c("a15","a285","gender","age","hunyin","lnincome_fam","edu","a35","v458","a36")])coroutputcolnames(coroutput)<-c("自评健康","互联网","性别","年龄","婚姻","收.
2021-12-29 12:48:58
263
原创 根据p值计算显著性
tvalue<-c(5.83, 4.22, 5.83, 4.24, 5.83, 5.12, 5.83, 6.86)sig<-ifelse(tvalue>=2.56,"显著***", ifelse(tvalue>=1.96, "显著**", ifelse(tvalue>=1.64,"显著*","不显著")))sig...
2021-12-29 12:08:03
2047
原创 随机森林与变量重要性检测
###随机森林筛选变量重要性排序!#(1)训练集、测试#将数据集分为训练集和测试集,比例为7:3library(openxlsx)train_sub = sample(nrow(newdata),7/10*nrow(newdata))train_data = newdata[train_sub,]test_data = newdata[-train_sub,]#(2)step_2library(pROC) #绘制ROC曲线library(randomForest)#数据预处理trai.
2021-12-29 11:57:10
2807
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人