r - glmmLasso error and warning -


i trying perform variable selection in generalized linear mixed model using glmmlasso, coming error , warning, can not resolve. dataset unbalanced, participants (ptno) having more samples others; no missing data. dependent variable binary, other variables (beside id variable ptno) continous. suspect generic happening, fail see , have not found solution in documentation or on web. code, adapted glmmlasso soccer example is:

glm8 <- glmmlasso(group~ndufv2_dctabl+gper1_dctabl+ esr1_dctabl+esr2_dctabl+klf12_dctabl+sp4_dctabl+sp1_dctabl+  pgam1_dctabl+ank3_dctabl+rasgrp1_dctabl+akt1_dctabl+nudt1_dctabl+                   polg_dctabl+   adarb1_dctabl+ogg_dctabl+ pde4b_dctabl+  gsk3b_dctabl+ apoe_dctabl+  mapk6_dctabl, rnd = list(ptno=~1),       family = poisson(link = log), data = stackdata, lambda=100,      control = list(print.iter=true,start=c(1,rep(0,29)),q.start=0.7))

the error message displayed below. specficially, not believe there nas in dataset , unsure meaning of warning regarding factor variable.

iteration 1 error in grad.lasso[b.is.0] <- score.beta[b.is.0] - lambda.b * sign(score.beta[b.is.0]) : nas not allowed in subscripted assignments in addition: warning message: in ops.factor(y, mu) : ‘-’ not meaningful factors

an abbreviated dataset containing necessary variables available in r format , can downladed here. hope can guided bit how go on analysis. please let me know if there wrong dataset or cannot download it. appreciated.

just follow on @kristofersen comment above. indeed start vector messes analysis up.

if run

glm8 <- glmmlasso(group~ndufv2_dctabl+gper1_dctabl+ esr1_dctabl+esr2_dctabl+klf12_dctabl+sp4_dctabl+sp1_dctabl+  pgam1_dctabl+ank3_dctabl+rasgrp1_dctabl+akt1_dctabl+nudt1_dctabl+                   polg_dctabl+   adarb1_dctabl+ogg_dctabl+ pde4b_dctabl+  gsk3b_dctabl+ apoe_dctabl+  mapk6_dctabl,                    rnd = list(ptno=~1),                    family = binomial(),                    data = stackdata,                    lambda=100,                        control = list(print.iter=true))

then fine , dandy (i.e., converges , produces solution). have copied example poisson regression , need tweak code situation. have no idea whether output makes sense.

quick note: ran binomial distribution in code above since outcome binary. if makes sense estimate relative risks poisson may reasonable (and converges), need recode outcome 2 groups defined 1 , 2 , mess poisson regression.

in other words 

stackdata$group <- stackdata$group-1

before run analysis.


Comments

Post a Comment

Popular posts from this blog

php - How to add and update images or image url in Volusion using Volusion API -

actually have store importing products store product data insert , update image or image url not insert , update. this xml code name datapro.txt. <products_joined> <productcode>3710_012t</productcode> <vendor_partno>eah5450silent/di/1gd3(lp)</vendor_partno> <productname>test product ta</productname> <hideproduct>n</hideproduct> <stockstatus>20</stockstatus> <lastmodified>1/5/2016 10:25:00 am</lastmodified> <lastmodby>2</lastmodby> <productweight>0.9</productweight> <productprice>100</productprice> <productmanufacturer>asus tek</productmanufacturer> <vendor_price>32.69</vendor_price> <numproductssharingstock>0</numproductssharingstock> <categoryids>107</categoryids> <producturl>http://tebkq.mvlce.servertrust.com/productdetails.asp?productcode=3710_012t</producturl>
Read more

javascript - IE9 error '$'is not defined -

Image
<!--[if lte ie 9]><html class="lte9"><![endif]--><!--[if gt ie 9]><html><![endif]--> <!--[if !ie]><!--> <html> <!--<![endif]--> blah blah <script src="https://ajax.googleapis.com/ajax/libs/jquery/1.12.1/jquery.min.js"></script> <script src="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.7/js/bootstrap.min.js"></script> <script src="{{ asset('js/mypage.js') }}" type="text/javascript"></script> </body> </html> when load page, ie9 brings error '$' not defined . how can solve error? +++ mypage.js $(function(){ //when load page, //ie9 tag not moving. $(window).load(function() { $('div.mask').hide(); }); // when click , willshow show. //for not ie 9 if (!$('html').hasclass('lte9')) { var = $('ul.about > li > span.activebtn');
Read more