similar to: Plotting series with no data in xyplot

Good Day All, I have a negative binomial model which I have developed using the MASS library. I now would like to develop some predictions from it. Running the predict.glm (stats library) using type="response" gives me a non-integer value which was rather puzzling. I would like to confirm that this is actually the mean predicted value of the probability mass function as opposed

Good Day, I have several ASCII data files that I would like to import into R. They all have a SAS import file which is used to bring the data into SAS and I am hoping to use this to bring the data into R. There are lots of variables involved and the ASCII data file is 2308 columns long so I would certainly prefer to figure out a smart way of converting the data to R. The ASCII data is a

Good Day All, I have a negative binomial model that I created using the function glm.nb() with the MASS library and I am performing a cross-validation using the function cv.glm() from the boot library. I am really interested in determining the performance of this model so I can have confidence (or not) when it might be applied elsewhere If I understand the cv.glm() procedure correctly, the

Hi All, Could someone help me decode what this error means ? > BIC(nb.80) Error in log(attr(object, "nobs")) : Non-numeric argument to mathematical function > BTW, nb.80 is a negative binomial glm model created using the MASS library with the call at the bottom of the message In the hopes of trying to figure this out I tried the following workaround but it did not work

Hi All, I am working on some ordinal logistic regresssions using LRM in the Design package. My response variable has three categories (1,2,3) and after using the creating my model and using a call to predict some values and I wanted to use a simple .5 cut-off to classify my probabilities into the categories. I had two questions: a) first, I am having trouble directly accessing the

Hola, Una opción más adecuada que representar max, min y media podría ser utilizar un boxplot y representar la distribución completa de las temperaturas que tienes en el día... Si te convence la idea, aquí tienes un ejemplo con el código disponible (bajo el punto A simple boxplot using ggplot2) http://www.sr.bham.ac.uk/~ajrs/R/r-gallery.html Saludos, Carlos Ortega www.qualityexcellence.es El

Comparto totalmente la sugerencia de Carlos Ortega. Si aún así quieres gráficos de velas o gráficos OHLC el mejor paquete que puedes usar es uno de temas financieros: http://www.quantmod.com/ Ten en cuenta que en los mercados bursátiles tiene algún sentido usar ese tipo de gráfico ya que en la apertura y mayormente en el cierre diario se concentran operaciones. Me imagino que esto no es así con

Dear All, I'm having trouble ploting multiple xYplots (Hmisc) within the same pg although I've done it before w/ xyplots. I've produced 4 similar plots e.g. rna.h<-xYplot(Cbind(RNA,Lo,Up)~HO|factor(MO,labels=c("April", "May","June")),data=diel.data,method="bars",type="b", keys="lines",pch=16,col=1,cap=.04,

Dear All, I'm having trouble ploting multiple xYplots (Hmisc) within the same pg although I've done it before w/ xyplots. I've produced 4 similar plots e.g. rna.h<-xYplot(Cbind(RNA,Lo,Up)~HO|factor(MO,labels=c("April", "May","June")),data=diel.data,method="bars",type="b",

The following becomes more and more a topic for an R mailing list, possibly R-devel. R for Windows crashes are important and should be discussed and debugged with the experts -- which are not reading ESS-help AFAIK. Martin >>>>> "ColCun" == Cunningham, Colin A <colin.a.cunningham@intel.com> writes: ColCun> In my very brief exposure to R, I have

This message was rejected by the mail hub at The University of Birmingham (bham.ac.uk) because it has an apparently executable attachment "submited.pif". The email rejected was not delivered to the following recipients: M.Spann at bham.ac.uk. Executable attachments are not being accepted at The University if Birmingham because they have been used by recent viruses such as that

I am trying to drop rows of a dataframe based on values of the column PID, but my strategy is not working. I hope someoen can tell me what I am doing incorrectly. # Values of PID column > jdata[,"PID"] [1] 16608 16613 16355 16378 16371 16280 16211 16169 16025 11595 15883 15682 15617 15615 15212 14862 16539 [18] 12063 16755 16720 16400 16257 16209 16200 16144 11598 13594 15419 15589

Dear R community, Since upgrading to R v.2.6.1 and re-installing package Hmisc (binary for Mac OS X v.3.4-3), I have been getting a error when trying to make xYplots: >plotcv<-c(34.88, 41.51, 45.81, 51.05, 51.66) >plotcv.se<-c(2.406551, 3.071291, 4.331407, 3.213873, 4.838150) >month<-c(6, 7, 8, 9, 10) >library(Hmisc) >xYplot(Cbind(plotcv, plotcv + plotcv.se,

Dear all, Compiling R-0.63 on my Alpha machine immediately encounters the following message. Your help is appreciated. Make: Don't know how to make ./src/*/*.[chfy]. Stop. Regards, -- Xiaoming Cai School of Geography .&______~*@*~______&. University of Birmingham

I tried to compile R-0.63.1 to R-0.63.3 on my DEC Alpha with the operating system of UNIX 4.0. All have the same problem after R.binary has been successfully compiled. The last few lines of the message are shown as below: ... make[4]: Entering directory `/home/ming/R-0.63.3/src/library/modreg/src' f77 -c bsplvd.f -o bsplvd.o f77 -c bvalue.f -o bvalue.o f77 -c bvalus.f -o bvalus.o f77

Hi all, I'm searching for a way to exclude outliers from my dataset while making xyplots. While plotting using pairs(), I exclude specific row in my data frame and save the settings as a variable which I later include as an argument: # Discard outliers and save settings as idx idx=with(fieldTrial0809, which(Pro>0 & Pro<0.95 & Fum>0 & Fum<0.4 & Mal>0.1 &

Dear all, Compiling R-0.64.0 on DEC osf4.0 has the following error message (earlier version was OK). cc -ieee_with_inexact -g -I../include -I../../src/include -c pnf.c -o pnf.o cc -ieee_with_inexact -g -I../include -I../../src/include -c pnt.c -o pnt.o cc: Error: pnt.c, line 83: In this statement, "1021" is not an lvalue, but occurs in a context that requires one. if (df >

I have tried compiling R on a SGI 6.0 and got the following: cc -o ../../bin/R.X11 apply.o arithmetic.o array.o attrib.o bind.o builtin.o character.o coerce.o colors.o complex.o context.o cov.o cum.o debug.o devPS.o devPicTeX.o deparse.o deriv.o devices.o dotcode.o dstruct.o duplicate.o envir.o errors.o eval.o format.o fourier.o g_her_glyph.o g_fontdb.o g_cntrlify.o g_alab_her.o gram.o

I have tried compiling R on a dec alpha running digital unix 4.0 and got the following: f77 -shared -o ts.so PPsum.o burg.o eureka.o filter.o pacf.o starma.o stl.o carray.o mburg.o myw.o qr.o -lUfor -lfor -lFutil -lm -lots -lm ld: Warning: Unresolved: R_alloc R_NaReal R_IsNA R_IsNaNorNA R_chk_calloc R_chk_free Rf_error vmaxget vmaxset dqrdc2_ dqrcf_ mkdir ../../../../library/ts/libs gmake[4]:

Actually I'm using the following versions: smbd -V output: "Version 2.2.1a" nmbd -V output: "Version 2.2.1a" Stefano -----Messaggio originale----- Da: Mark Cooke [mailto:mpc@star.sr.bham.ac.uk] Inviato: gioved? 31 luglio 2003 12.14 A: Stefano Del Negro Oggetto: Re: [Samba] smbtar Samba not liking the accented 'i' character ? Just a guess, Warning - don't