similar to: Thank you

Forgive my resending this post. To data I have received only one response (thank you Bert Gunter), and I still do not have an answer to my question. Respectfully, John Windows XP R 2.12.1 contrast package. I am trying to understand how to create contrasts for a model that contatains an interaction. I can get contrasts to work for a model without interaction, but not after adding the

I am running my code in a loop and it does not work but when I run it outside the loop I get the values I want. n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta <- 0 ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit))

I am new to R and I am trying to do a monte carlo simulation where I generate data and interject error then test various cut points; however, my output was garbage (at x equal zero, I did not get .50) I am basically testing the performance of classifiers. Here is the code: n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta

R 2.5 windows XP I am getting an error from glm() that I don't understand. Any help or suggestions would be appreciated. N.B. 1<=AAMTCAREJ<=327900 > summary(data$AAMTCAREJ) Min. 1st Qu. Median Mean 3rd Qu. Max. 1.0 404.3 1430.0 6567.0 5457.0 327900.0 > fitglm<-glm(AAMTCAREJ~sexcat+H_AGE+SmokeCat+InsuranceCat+MedicadeCat+ +

On Thu, Jul 30, 2009 at 05:44:05PM +0200, Andreas Neustifter wrote: > Well it does look not too good the text is getting to small in this > case, what do you think of this 128x128 version? It looks like an sprite coming from a Megaman game :). My humble opinion is that while the dragon looks nice when large, when small it just feels childish and unprofessional. My 2 cents. -- Felipe.

library(ROCR) n <- 1000 fitglm <- function(iteration,intercept,sigma,tau,beta){ x <- rnorm(n,0,sigma) ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit)) *int*<-coef(model)[1] *slope*<-coef(model)[2] # when add error you are suppose to get slightly bias slope. However when I change

Hello, I would like to know if dovecot/pigeonhole provide a way to do sieve redirection with sender rewriting scheme (srs)? I apologize in advance if that e-mail address isn't meant for assistance, but I've asked on the official irc channel on freenode and was told to direct my query to dovecot at dovecot.org, I hope you won't mind. Thanking you in advance Best Regards

At useR 2006 I mentioned that it would be nice to have a way to specify binomial links that involved free parameters and described some experience with a Gosset link involving a free degrees of freedom parameter, and a Tukey-lambda link with two free parameters. My implementation of this involved some rather kludgey modifications of binomial, make.link and glm that (essentially) added a

Hi, I am a brand new user of R and I have a really stupid problem: I am having troubles applying fft on 2D matrices. I import a data file of 256*256 elements with read.table (it is actually a grey scale image, the corresponding data file being generated with Matlab), and I find that the fft can be applied on one single column of data, but not on one row. (Error message: Error in fft(z, inverse)

" Does anyone know if there is a program for linux like ices that can dump live audio from an input to an mp3 file? Or, is there any software for windows that can capture the mic input and line in input and encode them seperately but both live? Thanks, Andy Carlson IT Manager City of Princeton www.princeton-il.com Office (815) 875-2631 ext. 253 Cell (815) 719-0023

Hi. Can anyone suggest a simple way to obtain in R a list of vector coordinates of the following form? The code below is Mathematica. In[5]:= Flatten[Table[{i,j,k},{i,3},{j,4},{k,5}], 2] Out[5]= {{1,1,1},{1,1,2},{1,1,3},{1,1,4},{1,1,5},{1,2,1},{1,2,2},{1,2,3},{1 ,2,4},{1,2, 5},{1,3,1},{1,3,2},{1,3,3},{1,3,4},{1,3,5},{1,4,1},{1,4,2},{1,4,3}, {1,4,

Hello everybody, I'm working on implementing the speex codec on ARM966E @ 96Mhz. For the encoding I found that 20ms of sample voice takes 13ms to be encoded. The goal of my project is to built an VoIP application based on ARM966E but with the actually performance it's not possible :( My question is what is approximately the speex encoding time with narrowband mode on ARM9 (ARM926 or

Hello, One year ago a bug was reported regarding pigeonhole vacation plugin https://www.mail-archive.com/dovecot at dovecot.org/msg78588.html a mail sent to User at domain.com instead of user at domain.com will generate a " vacation action: discarding vacation response for implicitly delivered message; no known (envelope) recipient address found in message headers (recipient=<user at

Here is an even more awkward property of tapply with factors. What we want to do is to determine the modal value of the first factor given the level of the second factor. We do that by applying table in a tapply call then finding the maximum count then ... Turns out that table gives unexpected results in tapply. > table(Machines$Machine) A B C 18 18 18 > tapply( Machines$Machine,

Hello! I just uploaded a photo on zroslina's DailyFlog page that I want you to see. Please come and see: http://www.mydailyflog.com/go/invite_register/zroslina/56282492&stc=18 Thanks! Roslina Zakaria ======================================== Got a digital camera? MyDailyFlog is a personal photo-blogging space where you can easily post your latest and greatest photos, and share them

HI, I would like to extract data in a specific way.  For example, the rainfall data 0,0,1.5,0,0, 3,1,2.5,0,0,0,0, 2.3,0,0,0, 2.1,1.4,0,0,0, 3,2,1,0,0,0... data_1: 1.5, 2.3               ( a single nonzero data between zeros data) data_2: 3.1, 2.5, 2.1,1.4   ( two nonzero data between zeros data) data_3: 3,1,2.5, 3,2,1       ( three nonzero data between zeros data) Thank you so much for any

Dear r-users, I have dailly rainfall data from year 1971 to 2000. I use aggregate to form monthly rainfall data.  What I don't understand is that the data for the year 2000 become on the top, instead of year 1971.  Here are some codes and output: agg_dt1     <- aggregate(x=dt1[,4],by=dt1[,c(1,2)],FUN=sum) > head(agg_dt1,20); tail(agg_dt1,20)    Tahun Bulan     x 1      0     1 398.6

One Question: if I make a lattice-plot (beware: under development!) with a special layout e.g. c(2,2,4): p1 <- xyplot(y~x|id,data=data,layout=c(2,2,3)) is there a way to print.trellis e.g. the second page of the trellis object "p1"? print.trellis(p1,page=2) #!wrong! I want this to automatically dev.print the pages (no, I don't want multiple postscript pages). If this

Hi R-users, I would like to use the copula package.? I? the package plus the mvtnorm and try to run the example given, but I got the following message: install.packages(repos=NULL,pkgs="c:\\Tinn-R\\copula_0.8-3.zip") norm.cop <- normalCopula(c(0.5, 0.6, 0.7), dim = 3, dispstr = "un") t.cop <- tCopula(c(0.5, 0.3), dim = 3, dispstr = "toep", df = 2, df.fixed =

Hello, I'm trying to use the tool virsh of libvirt. I want to manage via CLI a vm on a vmware hypervisor which is hosted on a remote host (Is that what you call the node??). I read a lot of thing on your site but I really don't understand how to do something... I saw that: >virsh edit $your-vm-name I think "$your-vm-name" is the xml file you're talking about? Is it