similar to: R 1.7.0 startup error: .addMethodFrom...

xyplot doesn't seem to want to label my x-axis with dates but instead puts the day-number for each date. begdate is the number of days since January 1, 1960 and was initially created by library(date) ... polls$begdate<-mdy.date(begmm,begdd,begyy) I create a new dataframe (pollstack) which includes begdate. In the process begdate seems to lose its date attribute so I redo it as: >

What is a simple way to convert a factor into a matrix of dummy variables? fm<-lm(y~f) where f is a factor takes care of this in the estimation. I'd like to save the result of expanding f into a matrix for later use. Thanks. Charles -- Charles H. Franklin Professor, Political Science University of Wisconsin, Madison franklin at polisci.wisc.edu chfrankl at wisc.edu

A newbie question: I'm trying to decide whether to run a maximum likelihood estimation in R or Stata and am wondering if the R mle routine is reasonably robust. I'm fairly certain that, with this data, in Stata I would get a lot of complaints about non-concave functions and unproductive steps attempted, but would eventually have a successful ML estimate. I believe that, with the

Is there any way to echo comments from an R source file into an SWeave->LaTeX document? Example: # Npop is population total # Npoph0..Npoph2 are stratum totals # Npoph is vector of stratum totals Npop<-sum(to2000) Npoph0<-sum(to2000[bg==0]) Npoph1<-sum(to2000[bg==1]) Npoph2<-sum(to2000[bg==2]) Npoph<-c(Npoph0,Npoph1,Npoph2) In the final LaTeX document, I'd like the

II am translating some gauss code into R, and gauss has a matrix product function called the horizontal direct product (*~), which is some sort of variant on the Kronecker product. For example if x is 2x2 and y is 2x2 the horizontal direct product, z, of x and y is defined (in the Gauss manual) as: row 1 = x11*y11 x11*y12 x12*y11 x12*y12 row 2 = x21*y21 x21*y22 x22*y21 x22*y22 Or in R

Hi there Within sampleSelection, I'm trying to calculate the marginal effects for variables that are present in both the selection and outcome models. For example, age might have a positive effect on probability of selection, but then a negative effect on the outcome variable. i.e. Model<-selection(participation~age, frequency~age, ...) Documentation elsewhere describes one method for

Hello, I'm using optim to program a set of mle regression procedures for non-normal disturbances. This is for teaching and expository purposes only. I've successfully programmed the normal, generalized gamma, gamma, weibull, exponential, and lognormal regression functions. And optim returns reasonable answers for all of these compared with the identical optimization problems in STATA and

Hi, With the help of Dirk Eddelbuettel, I was able to fix a few little problems that prevented the rimage (0.5-7) package from working under (at least) R 2.8 or 2.9 and g++ 4.2 or 4.3. I have uploaded a new release, 0.5-8, that fixes these problems. It have attempted to contact the package's maintainer, Tomomi TAKASHINA <t.takashina at computer.org>, to no avail. Kurt Hornik

FYI ---------- Forwarded message ---------- From: Bear Braumoeller <braumoeller.1 at polisci.osu.edu> Date: 2009/8/4 Subject: Re: Package boolean on CRAN To: "cgb at datanalytics.com" <cgb at datanalytics.com> Cc: Ronggui Huang <ronggui.huang at gmail.com>, Ben Goodrich <goodrich at fas.harvard.edu> Gentlemen, Thanks for your enquiries. ?Boolean is under

My problem is I have data with both categorial and numerical data, currently only the categorical number contains missing data, was wondering do I make a new dataframe containing only the categorical columns? How would you use Latent Class Model specifically poLCA to impute the missing data? http://www.sscnet.ucla.edu/polisci/faculty/lewis/pdf/poLCA-JSS-final.pdf The reason why I chose not to

Hi, With the help of Dirk Eddelbuettel, I was able to fix a few little problems that prevented the rimage (0.5-7) package from working under (at least) R 2.8 or 2.9 and g++ 4.2 or 4.3. I have uploaded a new release, 0.5-8, that fixes these problems. It have attempted to contact the package's maintainer, Tomomi TAKASHINA <t.takashina at computer.org>, to no avail. Kurt Hornik

I am trying to use the optimizing function genoud() with the snow package on a couple of i686 machines running Redhat Linux WS4 . I don't know anything about PVM or MPI, so I just followed the directions in snow and rgenoud for the simplest method and started a socket cluster. My function fn for genoud involves updating an input file for a separate numerical model with the latest parameter

I would rather use cleveland dot plots than bar charts to display my study results. I have not been able to find (or figure out) an R package that is capable of producing the publication quality dot charts Im looking for. I have either not been able to get error bars (lattice), cannot order the data display properly (latticeExtra), or cannot make adjustments to axes. Does anyone have a quick

Hello, I´ve encountered a couple hurdles after installing Ubuntu for the first time and working with R with RKWard. The Keyboard: I have to hit the quote button twice to get a quote or apostrophe! I tried adjusting keyboard settings but nothing changed. And the RKWard GUI won´t recognize quotes from my machines even though they show up on the RKWard text editor. On the other hand, RKWard will

Hello, What is your opinion on Ubuntu. I´ve encountered a couple hurdles here. The Keyboard: I have to hit the quote button twice to get a quote or apostrophe! I tried adjusting keyboard settings but nothing changed. And the RKWard GUI won´t recognize quotes from my machines even though they show up on the RKWard text editor. On the other hand, RKWard will quoted material whien I copy and

The Wine development release 1.3.36 is now available. What's new in this release (see below for details): - UDisks backend for the dynamic device support. - Polygons and geometric pens implemented in the DIB engine. - JScript now compiles entire functions. - Improved support for vertical fonts. - A few more scripts in UniScribe. - A bunch of ctype functions implemented in MSVCP.

I have search the archive and I could not find what I need so I will try to ask the question here. I read a table in (read.table) a <- read.table(.....) The table has column names like DayOfYear, Quantity, and Category. The values in the row for Category are strings (characters). I want to get all of the rows grouped by Category. The number of unique category names could be around 50. Say

Hi, I'm trying to get zapras do work, I had downloaded the pppd-source and the 2 patches. I succefull compiled and install the patched version of pppd, but got this error in message-log Jul 15 11:43:32 voip1 pppd[9296]: In file /etc/ppp/filters: unrecognized option 'active-filter' Jul 15 11:43:57 voip1 pppd[9299]: Plugin zaptel.so loaded. Jul 15 11:43:57 voip1 pppd[9299]: Zaptel

Hi, I have simple question about Skype. What are the methods of selecting packets which belongs to Skype?? I know about 7layer but I don''t belive that is only way. Is 7layer realy good and stable solution for routers which must handle more than 1000 users ? Thanks in advance Pozdrawiam Szymon Turkiewicz

Hello, I use HTB with linux-2.4.20-ac1 and iptables-1.2.7a (pathced for IMQ). I configured one HTB qdisc on an IMQ device. This qdisc has some classes (about 5) and all of the classes has qdisc pfifo. I sometimes got out from my linux (while I was working) and I found in logs: Dec 4 11:49:27 core kernel: HTB init, kernel part version 3.7 Dec 4 11:49:28 core kernel: HTB: quantum of class