similar to: linear regression with dates

Hi, I have been examining large data and need to do simple linear regression with the data which is grouped based on the values of a particular attribute. For instance, consider three columns : ID, x, y, and I need to regress x on y for each distinct value of ID. Specifically, for the set of data corresponding to each of the 4 values of ID (76,111,121,168) in the below data, I should invoke

Dear Rui, >From my understanding of time-dependent covariates (not an expert but have been working on a similar problem), it would appear that the coding of the status column is not correct. Unless you have observed an event at each interval you should only have status=1 for the last interval. In your example I see 3 in total. Also, I think that if "end" is proportional to your

Hi, I've lost my mind on it... I have to scatterplot two vectors, grouped by a third variable, with two different dimensions according to whether each cell line in the plot is sensitive or resistant to a given drug, and with a different color for each of 9 tissues of origin. Here's what I've done:

Greetings, I would like to plot three lines on the same figure, and I am lost. There is an answer to a similar thread… but I tried matplot and it is beyond me. An example of the data follows: Year EM IM BM 1983 9.1 16.8 -7.7 1984 12.0 18.0 -6.0 1985 13.6 19.1 -5.5 1986 12.4 17.3 -4.9 1987 14.6 20.3 -5.7 1988 20.6 23.3 -2.6 1989 25.0 27.2 -2.2 1990 28.4 30.2 -1.8 1991 33.3 31.2 2.1 1992 40.6

Hi all, I have a question and I do not know If I am at the right place to ask this question. But is there someone that has the formula of CKD-Epi in code in R? I have tried a lot of loops, but none of the approaches give me the right answer. Is there someone who has this formula coded? Thank you! [[alternative HTML version deleted]]

Hi,* *I would like to set the entries of an array that appear as "NA" to a certain integer, but did not find a way to do so. The related code is given below: A <- array(-1000 , dim=c(100) ); for(i in 1: 100 ){ A[i] =B[i,57] - B[i,5]; } A [ which( A == "NA") ] <- 900; The "NA"s still appear as the entries of the array A. Can you tell me what is

Hi all, I am trying to detect association between a covariate and a disease outcome using R. This covariate shows time-varying effect, I add a time-covariate interaction item to build Cox model as follows: COX <- coxph(as.formula("Surv(TIME,outcome)~eGFR_BASE+eGFR_BASE:TIME"),ori.data); coef exp(coef) e(coef) z p eGFR_BASE

I want to make such plot using the following data set plot Mean with Sd(+/-) by Dim group by Lac here is data Does anyone know how to plot using xyploy Thanks, Aimin Mean Sd Var Min Max Dim Lac 704 44.00000 NA NA 44.0 44.0 -30 3 703 45.92000 9.5484030 91.1720000 60.0 34.4 -30 4 702 57.40000 NA NA 57.4 57.4 -30

Hello, I have the following data.frame structure(list(Study = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,

I have an attached data.frame with a binary composite variable (comp) and a predictor (egfr). I am trying to plot a restricted cubic spline using: require (Hmisc) require (Design) rcspline.plot(x=egfr,y=comp,model="logistic",nk=5,knots=NULL, xlab="label1", ylab="Label2", plim=c(0,1), plotcl=TRUE, showknots=TRUE,

1. When you don't know (or are not sure) what an object is, str() is your friend. 2. My guess is that `lidnames' is a character vector with names. 3. If genes.txt has only only column, you might as well use: mygenes <- scan("genes.txt", what="") which reads the data into mygenes as a character vector. Then your command should work. [read.table()

Hello, It's an alternative to use SAS by function in R? I want to plot d histograms by plot.from example bellow: Thank you! plot d 1 1 16.3 2 1 25.0 3 1 57.8 4 1 17.0 5 2 10.8 13 2 96.4 17 3 76.0 18 3 32.0 19 3 11.0 20 3 11.0 24 3 106.0 25 3 12.5 21 4 19.3 22 4 12.0 26 4 15.0 27 5 99.3 32 7 11.0 36

Dear all, Here is my code which am using to combine 5th column from different data sets. Here is the function to do my job genesymbol.append.file <-NULL gene.column <- NULL readGeneSymbol <- function(files,genesymbol.column=5){ for(i in fnames){ temp <- read.table(fnames,header=T,sep="\t",stringsAsFactors=F,quote="\"")

Dear Community, I have been using two types of survival programs to analyse a data set. The first one is an R function called aftreg. The second one an STATA function called streg. Both of them include the same analyisis with a weibull distribution. Yet, results are very different. Shouldn't the results be the same? Kind regards, J -- View this message in context:

I need to install VICIDIAL (cause is a free outbound call solution) and vicidial wants the following perl modules: MD5 Digest::MD5 Digest::SHA1 readline Bundle::CPAN DBI (found the rpm perl-DBI) DBD::mysql (found the RPM perl-DBD-Mysql) Net::Telnet Time::HiRes Net::Server OLE::Storage_Lite Spreadsheet::ParseExcel I do not want to taint the installation if an RPM is available on the DVD or in

Hi all, I have a data set with missing value. I would like to estimate those missing value by using normal ratio method. Below is part of my data: AS BL Serdang Jhr Phg Target station 0 0.0 12.8 0.0 23.7 0.0 6 0.0 81.7 0.2 0.0 NA 0 1.5 60.9 0.0 0.0 15.5 1 13.0 56.8 17.5 32.8 6.4 4 3.0 66.4

Dear All, I'am looking for examples showing that correlation does not imply causality, the targeted audience consists of undergraduate students (their first year at the university but in the BioMathStat track). All practicals are under R. I was able to extract this from R datasets: ### begin data(sunspots) data(lynx) spots <- window(sunspots, freq = 1, start = 1880, end = 1900) lnx <-

> On Dec 8, 2015, at 9:03 AM, Tru Huynh <tru at centos.org> wrote: > > Hi, > > > On Tue, Dec 08, 2015 at 08:24:34AM -0700, Wes James wrote: >> > <please trim your replies> > ... >> >> nslookup works: >> >> Non-authoritative answer: >> Name: mirrorlist.centos.org >> Address: 108.61.16.227 >> Name:

I have a DomU as a firewall (using the Dom0 kernel). I hid the three nics from Dom0 and assigned them to the firewall VM. The firewall VM sees the three nics no problem. Eth0 comes up and grabs an IP automatically from my DHCP server, however eth1 and eth2 are configured with static IPs but they do not come up with an IP: Here is my ''ifconfig'': [root@fw01 ~]# ifconfig eth0

Hi, I installed the shorewall 1.3.8-2 debian package to my debian testing machine which serves as the gateway to the internet. Since I have two other machine connect to internet thru this gateway machine, I also downloaded the configuration guide for "basic two-interface firewall" and followed the instructions. When I try to start the shorewall I get the following message and can not