R help - Apr 2016 - R 3.2.5 is released

The 3.2.4-revised version turned out to give trouble for some of CRAN's
subsystems.

Accordingly, a rebadged version 3.2.5 is now released; it only differs in the
version number and a few clean-up items. If you have a working install of
3.2.4-revised there should be no reason to upgrade it.

You can get the source code from

http://cran.r-project.org/src/base/R-3/R-3.2.5.tar.gz

or wait for it to be mirrored at a CRAN site nearer to you.

For the R Core Team,

Peter Dalgaard

New md5 sums are

MD5 (AUTHORS) = eb97a5cd38acb1cfc6408988bffef765
MD5 (COPYING) = eb723b61539feef013de476e68b5c50a
MD5 (COPYING.LIB) = a6f89e2100d9b6cdffcea4f398e37343
MD5 (FAQ) = cd7bfa946b8650cb87807e94e46984c6
MD5 (INSTALL) = 3964b9119adeaab9ceb633773fc94aac
MD5 (NEWS) = 42d76ce7f8e80977d5043bca4234f4c9
MD5 (NEWS.0) = bfcd7c147251b5474d96848c6f57e5a8
MD5 (NEWS.1) = eb78c4d053ec9c32b815cf0c2ebea801
MD5 (NEWS.2) = e840d32b7ef7a7603455d30d6d54fda7
MD5 (R-latest.tar.gz) = 7b23ee70cfb383be3bd4360e3c71d8c3
MD5 (README) = aece1dfbd18c1760128c3787f5456af6
MD5 (RESOURCES) = 529223fd3ffef95731d0a87353108435
MD5 (THANKS) = ba00f6cc68a823e1741cfa6011f40ccb
MD5 (VERSION-INFO.dcf) = 661c20647b9f0b9fa980297c53438609
MD5 (R-3/R-3.2.5.tar.gz) = 7b23ee70cfb383be3bd4360e3c71d8c3


The relevant NEWS file entries (relative to 3.2.4) are

CHANGES IN R 3.2.5:

  BUG FIXES:

    ? format.POSIXlt() behaved incorrectly in R 3.2.4.  E.g. the output
      of format(as.POSIXlt(paste0(1940:2000,"-01-01"), tz =
"CET"),
      usetz = TRUE) ended in two "CEST" time formats.

    ? A typo in the Makefile for src/extra/xz prevented builds of
      liblzma.a. (Notice that this will become unbundled in 3.3.0.)



-- 
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Office: A 4.23
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com

Howdy everyone



I?m trying to get Odds ratio and OR confidence intervals using a probit model,
but I'm not getting.

 

Do you think you can help me?

 

I?m new with R L

 

naive                   =
summary(glm(pcr.data[,7]~boldBeta_individual+pcr.data$age,family=binomial(link=probit)))

naive_answer            = c(naive$coefficients[,1:3])                          
#naive estimates for

                                                                               
#alpha (first 4 collumns: intercept; beta_intercept, beta_slope and age) and

                                                                               
#and SE(last 4 collumns: intercept; beta_intercept, beta_slope and age)

 

OR.naive = exp(1.6*coef(naive))

 

(till here works, the problem is with the confidence interval)

 

I tried to get the Standard error from the variance, but I?m not sure if this
can be done as I?ve done.

 

 

Var_coef <- 1.6^2*var(coef(naive))

SE_coef <- Var_coef/sqrt(nsample)                    ########## I thi k this
is correct

 

OR.naive.inf <- exp(OR.naive - (1.96 * SE_coef))

OR.naive.sup <- exp(OR.naive + (1.96 * SE_coef))

 

if I used logit link I would get the CI with confint(na?ve) command, but with
probit I don't think so. Is there a way?

 

What should I do?

 

 

 


Atenciosamente,
Rosa Oliveira

-- 
____________________________________________________________________________



Rosa Celeste dos Santos Oliveira, 

E-mail: rosita21 at gmail.com <mailto:rosita21 at gmail.com>
Tlm: +351 939355143 
Linkedin: https://pt.linkedin.com/in/rosacsoliveira
<https://pt.linkedin.com/in/rosacsoliveira>
____________________________________________________________________________
"Many admire, few know"
Hippocrates