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
<https://scholar.google.com/citations?user=-U03XCgAAAAJ&hl=en&oi=ao>*;
*ResearchGate
<https://www.researchgate.net/profile/Md_Hossain107>*; *ORCID iD
<https://orcid.org/0000-0003-3593-6936>*
[[alternative HTML version deleted]]
______________________________________________
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.
