[R] How to correct my error message

[varin sacha]

Dear All,

Many thanks for your responses. I got it !

Best,









Le mardi 27 octobre 2020 à 21:16:50 UTC+1, Md. Moyazzem Hossain <hossainmm using 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.goslee using 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 using 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 using 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 using 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