similar to: feedback/question on function update()

Dear R developers, I have noticed a discrepancy between the coefficients returned by R's glm() for logistic regression and SAS's PROC LOGISTIC. I am using dist = binomial and link = logit for both R and SAS. I believe R uses IRLS whereas SAS uses Fisher's scoring, but the difference is something like 100 SE on the intercept. What accounts for such a huge difference? Thank you for

Hi all, So let's say I have a matrix, mdat and I only know the index number. How do I retrieve the column and row names? For example, > mdat <- matrix(c(1,2,3, 11,12,13), nrow = 2, ncol=3, byrow=TRUE, dimnames = list(c("row1", "row2"), c("C.1", "C.2", "C.3"))) > mdat[4] [1] 12 >

Dear R-help team, I am getting addicted to using R but keep on getting many challenges on the way especially on data management (data cleaning). I have been wanting to drop all the rows if there values are `NA' or have specific values like 1 or 2 or 3. mdat <- matrix(1:21, nrow = 7, ncol=3, byrow=TRUE, dimnames = list(c("row1",

Hi everyone. Long time listener, first-time caller here. I have a data set that's been melted with the excellent reshape package, but I can't seem to cast it the way I need to. Here's the melted data's structure: > str(mdat) 'data.frame': 6978 obs. of 4 variables: $ VehType : Factor w/ 2 levels "Car","Truck": 1 1 2 1 1 2 1 1 1 1 ... $ Year :

Hello there, I have a question regarding bar plots. I am trying to plot the data from the following matrix as a barplot - # input data mdat <- matrix(c(0.1,0.9,0.9,0.1,0.5,0.5,0.45,1-0.45,0.6,0.4,0.8,0.2), nrow = 6, ncol=2, byrow=TRUE, + dimnames = list(c("Mon", "Mon", "Tues", "Tues", "Thurs", "Friday"), +

Dear R help, Apologies for the less than informative subject line. I will do my best to describe my problem. Consider the following matrix: mdat <- matrix(c(1,0,1,1,1,0), nrow = 2, ncol=3, byrow=TRUE, dimnames = list(c("T1", "T2"), c("sp.1", "sp.2", "sp.3"))) mdat In my actual data I have time

Hi all R-users, If I have a matrix with numeric elements as follows, the values are rounded when I try to refer to a specifici element using [], the value is rounded. The same thing happens if the matrix is read from a file, the values are stored to the correct precision but then when I try to refer to a specific element (such as using [], it is rounded. How do I avoid this rounding?

Hello R users (and Hadley) I have another question about ggplot2 :-) (version 0.8) `dat` <- structure(list(D = c("a", "b", "c", "d"), G = c(1.51520888871520, 1.88812208268440, -6.60521862, 0.55968739), E = c(1.38888592256404, 1.39366168665589, 1.22509259382058, 1.36617701059296), I = c(6.92634958902857, 6.94416045215158, 13.2179488828556,

Hi, I was trying to make another legend for the rug plot. Sample code: library(ggplo2) ids <- paste('id_',1:3,sep='') before <- sample(9) after <- sample(1:10,9) dat <- as.matrix(cbind(before,after)) rownames(dat) <- rep(ids,3) position <- c(rep(10,3),rep(13,3),rep(19,3)) mdat <- cbind(melt(dat),position) ggplot(mdat, aes(position, value)) +

Hello, It looks to me as if the ! logical operator cannot be called when subsetting a data.frame. In the example below, the value column has two factor levels (but my typical datasets have more), and what I am trying to do is to exclude all lines for which the "value" is different from "A". I have got a syntax error message everytime I try to use the subset() function.

Hi R-experts, I try to find a way to transfer a matrix to a data.frame that is used as input of aov. can you give me advice for that? >mdat <- matrix(c(1,2,3, 11,12,13), nrow = 2, ncol=3, byrow=TRUE, dimnames = list(c("row1", "row2"), c("Col1", "Col2", "Col3"))) >mdat      Col1 Col2 Col3 row1    1    2    3 row2   11   12   13 ===>

Hi, I have been recently experimenting with the lattice package, which I must admit is just great! However, I'm sort of stuck in modifying certain parameters; Would appreciate some pointers on a couple of things: 1) Is it possible to change the font of the labels (say to computer modern) -- either in the Windows output or thru' Sweave (generating EPS/PDF)? 2) As you will notice, there

Hi, Can you please help me with this please? What I am trying to do is call a vector from R function and used in the new function So I create 4 functions with these arguments M11 <- function(TrainData,TestData,mdat,nsam) { ls <- list() I have few statments one of them is vectx <- c(,1,2,3,4,5,6,6) vectz <- c(12,34,5,6,78,9,90) and then................ ls(vectx=vtecx,vectz=vectz)

Hello Xen-Users, I am currently provisioning old Linux guests as HVM on a server, and giving clients access through VNC. The current setup uses the default Xen framebuffer as a VNC server (all configuration files have vnc=1 and a few other settings), and my clients are using vinagre (a.k.a. Remote Desktop Viewer) on Arch as their VNC client. Problem: Keyboard input is not consistently handled

Hello All, I am trying to color five subsets of points in a series of vertical strip charts by quintiles using the quantile function, but I have not been able to do this successfully, nor have I been able to find an example of this online. I would like to have the "E1", "E2", and "E3" columns plotted on a single plot (containing the three strip charts). I am

Bad scaling will waste a lot of everyone's time. I put the data in a data frame mdat, then library(nlmrt) mdat<-read.csv("muzzamil.csv", header=T) fmn <- nlxb(y~a * (x^b), data=mdat, start=c(a=1,b=1), trace=T) fm <- nls(y~a * (x^b), data=mdat, start=c(a=1,b=1), trace=T) fmn2 <- nlxb(y~a2 * ((x-1979)^b2), data=mdat, start=c(a2=1,b2=1), trace=T) fm2 <- nls(y~a2 *

dear all, i am anewcomer. i have a set of paired data, the values are great different from each other. i want to show them in histogram. first, how to draw each bar for corresponding paired data side by side; second, how to set the scale of y axis? one is up to 100-fold to the other. ifdoubleslash is used to omit some coordinates, the histogram will looks better. is there anyone that can help

Dear R-developers, I would like to suggest a 'method' slot for format.ftable() (see an adjusted 'format.ftable()' below, taken from the source of R-2.15.2). At the moment, format.ftable() contains several empty cells due to the way the row and column labels are printed. This creates problems (= unwanted empty columns/rows) when converting an ftable to a LaTeX table; see an

I just downloaded Debian Squeeze and installed Xen from packages, following the instructions here: http://wiki.debian.org/Xen I am able to get PV guests up and running easily, but HVM guests will not even boot. The HVM I am currently trying to run is Debian sarge, which I installed into a disk image and configured using qemu. I am running on an HP Proliant ML150G6, which has Intel-VT hardware

Hi, I have a problem with understanding what the function cast() from the package reshape is doing. In the example below I have a 2x2x2 array which I first melt and then cast to get the averages over the field 'strain' for every combination of the fields 'treatement' and 'gene': ------ > mdat <- melt(array(rnorm(8), dim=c(2,2,2),