<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Hi Suman,</div><div><br></div><div>about the predRisk function of PredictABEL, we think that </div><div><br class="Apple-interchange-newline"><blockquote type="cite">
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
<div bgcolor="#ffffff" text="#000000">
<font face="Cambria">function (riskModel, data, cID, filename) <br>
{<br>
if (any(class(riskModel) == "glm")) {<br>
predrisk <- predict(riskModel, type = "response")<br></font></div></blockquote><div><br></div>here you meant actually</div><div><span class="Apple-style-span" style="font-family: Cambria; "><br></span></div><div><span class="Apple-style-span" style="font-family: Cambria; ">predrisk <- predict(riskModel, newdata=data, type = "response")</span></div><div><span class="Apple-style-span" style="font-family: Cambria; "><br></span></div><div><span class="Apple-style-span" style="font-family: Cambria; ">which would allow to do prediction with estimated model on new data</span></div><div><br><blockquote type="cite"><div bgcolor="#ffffff" text="#000000"><font face="Cambria">
}<br>
else {<br>
stop("The argument 'riskModel' should be a (GLM)model")<br>
}<br>
if (!missing(data) && !missing(cID) &&
!missing(filename)) {<br>
tab <- cbind(ID = data[, cID], PredRisk = predrisk)<br>
write.table(tab, file = filename, row.names = FALSE, <br>
sep = "\t")<br>
}<br>
return(predrisk)<br>
}<br>
<br></font></div></blockquote><div><br></div>Should we fix that in the code in repo?<br><div><br></div></div>bw,<div>Lina & Yurii</div></body></html>