Displaying 1 result from an estimated 1 matches for "sigmab".
Did you mean:
sigma
2009 Feb 02
0
Using Information from the Stats4 package in base envir
...that are defined by a
given cluster label?
Thanks again for your help.
library(stats4)
#Define the Models
#Two simple linear models:
#Model A: respA = a0 + a1x + e(a), e(a) ~N(0, sigma^2(a))
#Model B: respB = b0 + b1x + e(b), e(b) ~N(0, sigma^2(b))
a0= 0; a1 = 1.5; sigmaA=4; b0=50; b1=-2; sigmaB=4; n=250; min=0; max=50;
#Generate Data from the Models
x1 = runif(n, min, max); #Explanatory variable
eA = rnorm(n, 0, sigmaA); eB = rnorm(n, 0, sigmaB);
respA = a0 +a1*x1 + eA; respB = b0 +b1*x1 + eB; #Responses
#K-Means Clustering
resp1 = c(respA, respB);#Put response into single vecto...