Displaying 3 results from an estimated 3 matches for "fit42".
Did you mean:
fit2
2004 Dec 20
2
problems with limma
...model.matrix(~0+ts)
> colnames(design) <- levels(ts)
> fit4 <- lmFit(all,design)
> cont.matrix <- makeContrasts(
+ Baseline = akrC - nzwC,
+ NZW_Smk = nzwT - nzwC,
+ AKR_Smk = akrT - akrC,
+ Diff = (akrT - akrC) - (nzwT - nzwC),
+ levels=design)
> fit42 <- contrasts.fit(fit4,cont.matrix)
> fit42 <- eBayes(fit42)
> #
> topTable(fit42,coef="Baseline",adjust="fdr",number=5)
M t P.Value B
3189 942.0993 13.57485 0.004062283 -4.528799
8607 2634.1826 11.23476 0.006913442 -4....
2004 Dec 21
0
Fwd: problems with limma
...esign) <- levels(ts)
>> fit4 <- lmFit(all,design)
>> cont.matrix <- makeContrasts(
> + Baseline = akrC - nzwC,
> + NZW_Smk = nzwT - nzwC,
> + AKR_Smk = akrT - akrC,
> + Diff = (akrT - akrC) - (nzwT - nzwC),
> + levels=design)
>> fit42 <- contrasts.fit(fit4,cont.matrix)
>> fit42 <- eBayes(fit42)
>> #
>> topTable(fit42,coef="Baseline",adjust="fdr",number=5)
> M t P.Value B
> 3189 942.0993 13.57485 0.004062283 -4.528799
> 8607 2634.182...
2001 Jun 01
1
nls works but not gnls
This works fine:
fit42<-nls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal),
data=df,
start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6),
na.action=na.omit)
But this, identical except using gnls, doesn't converge:
fit43<-gnls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal),
data=df,
start=c(...