Displaying 20 results from an estimated 929 matches for "b3".
Did you mean:
3b
2003 Nov 18
1
aov with Error and lme
...t;a2","a2","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3","a3")
factB<-c("B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2","B3","B1","B2...
2011 Jul 05
3
plotting survival curves (multiple curves on single graph)
...ch.
--Trey
The function and calls below use the data in this Excel file (feel
free to access):
https://docs.google.com/leaf?id=0B5zZGW2utJN0ZDk1NjA0ZjUtMWU0ZS00ZGQ3LWIxZTUtOWE0NGVmYWMxODJl&hl=en_US
## - plot Siler survival curve
##############################
silsurv<-function(a1,b1,a2,a3,b3)
{
sil=function(t)
{
h.t<-a1*exp(-b1*t)+a2+a3*exp(b3*t)
S.t<-exp(-a1/b1*(1-exp(-b1*t))-a2*t+a3/b3*(1-exp(b3*t)))
d.t<-S.t*h.t
#return(d.t)
return(S.t)
#return(h.t)
}
t<-seq(0,90,1)
plot(t,sil(t),ylim=c(0,1),type='...
2016 Feb 13
2
Code in headers
> On Feb 11, 2016, at 12:43 AM, via llvm-dev <Alexander G. Riccio> wrote:
>
> I don’t think that we can agree to abstract code guidelines without knowing what it means in practice for the codebase. If you’re interested in this, please include a diff that shows the impact to the headers, and we should also measure what happens to the performance of the generated compiler.
>
>
2012 Jan 12
2
Function accepted by optim but not mle2 (?)
Dear Sir/ Madam,
I'm having trouble de-bugging the following - which works perfectly
well with optim or optimx - but not with mle2.
I'd be really grateful if someone could show me what is wrong. Many
thanks in advance. JSC:
gompertz<- function (x,t=data)
{
a3<-x[1]
b3<-x[2]
shift<-data[1]
h.t<-a3*exp(b3*(t-shift))
S.t<-exp(a3/b3*(1-exp(b3*(t-shift))))
PF<-1-(exp(a3/b3*(1-exp(b3*(data[1]-shift)))))
QF<-1-(exp(a3/b3*(1-exp(b3*(data["length"(data)]-shift)))))
return(S.t*h.t/(QF-PF))
}
## GOMPERTZ - maximum lik...
2011 Jun 28
1
plotting survival curves with model parameters
...z-Makeham, with the addition of a
juvenile component that includes two parameters---one that describes
the initial infant mortality rate, and a negative exponential that
describes typical mortality decline over the juvenile period. The
entire hazard is expressed as
h(x) = a1*exp(-b1*x)+a2+a3*exp(b3*x)
I've had success in plotting the curves using the following function:
############################################################
siler<-function(x=c(0.1,0.5,0.001,0.003,0.05)) # model Siler parameters
{
sil=function(t)
{
a1 = x[1]
b1 = x[2]
a2 =...
2005 Oct 07
3
Converting PROC NLMIXED code to NLME
...convert the following NLMIXED code to NLME, but am
running into problems concerning 'Singularity in backsolve'. As I am new
to R/S-Plus, I thought I may be missing something in the NLME code.
NLMIXED
***********
proc nlmixed data=kidney.kidney;
parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43
varu=0.5;
eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;
hazard=eta+log(delta)+log(gamma)+(gamma-1)*log(rtime);
survivor=(-exp(eta))*delta*(rtime**gamma);
ll=(event*hazard)+survivor;
model rtime ~ general(ll);
random u~normal(0,varu) subject=patient out=frailty;
run;
NLME
********
k...
2016 Feb 16
2
WebKit B3 (was LLVM Weekly - #110, Feb 8th 2016)
> On Feb 15, 2016, at 4:25 PM, Philip Reames <listmail at philipreames.com> wrote:
>
> After reading https://webkit.org/blog/5852/introducing-the-b3-jit-compiler/ <https://webkit.org/blog/5852/introducing-the-b3-jit-compiler/>., I jotted down a couple of thoughts of my own here: http://www.philipreames.com/Blog/2016/02/15/quick-thoughts-on-webkits-b3/ <http://www.philipreames.com/Blog/2016/02/15/quick-thoughts-on-webkits-b3/>
Thank...
2012 May 13
2
Discrete choice model maximum likelihood estimation
...e to add
them as well.
I am not familiar with the syntax I have to write in the likelihood
functions, so I really doubt that they are true. Below I simplify the
problem and provide the code I've written:
# Probabilities for discrete choices for a=3, a=2 and a=1 respectively
P3 <- function(b3,b,t) {
P <- exp(b3+b*(t==1))/(1-exp(b3+b*(t==1)))
return(P)
}
P2 <- function(b2,b,t) {
P <- exp(b2+b*(t==1))/(1-exp(b2+b*(t==1)))
return(P)
}
P1 <- function(b1,b,t) {
P <- exp(b1+b*(t==1))/(1-exp(b1+b*(t==1)))
return(P)
}
# Likelihood functions for discrete choices for a=3,...
2003 Oct 04
2
mixed effects with nlme
...ME and the
last version of R. More concretely, I have a repeated measures design with
a single group and 2 experimental factors (say A and B) and my interest is
to compare additive and nonadditive models.
suj rv A B
1 s1 4 a1 b1
2 s1 5 a1 b2
3 s1 7 a1 b3
4 s1 1 a2 b1
5 s1 4 a2 b2
6 s1 2 a2 b3
7 s2 6 a1 b1
8 s2 8 a1 b2
9 s2 10 a1 b3
10 s2 3 a2 b1
11 s2 6 a2 b2
12 s2 6 a2 b3
13 s3 1 a1 b1
14 s3 6 a1 b2
15 s3 5 a1 b3
16 s3 3...
2012 Apr 19
3
How to "flatten" a multidimensional array into a dataframe?
Hi,
I have a three dimensional array, e.g.,
my.array = array(0, dim=c(2,3,4), dimnames=list( d1=c("A1","A2"),
d2=c("B1","B2","B3"), d3=c("C1","C2","C3","C4")) )
what I would like to get is then a dataframe:
d1 d2 d3 value
A1 B1 C1 0
A2 B1 C1 0
.
.
.
A2 B3 C4 0
I'm sure there is one function to do this transformation, I just don't
know which one.
Thanks for your help,...
2012 Mar 12
3
how to calculate a variance and covariance matrix for a vector
Hello,
I have a vector {a, b1, b2, b3, b4}. How can I calculate the following
matrix:
var(a) cov(a, b1) cov(a, b2) cov(a, b3) cov(a, b4)
cov(a, b1) var(b1) cov(a, b2) cov(a, b3) cov(a, b4)
...
...
cov(a, b1) cov(a, b2) cov(a, b3) cov(a, b4) var(b4)
I would very appreciate your inputs. Thank you very much.
Sincerely,
Jialin...
2012 Oct 18
7
summation coding
I would like to code the following in R: a1(b1+b2+b3) + a2(b1+b3+b4) +
a3(b1+b2+b4) + a4(b1+b2+b3)
or in summation notation: sum_{i=1, j\neq i}^{4} a_i * b_i
I realise this is the same as: sum_{i=1, j=1}^{4} a_i * b_i - sum_{i=j} a_i
* b_i
would appreciate some help.
Thank you.
--
View this message in context: http://r.789695.n4.nabble.com/sum...
2016 Feb 15
5
WebKit B3 (was LLVM Weekly - #110, Feb 8th 2016)
...; On Feb 9, 2016, at 9:55 AM, Rafael Espíndola via llvm-dev <llvm-dev at lists.llvm.org> wrote:
> >
> > JavaScriptCore's [FTL JIT](https://trac.webkit.org/wiki/FTLJIT <https://trac.webkit.org/wiki/FTLJIT>) is moving away
> > from using LLVM as its backend, towards [B3 (Bare Bones
> > Backend)](https://webkit.org/docs/b3/ <https://webkit.org/docs/b3/>). This includes its own [SSA
> > IR](https://webkit.org/docs/b3/intermediate- <https://webkit.org/docs/b3/intermediate-representation.html>representation.html <https://webkit.org/docs/b3/i...
2008 Jun 24
2
subscripting in data frames with NA
Dear list:
Given
> str(b3)
'data.frame': 159 obs. of 6 variables:
$ index_pollution : num 8.228 10.513 0.549 0.915 10.416 ...
$ position_descrip: chr "2" "2" "2" NA ...
$ position_geo : chr "3" "0" "3" "3" ...
$ institution :...
2007 Sep 21
1
Error using nls()
...547, 656, 628, 627), d1 = c(1, 1, 1, 1, 1, 1, 1.05,
1.05, 1.05, 1.15, 1.15, 1.15)), .Names = c("q1", "d1"), row.names =
as.integer(c(NA,
12)), class = "data.frame")
I tried to estimate a nonlinear regression model using nls(). I get
>
nls(log(q1)~b1+(b2/b3)*(d1^b3-1),data=x,start=list(b1=0,b2=1,b3=1),trace=TRUE)
246.132 : 0.0000 1.0000 1.0000 51.0907
Error in qr.solve(QR.B, cc) : singular matrix 'a' in solve
However, using:
> jjf <- function(x){z <- log(q1)-x[1]+(x[2]/x[3])*(d1^x[3]-1);sum(z*z)}
> optim(c(0,1,1),jjf)
rende...
2020 Jun 30
2
Samba AD + adblocking in bind9
...And I include the named.conf.ads in my named.conf.local:
*include "/etc/bind/named.conf.ads";include
"/var/lib/samba/bind-dns/named.conf";*
This works. When I try to click on an ad, I get redirected to localhost,
which is fine. Samba, complains however:
Jun 30 06:37:34 bubba-b3-two systemd[1]: Started Samba AD Daemon.
Jun 30 06:37:34 bubba-b3-two winbindd[3237]: [2020/06/30 06:37:34.807028,
0] ../lib/util/become_daemon.c:138(daemon_ready)
Jun 30 06:37:34 bubba-b3-two winbindd[3237]: daemon_ready: STATUS=daemon
'winbindd' finished starting up and ready to serve...
2003 Sep 07
0
µÚËĽ챱¾©³¯Ñô¹ú¼ÊÉÌÎñ½Ú
=BE=B4=C6=F4=D5=DF=A3=BA
=A1=B0=B5=DA=CB=C4=BD=EC=B1=B1=BE=A9=B3=AF=D1=F4=B9=FA=BC=CA=C9=CC=CE=F1=
=BD=DA=A1=B1=BD=AB=D3=DA2003=C4=EA10=D4=C2=D4=DA=B1=B1=BE=A9=BE=D9=B0=EC=A1=
=A3=BD=EC=CA=B1=A3=AC=B9=FA=BC=D2=C1=EC=B5=BC=C8=CB=A1=A2=CD=E2=B9=FA=D5=FE=
=D2=AA=BA=CD=C0=B4=D7=D4=C3=C0=A1=A2=B5=C2=A1=A2=B7=A8=A1=A2=B0=C4=B4=F3=C0=
=FB=D1=C7=B5=C8=B9=FA=B5=C4=C9=CC=...
2011 May 30
2
Basic question about three factor Anova
Read the data using scan():
#
# a1 a2 a3 a4
# ------------- ------------- ------------- -------------
# b1 b2 b3 b1 b2 b3 b1 b2 b3 b1 b2 b3
# --- --- --- --- --- --- --- --- --- --- --- ---
#
# c1:
# 4.1 4.6 3.7 4.9 5.2 4.7 5.0 6.1 5.5 3.9 4.4 3.7
# 4.3 4.9 3.9 4.6 5.6 4.7 5.4 6.2 5.9 3.3 4.3 3.9
# 4.5 4.2 4.1 5.3...
2016 Feb 16
2
WebKit B3 (was LLVM Weekly - #110, Feb 8th 2016)
...ew Trick wrote:
>>
>>> On Feb 15, 2016, at 4:25 PM, Philip Reames < <mailto:listmail at philipreames.com>listmail at philipreames.com <mailto:listmail at philipreames.com>> wrote:
>>>
>>> After reading https://webkit.org/blog/5852/introducing-the-b3-jit-compiler/ <https://webkit.org/blog/5852/introducing-the-b3-jit-compiler/>., I jotted down a couple of thoughts of my own here:http://www.philipreames.com/Blog/2016/02/15/quick-thoughts-on-webkits-b3/ <http://www.philipreames.com/Blog/2016/02/15/quick-thoughts-on-webkits-b3/>
>>...
2010 Feb 09
2
How can I rearange my dataframe
An embedded and charset-unspecified text was scrubbed...
Name: ???????????
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100209/db12d14b/attachment.pl>