[18853] in s-news-athena

home help back first fref pref prev next nref lref last post

[S] Query: warning messages

daemon@ATHENA.MIT.EDU (SOFIA Stephano)
Wed Aug 4 13:26:57 1999

From: cs0sso@cis.sunderland.ac.uk (SOFIA Stephano)
To: s-news@wubios.wustl.edu (S-News)
Message-Id: <1999Aug04.182142.1814.2895216@missgate.sunderland.ac.uk>
Mime-Version: 1.0
Content-Type: text/plain; charset="US-ASCII"
Date: Wed, 04 Aug 1999 18:24:05 +0100



Dear S-Plus users,
let me ask you two questions.
1.
Suppose that

G is a square matrix of dimension (n,n);
FTOT is a vector of length n;
V is a number.

I wrote a procedure where at each step ( starting from 2 to N, the length of 
the timeseries Y, and giving the initial step where it is required )
where, according to the calculations below reported,

W and COV are  arrays of matrices, each matrix a square matrix of simension 
(n,n);
m is an array of vectors;
R is a square matrix of dimension (n,n);
Q is a matrix of dimension (1,n);
A,a,e and f are vectors of dimension n.

W[k,  ,  ] <- (1/d - 1) * G %*% COV[k - 1,  ,  ] %*% t(G)
R[k] <- G %*% COV[k - 1,  ,  ] %*% t(G) +W[k,  ,  ]
Q[k] <- drop(t(FTOT) * R[k] * FTOT) + V
A[k] <- (R[k] * FTOT * 1)/Q[k]
a[k] <- G * m[k - 1]
e[k] <- Y[k] - t(FTOT) * a[k]
m[k, , ] <- a[k] + A[k] * e[k]
f[k] <- drop(t(FTOT) * a[k])
COV[k,  ,  ] <- R[k] - A[k] * Q[k] * t(A[k])

At the end of the calculations, I get a warning message for each of the 
previous components and at each step:
" Replacement length not a multiple of number of elements to replace "

What does it mean ? How can I solve these warning messages?

2.
How can I view all the calculations defined inside a function written by 
myself ( not only the final result, but all the steps which brought to this 
result ) ?

thanks for your help an support
Stefano

-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu.  To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message:  unsubscribe s-news

home help back first fref pref prev next nref lref last post