[R] Errore : requires numeric/complex matrix/vector arguments
John Kane
jrkrideau at inbox.com
Sun Oct 20 13:20:31 CEST 2013
It looks like you have posted in HTML and your code is basically unreadable.
Please repost in text format. Also I'd suggest reading one or both of the following lists for some hints on how to create a good question.
I think a key point here is to supply sample data using the dput() function as described in the links:
https://github.com/hadley/devtools/wiki/Reproducibility
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
Welcome to the R-help list
John Kane
Kingston ON Canada
> -----Original Message-----
> From: valentina81c at hotmail.it
> Sent: Sun, 20 Oct 2013 11:08:08 +0000
> To: r-help at r-project.org
> Subject: [R] Errore : requires numeric/complex matrix/vector arguments
>
> Dear R users,I'm a new user of R. I'm trying to do a LM test an there is
> this type of error: Error in t(mX) %*% mX : requires numeric/complex
> matrix/vector arguments.
> To be clear I write down the code in which mY ( 126,1 ) mX (126,1)
> mZ(126,1) are matrix.
>
> LMTEST <- function(mY, mX, mZ)#mY, mX, mZ must be matrices!#returns the
> LM test statistic and the degree of freedom{iT = dim(mY)[1]ip =
> dim(mY)[2]iDF = dim(mZ)[2]*ipmE = mY -
> mX%*%solve(t(mX)%*%mX)%*%t(mX)%*%mY
> the error starts from the above step (t(mX)%*%mX)%*%t(mX)%*%mY
> RSS0 = t(mE)%*%mEmXX = cbind(mX, mZ)mK = mE -
> mXX%*%solve(t(mXX)%*%mXX)%*%t(mXX)%*%mERSS1 = t(mK)%*%mKdTR =
> sum(diag(solve(RSS0)%*%RSS1))LM = iT*(ip-dTR)pval =
> 1-pchisq(LM,df=iDF)return( c(pval, LM, iDF) )}
> Any suggestion? Where is the problem? I am getting craxy!
> Valentina
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> 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.
____________________________________________________________
FREE 3D MARINE AQUARIUM SCREENSAVER - Watch dolphins, sharks & orcas on your desktop!
More information about the R-help
mailing list