for (i in 1:38) {
a<-as.data.frame(table(y10_t20[,i]))
names(a)[1]<-"cat_id"
names(a)[2]<-"cat_freq"
a$var.name<-var[i,]
b<-paste("a",i,sep="")
assign(b,as.data.frame(a))
rm(a)
}
for (i in 1:12) {
a<-as.data.frame(table(y10_t30[,i]))
names(a)[1]<-"cat_id"
names(a)[2]<-"cat_freq"
a$var.name<-var[i,]
b<-paste("a",i,sep="")
assign(b,as.data.frame(a))
rm(a)
}
t2<-rbind(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12)
var<-as.data.frame(colnames(y10_t40))
for (i in 1:5) {
a<-as.data.frame(table(y10_t40[,i]))
names(a)[1]<-"cat_id"
names(a)[2]<-"cat_freq"
a$var.name<-var[i,]
b<-paste("a",i,sep="")
assign(b,as.data.frame(a))
rm(a)
}
t3<-rbind(a1,a2,a3,a4,a5)
var<-as.data.frame(colnames(y10_t53))
for (i in 1:8) {
a<-as.data.frame(table(y10_t53[,i]))
names(a)[1]<-"cat_id"
names(a)[2]<-"cat_freq"
a$var.name<-var[i,]
b<-paste("a",i,sep="")
assign(b,as.data.frame(a))
rm(a)
}
library(foreign) library(plyr) library(reshape)
데이터 불러오기
setwd("/home/rserv/rserv-knitr/nps")
y10_t20<-read.dta("y10_t20.dta",convert.dates=T,convert.factors =F,missing.type=F,convert.underscore=F,warn.missing.labels=T) y10_t30<-read.dta("y10_t30.dta",convert.dates=T,convert.factors =F,missing.type=F,convert.underscore=F,warn.missing.labels=T) y10_t40<-read.dta("y10_t40.dta",convert.dates=T,convert.factors =F,missing.type=F,convert.underscore=F,warn.missing.labels=T) y10_t53<-read.dta("y10_t53.dta",convert.dates=T,convert.factors =F,missing.type=F,convert.underscore=F,warn.missing.labels=T) y10_hokey<-read.dta("y10_hokey.dta",convert.dates=T,convert.factors =F,missing.type=F,convert.underscore=F,warn.missing.labels=T)
변수 기본 정보 추출
var<-as.data.frame(colnames(y10_t20))
for (i in 1:38) { a<-as.data.frame(table(y10_t20[,i])) names(a)[1]<-"cat_id" names(a)[2]<-"cat_freq" a$var.name<-var[i,] b<-paste("a",i,sep="") assign(b,as.data.frame(a)) rm(a) }
t1<-rbind(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12,a13,a14,a15,a16,a17,a18,a19,a20,a21,a22,a23,a24,a25,a26,a27,a28,a29,a30,a31,a32,a33,a34,a35,a36,a37,a38)
var<-as.data.frame(colnames(y10_t30))
for (i in 1:12) { a<-as.data.frame(table(y10_t30[,i])) names(a)[1]<-"cat_id" names(a)[2]<-"cat_freq" a$var.name<-var[i,] b<-paste("a",i,sep="") assign(b,as.data.frame(a)) rm(a) }
t2<-rbind(a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11,a12)
var<-as.data.frame(colnames(y10_t40))
for (i in 1:5) { a<-as.data.frame(table(y10_t40[,i])) names(a)[1]<-"cat_id" names(a)[2]<-"cat_freq" a$var.name<-var[i,] b<-paste("a",i,sep="") assign(b,as.data.frame(a)) rm(a) }
t3<-rbind(a1,a2,a3,a4,a5)
var<-as.data.frame(colnames(y10_t53))
for (i in 1:8) { a<-as.data.frame(table(y10_t53[,i])) names(a)[1]<-"cat_id" names(a)[2]<-"cat_freq" a$var.name<-var[i,] b<-paste("a",i,sep="") assign(b,as.data.frame(a)) rm(a) }
t4<-rbind(a1,a2,a3,a4,a5,a6,a7,a8)
t<-rbind(t1,t2,t3,t4) rm(t1,t2,t3,t4)
for (i in 1:38){ rm(a[i]) }
write.dta(t,"var.t20.dta") rm(t)