Displaying 1 result from an estimated 1 matches for "v_tild".
Did you mean:
v_tilde
2010 Mar 09
1
error in a function
...det_S=det(R)
for (i in 1:t){
b=i*n
y_ausgew=y_tilde[a:b]#y_ausgew ist numeric
y_aus=matrix(y_ausgew)#Typen anpassen
X_ausgew=X[a:b,]#X_ausgew ist eine Matrix
x_tilde=X_ausgew-1/t*x_t #x_tilde ist auch eine Matrix
####X_aus=matrix(X_ausgew)
V1=R%*%(S%*%y_aus-x_tilde%*%beta)
V11=t(V1)
V2=V11%*%J_n
V_tilde=V2%*%V1
sumterm=sumterm+V_tilde
a=b+1
}
-((n*t)/2)*log(2*pi*sigm_sq)+t*(log(det_S)+log(det_R))-1/(2*sigm_sq)*sumterm
}
where beta is a matrix (2x1), X is a matrix (n*tx1) and y is numeric with
length nt.
My error message is: X_ausgew%*%beta - non-conformable arguments
Can anybody help me?
THA...