Displaying 20 results from an estimated 931 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,...
2025 Apr 01
1
HCL Tripp-Lite SMARTINT2200VS supported by tripplite_usb
...00 00
?56 32 30 38 33 58 58 0D
?3A 55 AA 0D 00 00 00 00
?55 FF FF 0D 00 00 00 00
?3A 53 AC 0D 00 00 00 00
?53 31 34 30 00 64 31 0D
?3A 44 BB 0D 00 00 00 00
?44 37 46 38 37 0D 00 00
?3A 4D B2 0D 00 00 00 00
?4D 37 39 38 33 0D 00 00
?3A 54 AB 0D 00 00 00 00
?54 37 35 31 46 33 30 0D
?3A 4C B3 0D 00 00 00 00
?4C 30 30 30 44 00 58 0D
?3A 44 BB 0D 00 00 00 00
?44 37 46 38 38 0D 00 00
?3A 53 AC 0D 00 00 00 00
?53 31 34 30 00 64 31 0D
?3A 44 BB 0D 00 00 00 00
?44 37 46 38 38 0D 00 00
?3A 4D B2 0D 00 00 00 00
?4D 37 39 38 33 0D 00 00
?3A 54 AB 0D 00 00 00 00
?54 37 35 31 46 33 30 0...
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/>
>>...