Displaying 2 results from an estimated 2 matches for "basisvals".
2012 Aug 24
0
diagonal matrix, array attributes and how to keep from setting an attribute on "NULL"
...mnames(Y)",list(rownames(Y),varnames,colnames(Y)))
Y2<-as.array(Y,dimnames=dimnames(Y))
assign("Y",Y2)
then rerun the analysis with the same result
Now traceback() gives:
7: as.array.default(Y)
6: as.array(Y)
5: kronecker(X, Y)
4: kronecker(X, Y)
3: diag(rep(1, nrep)) %x% basisvals$bvals.obs
2: LS.setup(pars, coefs, fn, basisvals, lambda, fd.obj, more, data,
weights, times, quadrature, eps = 1e-06, posproc, poslik,
discrete, names, sparse)
1: Profile.LS(fhn, data = data2, times = times, pars = pars, coefs = coefs,
lambda = lambda, out.meth = "nls&quo...
2012 Aug 08
1
dimnames in array
...seems to me that the analysis will run if I can make this array
have the correct length dimnames
I thank you for your help.
____________
*
fd.data11a = data11a
DEfd11a = Data2fd(fd.data11a, times, bbasis,
fdnames=list(NULL,NULL,varnames))
profile11a.obj=LS.setup(pars,DEfd11a$coefs,fhn,basisvals=bbasis,lambda,fd.obj=NULL,
more=NULL,data=data11a,weights=NULL,times=times,quadrature=NULL,eps=1e-6,
posproc=FALSE,poslik=FALSE,discrete=FALSE,names=NULL,sparse=FALSE)
lik = profile11a.obj$lik
proc = profile11a.obj$proc
qwts = rep(1/length(knots),length(knots))
qwts = qwts %*% t(lambda)...