similar to: XYPLOT - GROUPING WITH TWO CATEGORICAL VARIABLES

Hi everybody!! I don`t know if it is correct to make a question about R and Tinn-R in this mailing-list, but I need some help with a problem for which I don't find a solution anywhere. I have Tinn-R 2.3.5.2 in my notebook and I always used it with no problems, until a few days ago when this signs appears everytime I start the programme: "Key Violation. C:/Documents and

hi, people How can we compare two probit models brought out from the same data? Let me use the example used in "An Introduction to R". "Consider a small, artificial example, from Silvey (1970). On the Aegean island of Kalythos the male inhabitants suffer from a congenital eye disease, the effects of which become more marked with increasing age. Samples of islander males

Dear all, I would like to produce a map with information about the patenting activity in German districts, by coloring districts with different degrees of patenting activity in different colors. I work with the packages maptools, maps and spdep. The map data is read from an external .shp file (+ the corresponding .shx and .dbf files). Plotting a map with the IDs or the patenting indicator itself

Hello! I want to display four xyplots in the same window (sorry if I don't use the propers words, but my English is rude) I know that if I use plot(y~x) then x=(1:100) y=(1:100) par(mfrow=c(2,2)) plot(y~x) plot(y~x) plot(y~x) plot(y~x) would give what I want, but this is not the case because I need the functions of xyplot. So I did the following: plot.new()

I'm have a (minor) problem and a question. Problem: The following code creates 5 clusters of dots of different colors. However, I need the second call outside the data.frame call to get the colors to change for some reason. I assume there's an error in the data.frame() call, but I don't know what. Question: Is there is a way to cause the resulting graphic to rotate on it's

hello, I am using arima to evaluate a time series regression model. I am using categorical variables such as day of week(Su-Sa), month(Jan-Dec), and holiday status (1/0) as my independent variables. There is evidence of multiplicative interaction between holiday status and weekday and I would like to add an interaction term to my arima model, but I am not sure about best way to go about doing it.

Hello List Inhabitants: I don?t know what this operation is called or if there is a function that does it automatically, hence I seek your help! If I divide a large data collection tasks between two students, and I have a master list of samples, but one student records some of the values, and the other student the remainder, I need to get the two sets of student measurements into one column.

Hello Osama, and Hisham, At 1330GMT there was some malicious activity coming from your network IP 37.75.210.90. Please act accordingly. Things that may be of use "972599779558" N.

Hi, I just wanted to plot a boxplot with a nice curve going through it, I thought this would be a simple task but for some reason I can't get the two graphs on the same page accurately. Enclosed is the code showing the two plots seperately and together. I would have thought it should work if I could use boxplot() then plot() overlayed but it won't allow the command add=TRUE (which

i am trying to replicate the following graph using xyplot : attach(x) plot ( jitter(type), mortality, pch=16, xlim = c(0.25, 3.75)) lines ( c(1-0.375,1.375) , c ( median(mortality[type==1]), median(mortality[type==1])), lwd=5,col=2) lines ( c(2-0.375,2.375) , c ( median(mortality[type==2]), median(mortality[type==2])), lwd=5,col=2) lines ( c(3-0.375,3.375) , c ( median(mortality[type==3]),

Hello R Community, I have the following design question. I have a data set that looks like this (shortened for the sake of example). Gender Age M 70 F 65 M 70 Each row represents a person with an age/gender combination. We could put this data into a data frame. Now, I would like to do some actuarial analysis on this data set. To do so, I need to create and store

Dear R helpers,     Yesterday I had raised following query which was addressed by Mr Ellison. The query and the wonderful solution as provided by Mr. Ellison are as given below.    ## PROBLEM   I am calculating the 'Yield to Maturity' for the Bond with following characteristics.    Its a $1000 face value, 3 year bond with 10% annual coupon and is priced at 101. The yield to maturity can be

Estimados Queremos con el paquete FactoMineR hacer este tipo de tabla de mortalidad que lea los datos desde de una tabla csv Realizamos lo que viene en la ayuda y es muy interesante, sin embargo cuando mandamos a leer desde la tabla csv original de los autores no hace el análisis porque algo falta y no nos percatamos de que es. Adjunto tabla original Saludos cordiales #ESTO ES LO QUE

Hola. Como te dijo Carlos, el problema está en los nombres de las columnas y en los nombres de las filas. Cuando hice la importación (con dd<-read.csv('mortality.csv'), tuve problemas con las filas de nombre: - Malignant tumour of the larynx trachea bronchus and lungs - Malignant tumour of the lip pharynx and mouth - Other endocrinological metabolic and nutritional conditions

Dear R Users: I have the individual mortality rate and 95% CI of 100 hospitals, how to do the plot with the individual hospital in the Yaxis, and the mortality rate and 95% CI in the Xais and a overall mean as a reference line? Thanks and regards, Xin [[alternative HTML version deleted]]

Hi i am doing an environmental research The equation is as follow: gam(y1 ~ x1 + s(x2) + s(x3) + s(x4), family = gaussian, fit = true) I would like to obtain the beta coefficient and 95CI of x4 (or s(x4)), what should I do? Thanks, Lung -- View this message in context: http://r.789695.n4.nabble.com/mgcv-beta-coefficient-and-95-CI-tp3320491p3320491.html Sent from the R help mailing list

Hi all, I am not sure how to handle interactions with categorical predictors in the GAM models. For example what is the different between these bellow two models. Tests are indicating that they are different but their predictions are essentially the same. Thanks a bunch, > gam.1 <- gam(mortality.under.2~ maternal_age_c+ I(maternal_age_c^2)+ + s(birth_year,by=wealth) + +

Dear all, I'm trying reproduce an example of bioassays Yielding Concentration-Mortality Data particularly control - adjustment model from book Bioassay of Entomopathogenic Microbes and Nematodes chapter 7 with R. I used glm with family=binomial and link=probit, but I do not know how to implement parameter gamma (control mortality - mortality of the untreated control insect in this exaple)

Hello, all. I have co-opted a number of functions that can be used to plot the hazard/survival functions and associated density distribution for a Gompertz mortality model, given known parameters. The Gompertz hazard model has been shown to fit relatively well to the human adult lifespan. For example, if I wanted to plot the hazard (i.e., mortality) functions: pop1 <- function (t) {

I created a 2 D plan: > plot(range(-2.5,0.95),range(0.00,1.00),type="n",axes=TRUE) I made a projection of points with their coordonates (X,Y) in that plan > fullpoints=read.csv2("fullpoints.csv",h=T) > plot(fullpoints) The points are listed in that .csv file it is organized this way: 0,48875 0,142857143 0,409 0,142857143 0,45611 0,25 0,49833 0,222222222 0,61158