search for: lwt

...I am reading in text file data prepared in Access database by someone. One of the field contains timestamp data, how can I separate the timestamp data into two varaibles: date and time. Can I specify the field is in timestamp format when I first reading in ? My reading in data are as below: 449 LWT 22/10/2003 15:43:00 441 143 449 LWT 17/11/2003 15:25:00 421 169 449 LWT 14/11/2003 15:04:00 454 166 449 LWT 17/12/2003 5:55:00 428 199 449 LWT 7/12/2003 15:28:00 452 189 449 LWT 15/11/2003 15:20:00 457 167 Thanks in advance. Stella

...levels - "white", "black" and "other") in my logistic regression model, the correlation matrix of the coefficients gets messed up. I get something like: ----------------------------------------- Correlation of Coefficients: ( A L RACEb AGE , 1 LWT , 1 RACEblack 1 RACEother . . attr(,"legend") [1] 0 ` ' 0.3 `.' 0.6 `,' 0.8 `+' 0.9 `*' 0.95 `B' 1 ------------------------------------- I couldn't figure out how to interpret it. Here is the sequence of commands and the complet...

...#define function involving regression model with several covariates custom <- function(outcome, exposure, covar1, covar2, covar3){ model <- lm(outcome ~ exposure + covar1 + covar2 + covar3) expected <- predict(model) summary(expected) } library(MASS) attach(birthwt) custom(bwt, lwt, low, age, race) #Works when 3 covariates are specified custom(bwt,lwt,low,age) # Does not work with < or > 3 covariates varlist <- list(low,age,race) custom(bwt,lwt, varlist) #Does not work if covariates are included as a list Thanks very much for your help Tim -- Tim Elwell-Sutton...

...t Study # Described in Section 7.3 Page 230 # # Name Description Codes/Values # pair Pair 1 - 56 # low Low Birth Weight 1 = BWT<=2500g, # 0 = BWT>2500g # age Age of Mother Years # lwt Weight of Mother at Pounds # Last Menstrual Period # race Race 1 = White # 2 = Black # 3 = Other # smoke Smoking Status 0 = No,1 = Yes # During Pregnancy #...

...logistic regression in R. But it seems that I am having convergence problem. Am I doing anything wrong, or just the data is off? Here is my code in R - lbw <- read.table("http://www.biostat.jhsph.edu/~ririzarr/Teaching/754/lbw.dat" , header=TRUE) attach(lbw) lbw[1:2,] low age lwt race smoke ptl ht ui ftv bwt 1 0 19 182 2 0 0 0 1 0 2523 2 0 33 155 3 0 0 0 0 3 2551 #-----R implementation of logistic regression : gradient descent ------ sigmoid<-function(z) { 1/(1 + exp(-1*z)) } X<-cbind(age,lwt, smoke, ht, ui) #y<-low my_logistic<-function(X,y) { alph...

...anks, Patrick MY R CODE IS: Dataset <- read.table("<path>", header=TRUE, sep="", na.strings="NA", dec=".", strip.white=TRUE) Dataset$race <- factor(Dataset$race, levels=c('other','black','white')) GLM.1 <- glm(low ~ lwt + ptl + ht + race + smoke , family=binomial(logit), data=Dataset) summary(GLM.1) MY SAS CODE IS: PROC LOGISTIC descending DATA=p2; class race (ref='other'); MODEL LOW = lwt ptl ht race smoke / lackfit parmlabel expb link=logit; RUN; MY R OUTPUT IS: Coefficients: Estima...

...he values for PostProbMod(BF), one can say that higher posterior probability indicates the better model. Once again, thank you very much for the time you have given. I am looking forward for your reply. Regards, Deb R Code: library(MCMCpack) data(birthwt) model1 <- MCMCregress(bwt~age+lwt+as.factor(race) + smoke + ht, data=birthwt, b0=c(2700, 0, 0, -500, -500, -500, -500), B0=c(1e-6, .01, .01, 1.6e-5, 1.6e-5, 1.6e-5, 1.6e-5), c0=10, d0=4500000, marginal.likelihood="Chib95", mcmc=50000) model2 <- MCMCregress(bwt~age+lwt+as.factor(race) + smoke, data=birthwt, b0=c(2700, 0...

...EV SD condEV cond SD model 1 model 2 model 3 model 4 model 5 Intercept 100 0.1841 1.2204 0.184 1.220 1.017 1.175 -0.853 -1.057 0.532 age 17.8 -0.0113 0.0285 -0.063 0.036 . . . . -0.071 lwt 50.0 -0.0079 0.0093 -0.016 0.007 -0.017 - 0.017 . . . smokeTRUE 9.5 0.0469 0.1798 0.496 0.345 . . . . . ptdTRUE 99.4 1.5161 0.4751 1.526 0.461 1.407 1.596 1.732 1.46...

Fellow R Users: I'm not extremely familiar with lda or R programming, but a recent editorial review of a manuscript submission has prompted a crash cousre. I am on this forum hoping I could solicit some much needed advice for deriving a classification equation. I have used three basic measurements in lda to predict two groups: male and female. I have a working model, low Wilk's lambda,

...to function The example listed on the help for LME worked without problem. I don't know if my last error reports got out of Mexico or not so here they are again. 3. The following code on page 236 of Venables and Ripley 1997 birthwt.step2 <- step(birthwt.glm,~.^2+I(scale(age)^2)+I(scale(lwt)^2),trace=F) produced the following warning message 8 times Warning: ignored non function "scale" 4. The following code on page 358 of Venables and Venables 1997 plot(survfit(leuk.coxs),lty=2:3,lwd=3,add=T,log=T) produced the following warning message Warning: parameter "add&...

The birth weight example from ?stepAIC in package MASS runs well as indeed it should. However when I change stepAIC() calls to step() calls I get warning messages that I don't understand, although the output is similar. Warning messages: 1: In model.response(m, "numeric") : using type="numeric" with a factor response will be ignored (and three more the same.) Checked

...rom my humble point of view) is that all loading is done statically: Loading the R interpreter is done statically. Loading the xts library is done statically... etc... See below. Hopefully, one may one day consider R to be a statically-typed, type-inferred, compiled, statistical language, with Lwt-style multithreading. :) > yziquel at seldon:~/git/ocamlr-xts$ ocamlbuild xts.cmo > Finished, 1 target (0 cached) in 00:00:00. > + ocamlfind ocamlc -c -package R.interpreter -o xts.cmo xts.ml > File "xts.ml", line 37, characters 4-9: > Warning P: this pattern-matching is...

...Ability to inspect (read-only) the inner structure of R values, which is quite convenient: you get to know rather quickly what a given piece of R code returns, which you need to know to type R code statically in order to bind it to OCaml. - Not thread-safe at all. At least, not more than R is... Lwt-style multithreading of R code could be possible, modulo some simple and deep (i.e. below R API) changes in the R evaluation functions. POSIX threading a single R thread with multiple OCaml threads is not yet possible, but is within reach. - Doesn't interact well the R "Defaults"...

Hi everyone, I have some question about the Xenopsd. I’ve read the wiki document but cannot understand where to use it, Is this different with XAPI ? What’s the main work for? Can anyone describe it more clearly about it? Thanks a lot!! _______________________________________________ Xen-users mailing list Xen-users@lists.xen.org http://lists.xen.org/xen-users

Hey guys, I was playing around with the command "rake db_schema_dump" and kept getting this error "no such file or directory - /tmp/mysql.sock". I tried looking around to find this mysql.sock file; but no luck so far. Question: What is mysql.sock file? (how does it get created and how does it get used) I installed a MySQL DB and could not find this file. Could someone help

The following series of patches fill in most of the gaps in the OCaml bindings to libxl, to make them useful for clients such as xapi/xenopsd (from XCP). There are a number of bugfixes to the existing bindings as well. I have an experimental version of xenopsd that successfully uses the new bindings. An earlier version of the first half of the series was submitted to the last by Ian Campbell on

Friends, I have updated my Asterisk 1.2 presentation with the latest information. It is still available in the same place as before: http://www.astricon.net/asterisk1-2/ Please continue to test the beta of Asterisk 1.2, available at ftp.digium.com. We need all the feedback we can get. If you are a developer and have some time for community work, please check in with the bug tracker and help us

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking