This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

library(ggplot2)
library(kml)
## Loading required package: clv
## Loading required package: cluster
## Loading required package: class
## Loading required package: longitudinalData
## Loading required package: rgl
## Loading required package: misc3d
library(kml3d)
## 
## Attaching package: 'kml3d'
## The following object is masked from 'package:rgl':
## 
##     plot3d
library(kmlShape)
## Loading required package: lattice
## 
## Attaching package: 'kmlShape'
## The following object is masked from 'package:longitudinalData':
## 
##     distFrechet
library(cluster)
library(caret)

set.seed(1)
simulasi<-rbind()
n<-1000
for(i in 1:n){
  #panel data structure set
  d.panel<-data.frame()
  d.panel<-data.frame(cluster=rep(1:3,each=100),cross=rep(1:10,each=10),year=rep(1:10,30))
  X<-as.vector(d.panel[,3])
  e=as.vector(rnorm(300,sd=1))
  alpha=as.vector(c(10+rnorm(100,sd=2),12+rnorm(100,sd=2),14+rnorm(100,sd=2)))
  beta1=as.vector(c(2+rnorm(100,sd=1),2+rnorm(100,sd=1),2+rnorm(100,sd=1)))
  y=alpha+X*beta1+e
  d.panel<-data.frame(d.panel,y)
  
  d.panel<-data.frame(d.panel)
  d.panelt1 <- data.frame(d.panel[d.panel$year==1,])
  d.panelt2 <- data.frame(d.panel[d.panel$year==2,])
  d.panelt3 <- data.frame(d.panel[d.panel$year==3,])
  d.panelt4 <- data.frame(d.panel[d.panel$year==4,])
  d.panelt5 <- data.frame(d.panel[d.panel$year==5,])
  d.panelt6 <- data.frame(d.panel[d.panel$year==6,])
  d.panelt7 <- data.frame(d.panel[d.panel$year==7,])
  d.panelt8 <- data.frame(d.panel[d.panel$year==8,])
  d.panelt9 <- data.frame(d.panel[d.panel$year==9,])
  d.panelt10 <- data.frame(d.panel[d.panel$year==10,])
  oke<-data.frame(d.panelt1[,4],d.panelt2[,4],d.panelt3[,4],d.panelt4[,4],d.panelt5[,4],
                  d.panelt6[,4],d.panelt7[,4],d.panelt8[,4],d.panelt9[,4],d.panelt10[,4])
  colnames(oke)<-c("y11","y12","y13","y14","y15","y16","y17","y18","y19","y110")
}
oke
##          y11      y12       y13       y14      y15      y16      y17      y18
## 1  14.715174 17.92869 25.394384 23.916130 12.47675 29.66773 18.36732 18.49975
## 2  13.498882 10.80070 17.538744 15.712540 23.78111 21.92593 21.30315 19.40513
## 3  12.488915 10.66131 11.282674 21.305503 18.30084 28.82932 22.06497 25.54742
## 4  14.075043 14.75516 19.392112 15.592272 17.40805 28.90614 30.95417 28.11871
## 5   8.548412 17.04356  9.554789  8.736580 10.79209 30.24055 26.28645 36.04729
## 6   7.389169 16.01228 22.092965 20.735905 26.58861 16.02431 27.08116 26.42683
## 7  11.051012 12.64751 14.823045  9.973699 14.70083 13.56023 23.08052 12.19427
## 8  11.719933 14.05639 14.198547 20.120875 26.16792 18.43215 12.78164 25.62624
## 9  11.519021 11.46891 20.428358 17.329255 27.22948 16.30828 26.67708  9.45401
## 10 15.424577 16.01499 22.983456 18.650149 20.77035 22.07744 25.32504 31.67304
## 11 13.516461 21.43079 20.533969 25.221393 28.07269 38.78115 21.30997 39.07829
## 12 14.479105 19.24472 15.026730 20.508139 12.90112 24.89585 31.42319 25.44748
## 13 10.893892 13.11920 20.514423 22.895480 22.64957 35.73405 20.23020 19.72730
## 14  9.460101 17.41591 20.036879 15.422936 23.25562 31.00360 34.81692 34.10495
## 15 10.248706 16.99258 25.680529 13.872659 22.02656 29.47715 33.53741 16.04388
## 16 14.559235 16.08600 15.616710 21.297029 27.71840 20.62518 20.32682 33.07963
## 17 13.941569 13.96556 22.107322 14.384520 20.31275 32.03615 26.51898 29.55617
## 18 15.152581 15.73721 17.729138 20.728245 22.20652 24.95933 29.02245 23.09111
## 19 15.004644 13.44691 20.422730 21.313838 24.39320 25.62924 33.78503 40.89054
## 20 16.039961 13.49952 12.764131 18.340682 24.23891 30.05976 34.40676 33.88523
## 21 19.585809 22.44682 26.088028 22.563872 18.89469 40.99471 23.80509 43.44125
## 22 15.729545 22.81671 13.790146 22.604722 19.68797 23.44624 34.69984 25.13182
## 23 17.501621 19.93778 13.724458 19.991605 16.49664 32.20967 21.51102 32.53524
## 24 14.210130 18.64602 21.551967 23.873884 28.88993 19.07620 35.51937 45.46708
## 25 18.568086 16.31044 20.326163 18.403762 20.32967 36.34886 26.65905 25.00599
## 26 19.118316 18.08538 25.163586 23.500991 17.16370 30.41060 33.11044 14.87724
## 27 17.522668 18.56586 22.463968 19.081496 18.93491 21.10958 24.55480 28.68891
## 28 19.742847 19.23580 25.519549 24.313738 32.53539 15.39508 40.03380 27.20394
## 29 19.792783 22.95068 18.804485 19.917914 21.48793 23.11183 27.77735 33.32246
## 30 17.893984 18.30617 22.899610 13.874222 20.11236 22.19428 25.36364 37.55795
##         y19     y110
## 1  31.19634 28.52012
## 2  40.48840 39.07958
## 3  13.87049 39.15557
## 4  33.96012 44.81182
## 5  24.24370 25.67025
## 6  12.31420 26.86798
## 7  14.22048 35.82603
## 8  31.76034 27.10476
## 9  15.65306 28.18984
## 10 26.69949 28.50302
## 11 26.84726 37.37487
## 12 25.85217 23.16298
## 13 21.25248 40.78630
## 14 26.48077 25.89678
## 15 36.34826 18.15461
## 16 24.42247 39.00619
## 17 34.66268 28.83740
## 18 38.10208 40.29096
## 19 33.27764 35.98121
## 20 24.05786 24.79838
## 21 29.38939 38.34083
## 22 17.63825 28.45934
## 23 31.35596 49.04388
## 24 28.33440 35.04350
## 25 25.61988 45.58422
## 26 25.39116 32.01491
## 27 42.98126 29.72399
## 28 43.06032 11.81876
## 29 24.41738 25.18024
## 30 19.58624 25.08810

define base for the graphs and store in object 'p'

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.