I am new modeling and I am trying to analyse my data using the package Double
Hierarchical Generalized Linear Models. However, I always get the same sort of
error:
Error in z %*% v_h : non-conformable arguments.
I think that this error indicates that the matrices are created with
non-calculable dimensions, but I don't know how to solve it.
I am creating my model with two fixed and two random effects. with this format:
rm(list=ls(all=TRUE))
library("dhglm")
data=read.csv("log_data.csv", header=TRUE)
phi<-matrix(1,n<-nrow(data),1)
lambda<-matrix(1,n<-nrow(data),1)
data<-cbind(data,phi,lambda)
sr <- as.numeric(data$SR) ##Dependent variable
id <- as.numeric(data$ID) ##individual number
obs <- as.numeric(data$Obs) ##Observation number
o2r <- as.numeric(data$O2R) ##Fixed factor 1 (oxygen consumption 1)
o2a <-as.numeric(data$O2A) ##Fixed factor 2 (oxygen consumption 2)
lambda <-as.numeric(data$lambda)
phi<-as.numeric(data$phi)
model_mu<-DHGLMMODELING(Model="mean",
Link="log",
LinPred= sr~o2a+ o2r + (1|id) + (1|obs),
RandDist= c("inverse-gamma", "gamma"),
LinPredRandVariance=lambda~ 1 +(1|id) + (1|obs) ,
RandDistRandVariance=c("gaussian", "gaussian"))
model_phi<-DHGLMMODELING(Model="dispersion",
Link="log",
LinPred = phi~o2a+ o2r + (1|id) + (1|obs),
RandDist= c("gaussian", "gaussian"))
res_glm<-dhglmfit(RespDist="gaussian",
DataMain=data,
MeanModel=model_mu,
DispersionModel=model_phi,
PhiFix=NULL,
LamFix=NULL,mord=1,dord=1,Maxiter=200,convergence=1e-06)
Mariana Velasque
PhD Student
Marine Biology & Ecology Research Centre | School of Marine Science and
Engineering | Plymouth University | 616 Davy | Drakes Circus | Plymouth | PL4
8AA | UK | +44 1752 5 84598 (Shared phone) |
mariana.velasquebor...@plymouth.ac.uk
<https://webmail.plymouth.ac.uk/owa/redir.aspx?C=0UUp9Pz8dkK1n_4m1deYkaEPZLevftEIjV4BgdiwFBa9Smiq3Nzoxl6fFxk5PwzRXWIRYoEPm8k.&URL=mailto%3amariana.velasqueborges%40plymouth.ac.uk>
| marianavelas...@gmail.com
<https://webmail.plymouth.ac.uk/owa/redir.aspx?C=0UUp9Pz8dkK1n_4m1deYkaEPZLevftEIjV4BgdiwFBa9Smiq3Nzoxl6fFxk5PwzRXWIRYoEPm8k.&URL=mailto%3amarianavelasque%40gmail.com>
Doutoranda
Centro de Pesquisa em Biologia Marinha & Ecologia | Faculdade de Ciências e
Engenharia Marinha | Universidade de Plymouth | 616 Davy | Drakes Circus |
Plymouth | PL4 8AA | UK | +44 1752 5 84598 (Telefone compartilhado)
|mariana.velasquebor...@plymouth.ac.uk
<https://webmail.plymouth.ac.uk/owa/redir.aspx?C=0UUp9Pz8dkK1n_4m1deYkaEPZLevftEIjV4BgdiwFBa9Smiq3Nzoxl6fFxk5PwzRXWIRYoEPm8k.&URL=mailto%3a%7cmariana.velasqueborges%40plymouth.ac.uk>
|marianavelas...@gmail.com
<https://webmail.plymouth.ac.uk/owa/redir.aspx?C=0UUp9Pz8dkK1n_4m1deYkaEPZLevftEIjV4BgdiwFBa9Smiq3Nzoxl6fFxk5PwzRXWIRYoEPm8k.&URL=mailto%3amarianavelasque%40gmail.com>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
