similar to: GAP pointer

On Fri, 5 Aug 2005 Julia Reid wrote: > Subject: [R] GAP pointer > > 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

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

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 +

I'm getting 2 errors which are probably related, the second of which I haven't actually seen an issue with yet (I just found it looking at the error log) The issue is that I cannot change my subscriptions through Thunderbird (version 31.0). I created the subscriptions through Thunderbird, but now it won't let me change them. dovecot responds with sutff like this when I sniff

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 On Sun, 3 Aug 2014, Joe Reid wrote: > The time stamp on the .subscriptions.lock file gets updated every time I try > to make a change. To test I edited .subscriptions by hand and forced > Thunderbird to try to subscribe something and I got access to the folder I > added, so I presume it is reading the file correctly. > >

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: - - - - - - - - - - -

I have a time series of about 1500 measurements. There are sporadic data gaps. I would like to plot the time series with type "line". I only want a line to connect the periods with data; I don't want a line to connect the points across a data gap. Is there a function or recommended method to deal with this? For example there are 2000 days days 1:1000 each have one

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

rPi3 running current Buster.  I get this far: # upsc networkups Init SSL without certificate database battery.charge: 100 battery.charge.warning: 30 battery.runtime: 2891 battery.type: PbAC device.mfr: Vertiv Co device.model: Liebert PST5 device.type: ups driver.name: usbhid-ups driver.parameter.pollfreq: 30 driver.parameter.pollinterval: 2 driver.parameter.port: auto driver.parameter.productid:

Dear All: lme(sbp ~ cengirth, data = family, random= ~ 1 | familyid) converges but lme(sbp ~ cengirth, data = family, random= ~ cengirth | familyid) does not. I get the following message: Error in lme.formula(sbp ~ cengirth, data = family, random = ~cengirth | : iteration limit reached without convergence (9) The data has 488 rows and 154 familyid levels. For

Hi, I need to subscribe and use several Polycom-Kirk 5020 handsets along non-Polycom-Kirk handsets on a one-cell Polycom-Kirk 600/3 base station. Has anyone tried this ? Which values did you pick for Subscription mode (with or without Account Code) and IPEI ? Regards -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi, I am using the gap R package to do the Hardy Weinberg Case Control test for many SNP. I am not sure what the values initial1 and initial2 should be for the test. I tried values but they failed. I emailed the author but to no avail. There seems to be some documentation that is deleted at the top, if anyone can direct me how to get this I will be grateful. -- Thanks, Jim. [[alternative HTML