similar to: Bug in %in% (match)

hi Michael, am sorry for the incomplete reply . csv file data havinmg like this , >getSymbols("IBM") >weekly_data = to.weekly(week_name) >dataframe=data.frame(weekly_data) >outputfile_name="F:\\R-programs\\Outputfile.csv" >write.table(dataframe, file =outputfile_name,sep = ",",col.names =TRUE,row.names = T) > datafrom_table <-

Hi I have a called fil consisting of the following strings. > fil [1] " 102.2 639" " 104.2 224" " 105.1 1159" " 107.1 1148" " 108.1 1376" [6] " 109.2 1092" " 111.2 1238" " 112.2 349" " 113.1 1204" " 114.1 537" [11] " 115.0 303" " 116.1 490"

Hi, I have a question about estimating IRT models by using nlme, not just rasch model, but also other models. Behavior Research Methods <http://www.springerlink.com/content/1554-351x/> Volume 37, Number 2 <http://www.springerlink.com/content/1554-351x/37/2/>, 202-218, DOI: 10.3758/BF03192688 Using SAS PROC NLMIXED to fit item response theory models (2005). Ching-Fan

Hi all, Hope someone could help me. I am trying to run automatically the conversion of an Rwn file to a tex file. I am using windows 7, and cygwin. I tried to run automatically the Sweave.sh script, in its the most recent version available at R webpage: http://cran.r-project.org/contrib/extra/scripts/Sweave.sh Unfortunately, I got this error message: =========================== Raquel at

Dear R users,  I hope you can help me as you have done before.  I was wondering if I can make some of the points on my graph a different colour and symbol? Some of my soil samples are enriched and I wanted them to be shown with a red triangle and the samples which are not enriched I wanted to show with a blue square. I can't seem to find a way of doing this if there is a way?  My script:  

I'm running nut 2.0.1 with an MGE pulsar evolution 3000 with the mge-utalk (serial) driver on NetBSD 3.1 . I had a very strange occurance in my basement datacenter this early AM. the following was logged by upslog: 20080219 023452 100.0 123.0 049.0 [OL CHRG] NA 60.00 20080219 023522 100.0 122.8 049.0 [OL CHRG] NA 60.00 20080219 023552 098.0 123.0 049.0 [OL DISCHRG] NA 60.00 20080219 023622

Colleagues, I am fitting an Emax model using nls. The code is: START <- list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL <- list(maxiter=1000, warnOnly=T) #FORMULA <- as.formula(YVAR ~ EMAX - EFFECT * XVAR^GAMMA / (XVAR^GAMMA + C50^GAMMA)) ## alternate version of formula FORMULA <- as.formula(YVAR ~ EMAX - EFFECT / (1 + (C50/XVAR)^GAMMA)) FIT <-

Hi, I would like to overlap the cdf curve for observed and generated data  Here is my code: plot(cdf,main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", xaxs="i",yaxs="i",col=c("black","red"), lty=c(1,1),ylab="Cumulative probability", xlim=c(0,800),lwd=1) lines(ecdf(datobs))

Hi r-users, I would like to overlap 2 ecdf plots.  I tried this below and it gives me two plots of ecdf but just both just in black. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs)) lines(ecdf(gam_sum_gen)) Then I try to add colors etc and also the legend but fail. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i")

Hi, I would like to plot time against rainfall data (data is at the end) using xyplot. The basic code looks like this: xyplot(rainfall~time, type="a") When I do this, the graph looks ok except that the x-axis has too many values. I would just like to display the years and not the months on the x-axis. I've been fiddling around with 'scales', and read previous posts about

Hello, dear R-community. This is a question about TukeyHSD between factor combinations of a Three-Way ANOVA, which is - since it is a multi measure ANOVA - not a simple ANOVA but a Generalised Linear Mixed Model (GLMM), calculated with "lmer". > growth <- groupedData(length~meas|box_id,outer=~spec*comp*water,data=all.spec) > model <-

Hi R-users, I'm new to R and I'm trying to make a barplot combined with two lines (refering to secondary y-axis). Bars should represent the number of transfused patients by age class and sex and lines should represent the amount of blood units given in age classes. I have now successfully made a barplot and used par(new=TRUE) to plot another empty graph at the top of the barplot.

Hi: just got an APC UPS on hand. it's a strange new one, not smart-ups or back-ups. it's APC Easy UPS: https://www.schneider-electric.com/en/product-range/65883-easy-ups-on-line/ I don't know what's different. this series are all online UPS and they are cheap. it comes with serial and usb interface. with serial interface I can use apcsmart to connect. but with usb

Hi all, How can we use ks.test() to evaluate goodness of fit of mixtures of distributions? For example I have the following dataset: > x [1] 176.1 176.8 259.6 171.6 90.0 234.3 145.7 113.7 105.9 176.2 168.9 136.1 [13] 109.2 110.3 164.3 117.7 131.3 163.7 200.4 196.4 196.2 168.6 190.4 127.5 [25] 136.0 114.2 112.0 91.9 333.4 295.5 172.0 293.3 91.7 289.7 118.8 55.1 [37] 161.9 233.9 197.7

I'm trying to extend some work in the car and heplots packages that requires getting a table of multivariate means for one (or later, more) terms in an mlm object. I can do this for concrete examples, using aggregate(), but can't figure out how to generalize it. I want to return a result that has the factor-level combinations as rownames, and the means as the body of the table

Hi, I have a few questions about glht() and the interpretation of output from Tukey's in multcomp package for lme() model. The main issue is that I noticed that a plot that I produced with code letters seem to contradict the graph itself. I provide data and code below. I end with my questions. A few things about data set. "LMA.vcp" is continuous response variable.

Hello. I'm a new user of R, and I have a question regarding the use of aov and lm-functions. I'm doing a fractional factorial experiment at our production site, and I need to familiarize myself with the analysis before I conduct the experiment. I've been working my way through the examples provided at http://www.itl.nist.gov/div898/handbook/pri/section4/pri472.htm

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Has anyone used nut with a XPCC UPS, in particular one of the XVRT units? Their supplied monitoring tool "works", but I'd really prefer to use nut if possible. I've systematically gone through each of the drivers in /lib/nut and none seem to talk with the unit properly. I'll be talking with a XPCC support person tomorrow, though

Hi to everybody. Some days ago, I bought an UPS by Forza Power Technologies, model SL-1001. This Forza SL-1001 is line-interactive, has 14,4 A/h charge and is rated at 1000VA. Before, I had a TrippLite Omni 500 (or something like that), a simple contact-closure UPS, and I used nut with the simple genericupsdriver, so I'm new with this type of UPSs. NUT doesn't list Forza as a

I apologize that I've asked a similar question before, but being new to R I don't think I did a very good job of formating the question. I've included a text file since the date set is somewhat large. What I have is a huge string of numbers in a text file. The numbers are all separated by comma's and the groups are separated by a semicolon. What I would like to do is read each