动态代码.doc

蚄#ImporttheStrangeDataset葿rawdata=readLines("")肇#Gettheusefulwords螆data=unlist(strsplit(rawdata,"\\[\\{\"|\":\"|\",\"|\":\\[\\[|\\],\\[|\\]\\],\"|\\]\\]\\},\\{\"|\\]\\]\\}\\]"))螁data=data[which(data!="")]膀n=length(data)袅#Setindex袆nameind=which(data=="name")膁regind=which(data=="region")蚈incind=which(data=="e")袈popind=which(data=="population")羆lifind=which(data=="lifeExpectancy")薂endpoint=nameind-1莀endpoint=endpoint[2:180]蚇endpoint[180]=45999肆name=data[nameind+1]羃region=data[regind+1]螈#=(name=rep(0,210),region=rep(0,210),year=1800:2009,e=rep(0,210),pop=rep(0,210),life=rep(0,210))膅#ma莄splt=function(x)return((unlist(strsplit(x,","))))薀#Linearinterpolationfunction,forthosezerodata葿itpl=function(a)芅{薁 ind=which(a>0)节 if(ind[1]>1)膈 a[1:(ind[1]-1)]=rep(a[ind[1]],ind[1]-1)莅 n=length(ind)羂 if(ind[n]<length(a))蚀 a[(ind[n]+1):length(a)]=rep(a[ind[n]],length(a)-ind[n])羇 for(iin1:(n-1))莅 a[ind[i]:ind[i+1]]=rep(a[ind[i]],ind[i+1]-ind[i]+1)+(seq(ind[i],ind[i+1],1)-ind[i])*(a[ind[i+1]]-a[ind[i]])/(ind[i+1]-ind[i])薁 return(a)虿}羃#Dragdatainformationfromwords蚃tbls=NULL羁for(iin1:180)肇{羆 dfr=initdfr螃 dfr$name=rep(name[i],210)肈 dfr$region=rep(region[i],210)蝿 for(jin(incind[i]+1):(popind[i]-1))螅 {袃 tmp=splt(data[j])葿 ii=which(dfr$year==tmp[1])芇 dfr$e[ii]=tmp[2]薄 }羂 for(jin(popind[i]+1):(lifind[i]-1))袀 {罿 tmp=splt(data[j])薇 ii=which(df