similar to: R-help Digest, Vol 30, Issue 6

I am trying to do a simple segregation analysis using the GAP package. I have the documentation for pointer but I desperately need an example so that I can see how to format the datfile and the jobfile. For each individual, I have FamilyId, SubjectId, FatherId, MotherId, and AffectedStatus (0/1). I would like to obtain the likelihood ratio statistic for transmission. I would greatly

Hi everyone I was using spline to interpolate single or two consecutive missing data points in time series. However, when it comes to longer gaps in the data the spline function generate new data for both my known and unknown data (see below). Aside from not understanding why this happens, I thought thought I might try function "aspline" in library (akima). However, I cannot install or

Hello everyone, I have a .csv file with the following format: uniqueID SubjectID Distance_miles Tag 1 1001 5.5 3 2 1001 7 1 3 1001 6.5 1 4 1001 5 1 5 1002

Hi, I am quite new to R and trying to analyze the following data. I have 28 controls and 25 patients. I measured X values of 4 different locations (A,B,C,D) in the brain image of each subject. And X ranges from 0 to 1. I think "control or patient" is a between subject factor and location is a within subject factor. So, controls: 28 patients: 25 (unbalanced data set) respone measure:

Last week (Wed 9/10/2003, "regression questions") I posted a question regarding the use of gnls( ) and its dissimilarity to the syntax that nls( ) will accept. No one replied, so I partly answered my own question by constructing indicator variables for use in gnls( ). The code I used to construct the indicators is at the end of this email. I do have a nagging, unanswered

Dear list, If have some repeated measurement data which looks something like: time <- rep(1:5 , each=2*4) groups <- rep(c("Case", "Control"), each=4) subjects <- factor(rep(1:(2*4), 5)) responses <- time + rnorm(5*2*4) + as.integer(factor(groups)) data <- data.frame(responses, time, groups, subjects) Now I want to plot each subject in a separate panel:

Muchas gracias, Carlos. Sobre todo por el sqldf, que seguro me ahorre mucho tiempo. La opción que se plante de primero d <- table(tips$day) y luego dim(d) me parece menos eficiente y cómo que directamente sqldf("select count(distinct day) from tips"), pero supongo que esos son gustos! También son "cómodas" las líneas: aggregate(subjectid ~ cond, data = dat, FUN = function(x)

Estoy en 3.3.0 y "sqldf" lo instala sin problemas... El 23 de junio de 2016, 1:23, Mauricio Monsalvo <m.monsalvo en gmail.com> escribió: > Malas nuevas para mi: > package ?sqldb? is not available (for R version 3.3.0) > ¿Puedo hacer algo más que esperar? No me voy a "bajar" de versión de R. > > El 22 de junio de 2016, 20:02, Mauricio Monsalvo

Hi I have a dataframe in which some subjects appear in more than one row. I want to extract the subject-rows which have the minimum date per subject. I tried the following aggregate function. attach(dataframe.xy) aggregate(Date,list(SubjectID),min) Unfortunately, the format of the Date-column changes to numeric, when I'm applying this function. How can I preserve the date format? Thanks

With pleasure we announce the release of tinc version 1.0.24. Here is a summary of the changes: * Various compiler hardening flags are enabled by default. * Updated support for Solaris, allowing switch mode on Solaris 11. * Configuration will now also be read from a conf.d directory. * Various updates to the documentation. * Tinc now forces glibc to reload /etc/resolv.conf after it

With pleasure we announce the release of tinc version 1.0.24. Here is a summary of the changes: * Various compiler hardening flags are enabled by default. * Updated support for Solaris, allowing switch mode on Solaris 11. * Configuration will now also be read from a conf.d directory. * Various updates to the documentation. * Tinc now forces glibc to reload /etc/resolv.conf after it

Hello, I have been searching for ways to obtain these for combinations of fixed factors and levels other than the 'baseline' group (contrasts coded all 0's) from a mixed-effects model in lme. I've modelled the continuous variable y as a function of a continuous covariate x, and fixed factors A, B, and C. The fixed factors have two levels each and I'd like to know whether

Hi everybody, I'm fairly new to survival analysis with R and have some questions how to apply and interpret the coxph and related functions: I have time-dependent covariates with several measurements per subject with constant delta t. The covariates change in each time step. I fitted the following model: fit <- coxph(Surv(start, stop, event) ~ ratePO + rateC + BLamp + BLP80 +

Hi, I want to mirror NetBSD on this server: ftp://ftp.univ-orleans.fr/logiciel-libre/netbsd-ftp/ The problem is I dont want to mirror old versions (1.5.x) There are 1.5* branches as subdirs of 'packages' and others at the root of the server called 'NetBSD-release-1-5' and 'NetBSD-1.5.3'. I want to avoid to sync them. The command I use to sync is: - - - - - - - - - - -

Hello Steve, Hello OpenSSH-portable developers, I am building OpenSSH for our (EBCDIC-based) BS2000 mainframe operating system, and I noticed you do the same for OS/390. Because my initial ssh port was based on IBM's OSS port (ssh-1.2.2 or some such), I thought it was fair enough to help with a little co-operation; we might come up with a unified EBCDIC patch which could be contributed to

I am trying to plot where data points from a give subject are connected by a line. Each subject is represented by a single row of data. Each subject can have up to five observations. The first five columns of mydata give the time of observation, columns 6-10 give the values at each time point. Some subjects have all data, some are missing values. The code I wrote to draw the plot is listed below.

So far I have figured out that the following line reads our time series files into R OK. dtLs$dta <- read.table("C:/TryRRead/datFiles/JFeqfi4h.rta", header = TRUE, sep = ",", colClasses = "character") But I have to remove a main-title line so that the first line is the column titles line. This leads to having two sets of data files around when we would rather

Hola Necesito algo simple pero no logro resolverlo en R: contar elementos distintos (unique) pero agrupados. En SQL sería: SELECT Camara, count(DISTINCT(Laboratorio)) as Lab FROM [datos] GROUP BY Camara. Puedo hacerlo para toda la matriz de datos con length(unique(datos$Laboratorio)) pero no logro hacerlo agrupado por una variable (factor) dada. Estuve intentando con el paquete data.table pero me

On Sun, Oct 29, 2017 at 11:51:55PM +0100, Fengguang Wu wrote: >Hi Linus, > >Up to now we see the below boot error/warnings when testing v4.14-rc6. > >They hit the RC release mainly due to various imperfections in 0day's >auto bisection. So I manually list them here and CC the likely easy to >debug ones to the corresponding maintainers in the followup emails. >

On Sun, Oct 29, 2017 at 11:51:55PM +0100, Fengguang Wu wrote: >Hi Linus, > >Up to now we see the below boot error/warnings when testing v4.14-rc6. > >They hit the RC release mainly due to various imperfections in 0day's >auto bisection. So I manually list them here and CC the likely easy to >debug ones to the corresponding maintainers in the followup emails. >