Sunday, February 18, 2018

Lowess smooth in logistic regression using R based on Applied Logistic Regression, 3th ed, HOSMER

To to determine linearity in the logit for a continuous variable, one can use lowess smoothed scatterplot as described in Applied Logistic Regression, 3rd Edition, in chapter 4.2.1 Methods to Examine the Scale of a Continuous Covariate in the Logit.

In the book this is illustrated using Scale_Example data from the book and STATA. I dont have stata, but this can be replicated with a good agreement in R.

For this, I will use logitloess function described here.

The data sets used in the book can be obtained here: # read in Scale_Example.txt that contains data for the example
mydata = read.table("~/Applied_Logistic_Regression_by_Hosmer_Lemeshow_3th_ed_2013/Scale Example/Scale_Example.txt")

# seperate data into x and y

# plot lowess graph

The resulting plot is very similar to the one provided in the book (Figure 4.1). Its not exactly same, but this could be due to differences in the implementation of the logitloess methods between R and STATA. Nevertheless, the conclusions from the plot obtained in R are same as those in the book.

Sunday, February 11, 2018

SublimeText3: Gcc build system and termiantor

create new build system subl3 ~/.config/sublime-text-3/Packages/User/gcc.sublime-build with the contents of { "shell_cmd": "g++ \"${file}\" -o \"${file_path}/${file_base_name}\"", "file_regex": "^(..[^:]*):([0-9]+):?([0-9]+)?:? (.*)$", "working_dir": "${file_path}", "selector": "source.c, source.c++", "variants": [ { "name": "Run", "shell_cmd": "terminator -e 'gcc \"${file}\" -o \"${file_path}/${file_base_name}\" && \"${file_path}/${file_base_name}\";echo;echo; echo Press ENTER to continue; read line;exit; exec bash'" } ] } Based on

Thursday, February 08, 2018

Arch: dell inspiron 1545 Broadcom Limited BCM4312

Have to install this driver:
sudo pacman -S linux-headers broadcom-wl-dkms

Then after reboot wlp12s0 device was found: ip link
1: lo: mtu 65536 qdisc noqueue state UNKNOWN mode DEFAULT group default qlen 1000
link/loopback 00:00:00:00:00:00 brd 00:00:00:00:00:00
2: enp9s0: mtu 1500 qdisc fq_codel state UP mode DEFAULT group default qlen 1000
link/ether 00:25:64:47:bd:bb brd ff:ff:ff:ff:ff:ff
3: wlp12s0: mtu 1500 qdisc fq_codel state UNKNOWN mode DEFAULT group default qlen
link/ether 00:22:5f:d3:14:d2 brd ff:ff:ff:ff:ff:ff

Tuesday, February 06, 2018

Arch: Wine and GPower

GPower 4.1 installs in wine 3.0 and runs, but freezes after executing calculations due to some problems with drawing the results's plot. To fix the freezing, need to install gdiplus and lib32-gnutls:
winetricks vcrun6 lib32-gnutls gdiplus
This fixes the freezes, although the plots can't be still seen.

Friday, February 02, 2018