similar to: "Longitudinal" with binary covariates and outcome

Hi again! Following up my previous posting below (to which no response as yet), I have located a report which situates this type of question in a longitudinal modelling context. http://www4.stat.ncsu.edu/~dzhang2/paper/glm.ps Generalized Linear Models with Longitudinal Covariates Daowen Zhang & Xihong Lin (This work seems to originally date from around 1999). They consider an outcome Y,

Change the assembly code to use only relative references of symbols for the kernel to be PIE compatible. Position Independent Executable (PIE) support will allow to extended the KASLR randomization range below the -2G memory limit. Signed-off-by: Thomas Garnier <thgarnie at google.com> --- arch/x86/crypto/aes-x86_64-asm_64.S | 45 ++++++++----- arch/x86/crypto/aesni-intel_asm.S

Hi Folks, This tip is probably lurking somewhere already, but I've just discovered it the hard way, so it is probably worth passing on for the benefit of those who might otherwise hack their way along the same path. Say (for example) you want to do a logistic regression of a binary response Y on variables X1, X2, X3, X4: GLM <- glm(Y ~ X1 + X2 + X3 + X4) Say there are 1000 cases in the

Dear R-users, actually i try to parse some state protocols for my work. i an easy stetting the code below works fine, if states are reached only once. in harder settings it could be possible that one state gets visited more times. in this case for me its interesting to see how much waiting time lies between to states on the whole. by the way i didn't use R as a parsing tool so far, so

Hello. I'd like to compute the associate vector of distances between leaves in a binary non-rooted tree. The definition of a distance between two leaves in a binary non-rooted tree is the number of edges in the path joining the two leaves. I've tried the ape package but I'm unable to find this vector. For example, using rtree(5,rooted=F) I've obtained the following tree: $edge

Dear All user, Hello, I'm a student and I have some trouble with the experimental (columns-experiments) design of my project. I use a randomized block design with 4 treatments including a control. For each treatment, I use 3 replicates and 3 blocks. The treatments are: -T1 = COD (300 mg/Lit) COD=chemical oxygen demand -T2 = COD (200 mg/Lit) -T3 = COD (100 mg/Lit) -T4 = COD (0 mg/Lit) as

Dear List: I'm trying to use the boot function to estimate some standard errors. I actually programmed a bootstrap using some homebrew code and it worked fine. But, I am trying to use the more efficient boot function. I have placed some sample data for replication of my problem at the bottom of this email. For the sample problem, I have 10 subjects each with 5 observations Y_t = (t_1, ...,

Hi, I have data of the form: t x y trip t1 x1+e y1+e A t2 x2+e y2+e A t3 x3+e y3+e B t4 x4+e y4+e B t5 x5+e y5+e B ... ... ... ... where t is time and x/y are positions in space and e is a random error term. Trips A/B/C/etc are entirely independent

I am currently using the relaimpo package to estimate the relative importance of regressors (N= 4000): > m1 <- lm(y ~ x1+x2+x3+x4+x5+, data=data) > calc.relimp(m1, rela=TRUE) > m2=boot.relimp(m1, boot = 500, rela=TRUE, type="lmg") > booteval.relimp(m2) > plot(booteval.relimp(m2)) In a new dataset with 3 measurement points (0,6,12 weeks), I want to perform a similar

Database structure as follows: t1 -> t2 -> t3 -> t4 select * from t1; +----+-------+-------+ | id | t2_id | t3_id | +----+-------+-------+ | 1 | 1 | null | | 2 | null | 1 | | 3 | null | 2 | +----+-------+-------+ select * from t2; +----+-------+ | id | title | +----+-------+ | 1 | "hat" | +----+-------+ select * from t3; +----+-------+ | id | t4_id |

Hi, I've noticed that LLVM tends to generate suboptimal code and spill an excessive amount of registers in large functions, such as in those that are automatically generated by FFTW. LLVM generates good code for a function that computes an 8-point complex FFT, but from 16-point upwards, icc or gcc generates much better code. Here is an example of a sequence of instructions from a 32-point

The code below is intended to analyse a textbook example of a balanced incomplete block design: # # Data taken from pp. 219-230 in # Cox, D.R. (1958) Planning of Experiments. John Wiley and Son, Inc. New York. 308 pp. # day <- factor(rep(1:10, each = 3)) T <-

I have big .csv file. I would like to filter that file into a new table. For example, I have .csv file as below: f1 f2 f3 f4 f5 f6 f7 f9 f10 f11 t1 1 0 1 0 1 0 0 0 0 1 t2 1 0 0 0 0 1 1 1 1 1 t3 0 0 0 0 0 0 0 0 0 0 t4 1 0 0 0 1 0 0 0 0 0 t5 0 0 0 0 0 0 0 0 0 0 t6 0 0 0 0 0 0

# read.table is NOT part of the data.table package #library(data.table) DFM <- read.table( text= 'obs start end 1 2/1/2015 1/1/2017 2 4/11/2010 1/1/2011 3 1/4/2006 5/3/2007 4 10/1/2007 1/1/2008 5 6/1/2011 1/1/2012 6 10/5/2004 12/1/2004 ',header = TRUE, stringsAsFactors = FALSE) # cleaner way to compute D DFM$start <- as.Date( DFM$start, format="%m/%d/%Y" ) DFM$end

Hello, I've been redirected from cfe-dev, as code optimizations in clang are done in llvm layer. I'm investigating how optimized code clang generates, and have come across such an example: I have two procedures: void exec0(const int *X, const int *Y, int *res, const int N) { int t1[N],t2[N],t3[N],t4[N],t5[N],t6[N]; for(int i = 0; i < N; i++) { t1[i] = X[i]+Y[i];

Just to avoid possible confusion, let me correct a typo (at step [2] in the example below). Apologies! -----FW: <XFMail.111023084327.ted.harding at wlandres.net>----- Date: Sun, 23 Oct 2011 08:43:27 +0100 (BST) Sender: r-help-bounces at r-project.org From: (Ted Harding) <ted.harding at wlandres.net> To: r-help at r-project.org Subject: Re: [R] symmetric matrix multiplication On

(I'm sending to r-bugs because rpart is one of the recommended packages and is always installed. I'm also sending it directly to Dr. Ripley, as the maintainer.) rpart working as a classifier does not work (produces no splits) when the class indicator has no instances of one of the factor levels, as long as the factor level is not the final level. I have at least a partial fix, which I

Since the number of choices is small (6), how about this? Starting with Jeff's initial DFM: DFM <- structure(list(obs = 1:6, start = structure(c(16467, 14710, 13152, 13787, 15126, 12696), class = "Date"), end = structure(c(17167, 14975, 13636, 13879, 15340, 12753), class = "Date"), D = c(700, 265, 484, 92, 214, 57), bin = structure(c(6L, 3L, 5L, 1L, 3L, 1L), .Label

Ii is difficult to provide useful help, because you have failed to read and follow the posting guide. In particular: 1. Plain text, not HTML. 2. Use dput() or provide code to create your example. Text printouts such as that which you gave require some work to wrangle into into an example that we can test. Specifically: 3. Have you gone through any R tutorials?-- it sure doesn't look like

Hi all, I have a data set with time interval and depending on the interval I want to create 5 more variables . Sample data below obs, Start, End 1,2/1/2015, 1/1/2017 2,4/11/2010, 1/1/2011 3,1/4/2006, 5/3/2007 4,10/1/2007, 1/1/2008 5,6/1/2011, 1/1/2012 6,10/15/2004,12/1/2004 First, I want get interval between the start date and end dates (End-start). obs, Start , end, datediff