similar to: Should NA's in summary() output always be reported???

What is the correct form to write statement meaning: if (a==1) {b=2; c=3}; else {b=0; c=0}; Thank you Jue Wang, Biostatistician Contracted Position for Preclinical & Research Biostatistics PrO Unlimited (908) 231-3022

Dear all, I have 2 matrices, one is a 8x3 matrix, called X; the other matrix is a 3x3 indicator matrix with the diagonal element as 0 or 1. when a variable is included in the model, the corresponding diagonal element is 1, otherwise, it is 0. Let A be a set of matrices that contain the possible indicator matrix. e.g., A= [A1, A2, A3], where A1= [1,0,0;0,0,0,0,0,0], A2 =[1,0,0;0,1,0,0,0,0], A3

All, I am trying to include a day of week variable (1-7) in in a regression model. I would like to have the day of week treated as a categorical variable rather than a number the code looks like lm( dep ~ WKDY) I know this is a basic question, but help would be appreciated thanks spencer [[alternative HTML version deleted]]

Hi, I have a matrix with few hundred by a few hundred and want to extract rows, e.g. (0,1,0), or (1.3, 2.4, 4). I made an example: > A = rep(rep(c(0,1),4),2) > B = rep(rep(c(0,1),each=4),2) > C = rep(rep(rep(c(0,1),2),each=2),2) > X = data.frame(A,B,C,rnorm(16)) A B C [1,] 0 0 0 0.72767870 [2,] 1 0 0 -0.09673298 [3,] 0 0 1 0.51192790 [4,] 1 0 1 -1.54350441 [5,] 0 1 0

Dear Rusers, The default behavior in R when performing a regression model with missing values is to exclude any case that contains a missing value? How could i set the bahavior that R deal with missing values? e.g.: exclude cases listwise exclude cases pairwise replace with mean Thanks very much! -- Kind Regards, Zhi Jie,Zhang ,PHD Department of Epidemiology School of Public Health Fudan

Hi all: I want to fit a multinomial logit model with VGAM package, however I cannot find a way to check the residuals since residuals(my_model) and resid(my_model) does not work. Any suggestions? Thanks in advance. Alex Brito [[alternative HTML version deleted]]

Please, I??d like to store this sample matrix as a new object. How can I do this ? pulse <- c(67, 67, 68, 68, 68, 69, 69, 69, 69, 69, 70, 70, 70, 70, 71, 71, 72, 72, 73, 74) m <- NULL x <- 0 for (i in 1:5) { x <- sample(pulse,3) m <- mean(x) cat(x,m,"\n") } Thanks, Mauricio

R Help: I have read a number of tables into R with identical headings and I would now like to make a single table that has all the data appended under this single heading line. for example: t1 <- read.csv("f1",header=TRUE) t2 <- read.csv("f2",header=TRUE) all <- c(t1,t2) #all is now twice as wide as t1 or t2 with the same number of row!!!! #I need to know how

I am using R in a Microsoft Windows environment. I have a dataset called “mp1b”. I have a variable called h. h can take a value from -1 to 5. If h <1, I want to create a new dataset called mp2 that is the same as mp1b: “mp2<-mp1b” If h > 0, I want to set create a dataset mp2, where I limit the original dataset to those where mp1b$group = =h. similar to:

Dear useRs? I have the following problem! I have a function that calls one or more functions, depending on the input parameters. I am searching for the fastest way to select and execute the selected functions and return their results in a list. The number of possible functions is 10, however usually only 2 are selected (although sometimes more, even all). For examples, if I have function

Greetings all: OK, this is bugging the @#@%* out of me. I know the answer is simple and straightforward but for the life of me I cannot find it in the documentation, in the archives, or in my notes (because I know I've encountered this in the past). My problem is: I have a data frame with columns A, B, C, D, and E. A, B, and E are factors and C and D are numeric. I need a new data frame with

>From Tom: The subject is debugging a program written in the R language,under Windows. (Sorry, but I do not know either the Apple OS or *nix.) A computer program will usually not work on the first try, if only because of the risk of typos. Instead, it must be debugged. Roughly, here is the sequence: (1) One codes a program using the R language, and stores it on the hard drive, using

I have a .txt file obtained by saving a data frame in which the first four columns are factors (but represented as 1,2,3 etc). The first four lines are "Pig" "Evit" "Cu" "Litter" "Start" "Weight" "Feed" "Time" "4601" "1" "1" "1" 26.5 26.5 NA 1 "4601" "1"

Dear all, In SAS (GLM and MIXED) estimable functions (linear functions of the parameters) can be specified in the ESTIMATE and CONTRAST statements. Has anyone written a similar "utility" for use in connection with lm? Thanks in advance S?ren H?jsgaard -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello again. Now I have the following: > xx [1] "Solve(1 - R/100==(2*y)/10000,R)" > yacas(xx) expression(list(R == 100 * (1 - 2 * y/10000))) > I would like to put in a value for y and obtain R. I've tried more stuff with eval and Eval, but no luck yet. Any suggestions would be much appreciated. Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer

Hi, I am trying to index the rows of a model matrix where, due to missing values in the dataset used for the regression, the row numbers sometimes "jump" like this: > X[17:19,] (Intercept) PCINIFM 17 1 36.76471 19 1 13.23529 20 1 32.35294 So, the problem is if I select for example X[20, ], I actually get the 21th observation. I've tried to work

R CMD check gives an error when a vignette entry %\VignetteDepends{} has any extra space at the end. Maybe this is by design, but I found it confusing. Changing the function "vignetteMetaRE" in tools package could solve the problem. I.e. vignetteMetaRE <- function(tag) paste("[[:space:]]*%+[[:space:]]*\\\\Vignette", tag,"\\{([^}]*)\\}", sep = "")

Dear all, In S+, there are functions like create.ole.object call.ole.method release.ole.object for communicating with other programs which work as a COM server (on Windows). Is it possible to do something similar in R (I've studied the 'connections' facilities, but they do not seem to work). ========================================== S?ren H?jsgaard, PhD, Senior Scientist

Dear all, I have used the Scompile package (by Mathhew Calder) on previous releases of R (on Windows), and it has worked well. However, R.1.3.1 crashes when dyn.load'ing C-functions created by Scompile. >From the last R-newsletter it seems that there has been some changes to how external routines are being included in version 1.3.1. I won't claim to understand the details, but I was

Greetings! Is there any way to to lockout a certain machine identified by its ip address after a fix number of bad logon attempts? I'm aware of account policy option "bad lockout attempts" but that would block the whole account and not just the machine. Have a sunny day! Ferdinand