Displaying 3 results from an estimated 3 matches for "ctst".
Did you mean:
cost
2007 Apr 25
1
for loops
Hello everybody
I'm very new at using R so probably this is a very stupid question.
I have a matrix of "p" columns and I have to calculate for each of them the "two sample t-statistic" and p-value and to save the results into two different vectors.
I have divided my matrix into two submatrices: submatrix A containing the first "n1" rows (p columns) and submatrix B
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data?
...;-which(air$stratn==lest[i])
for (j in 1:length(a))
{
air$nctrl[a[j]]<-sum(air$mi[a[-j]])
}
}
#create cases and controls
cami<-rep(1,sum(mi))
ctmi<-rep(0,sum(air$nctrl))
capm<-rep(pm10, mi)
ctpm<-rep(pm10, air$nctrl)
cast<-rep(air$stratn, mi)
ctst<-rep(air$stratn, air$nctrl)
cady<-rep(dates, mi)
ctdy<-rep(dates, air$nctrl)
cases<-c(cami, ctmi)
stranu<-c(cast, ctst)
days<-c(cady, ctdy)
pmva<-c(capm, ctpm)
air2<-data.frame(cases, days, pmva,stranu)
air.cl<-clogit(cases~pmva + strata(stranu), method=&...
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data? Unformatted text-version, with an additional note
...gth(lest))
{
?? a<-which(air$stratn==lest[i])
?? for (j in 1:length(a))
?? {
????? air$nctrl[a[j]]<-sum(air$mi[a[-j]])
?? }
}
#create cases and controls
cami<-rep(1,sum(mi))
ctmi<-rep(0,sum(air$nctrl))
capm<-rep(pm10, mi)
ctpm<-rep(pm10, air$nctrl)
cast<-rep(air$stratn, mi)
ctst<-rep(air$stratn, air$nctrl)
cady<-rep(dates, mi)
ctdy<-rep(dates, air$nctrl)
cases<-c(cami, ctmi)
stranu<-c(cast, ctst)
days<-c(cady, ctdy)
pmva<-c(capm, ctpm)
air2<-data.frame(cases, days, pmva,stranu)
air.cl<-clogit(cases~pmva + strata(stranu), method="approxima...