Dear All, Many thanks for your responses. I got it !
Best, Le mardi 27 octobre 2020 à 21:16:50 UTC+1, Md. Moyazzem Hossain <hossai...@juniv.edu> a écrit : Dear Varin, I think the following code will solve your problem. n <- 60 b <- runif(n, 0, 5) a <- runif(n, 0, 5) z1 <- data.frame(x0=1:57, x1=rnorm(n*0.95,2,3)) z2 <- data.frame(x0=58:60, x1=rnorm(n*0.05,2,9)) combined=rbind(z1,z2) z=combined[,2] y_model <- 0.1 * b - 0.5 * z - a + 10 y_obs <- y_model +c( rnorm(n*0.95, 0, 0.1), rnorm(n*0.05, 0, 0.5) ) df<-data.frame(b,a,z,y_obs) Thanks. Md On Tue, Oct 27, 2020 at 7:21 PM Sarah Goslee <sarah.gos...@gmail.com> wrote: > Hi, > > a is of length 60. > b is of length 60. > z is of length 57. > > What do you expect to have happen when you create y_model ? What > happens to those other 3 observations? > > Sarah > > On Tue, Oct 27, 2020 at 3:07 PM varin sacha via R-help > <r-help@r-project.org> wrote: >> >> Dear R-experts, >> >> Here below my R code. The warning message is not a problem to me but there >> is an error message more problematic. I understand the error message but I >> don't know if it is possible to correct the error and if yes, how to correct >> it. >> >> Many thanks. >> >> >> n <- 60 >> b <- runif(n, 0, 5) >> a <- runif(n, 0, 5) >> z <- rnorm(n*0.95,2,3) + rnorm(n*0.05,2,9) >> y_model <- 0.1 * b - 0.5 * z - a + 10 >> y_obs <- y_model +c( rnorm(n*0.95, 0, 0.1), rnorm(n*0.05, 0, 0.5) ) >> df<-data.frame(b,a,z,y_obs) >> >> ______________________________________________ >> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > > > -- > Sarah Goslee (she/her) > http://www.numberwright.com > > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > -- Best Regards, Md. Moyazzem Hossain Associate Professor Department of Statistics Jahangirnagar University Savar, Dhaka-1342 Bangladesh Website: http://www.juniv.edu/teachers/hossainmm Research: Google Scholar; ResearchGate; ORCID iD ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.