similar to: Re: Universal legend in plot

Dear R-users: I asked a question on how I can have a universal legend in a plot and received the following result. I tried using "layout" but I can't seem to work on the "empty" plot (where I have to have the legend). I tried "oma" but I couldn't improve the quality of the plot, and that I didn't know how to specify all the line types using the

Dear R users: I am not entirely convinced that clogit gives me the correct result when I use pspline() and maybe you could help correct me here. When I add a constant to my covariate I expect only the intercept to change, but not the coefficients. This is true (in clogit) when I assume a linear in the logit model, but the same does not happen when I use pspline(). If I did something similar

Dear R users, I would like to sample from a mixture distribution p1*f(x1)+p2*f(x2). I usually sample variates from both distributions and weight them with their respective probabilities, but someone told me that was wrong. What is the correct way? Vumani

Dear R users: I have the following matrix. 0 1 1 0 1 0 2 1 0 3 0 0 I would like to spread the matrix such that whenever the row sum is greater than 1 the row is repeated the number of times given by the row sum. Furthermore I would like to split the following cases: 0 1 1 such that it map to the following matrix 0 1 0 0 0 1 such that each row adds up to 1. I have no problems with cases

Dear R users; I am trying to write a small program which reads in a data set, and selects observations from certain years before the analysis. I have a problem including the selection criteria in the header of the program. Here is the problem; dataFIT<-function(MODEL, MARGINS, yearConsidered){ library(foreign CovPaper<-read.spss("C:/Data/data.sav") NewData <-

Dear R users, I have a matrix with 2 columns with the variables: "daughter index", "mother index". I would like to draw a network digraph using this data, where each daughter is connected to a mother and between the connections inlcude a circle with the information on the indices ("daughter index", "mother index"): i.e. something similar to graphs

it is easy to fit a distribution using fitdistr poisdata <- rpois(n = 100, lambda = 2) poismle <- fitdistr(poisdata, "Poisson") poismle but i would like to know whether its possible to get an identical result using glm. I use poistab <- data.frame(table(poisdata)) colnames(poistab) <- c("width","freq"); poistab[,"width"] <-

Dear R users: Is there any R analogue for the S+ function "na.gam.replace". I would like to make an interaction of a categorical and smooth continuous covariate. Thanks. Vumani Dlamini Central Statistical Office Swaziland

Dear R users, Is it possible to force a hard return in an expression. I tried including one for a legend using an expression and it didn't work. Here is my code, legend(12, 0.10, c("Fitted density", as.expression(substitute(paste("Weibull (", alpha==shapeU," , ", beta==scaleU,") Mean = ",meanU," days"), list(shapeU =

Dear all; I am trying to use HMISC imputation function to perform multiple imputations on my data and I keep on getting errors for the code given in the help files. When using "aregImpute" the error is; >f <- aregImpute(~y + x1 + x2 + x3, n.impute=100) Loading required package: acepack Iteration:1 Error in .Fortran("wclosepw", as.double(w), as.double(x),

Dear R users; I would like to convert a series of vectors to matrices in the following way; (2,1,1) to a matrix 1 0 0 1 0 0 0 1 0 0 0 1 The idea is that the column sum of the matrix should be equal to the elements of the vector. Thanks. Vumani

Dear R users, I'm having trouble using aggregate() and would greatly appreciate your advice. I am using R 1.5.1 on Windows 2000. I want to call my function in the following way extract.data(x=dat[, "Age", "Year"]) where extract.data() uses aggregate() to count the number of cases for each combination of "Age" and "Year". I've defined

One of the CRAN mirror links for Australia is not correct, and I'm wondering how to have it updated. The incorrect address is http://cran.mirror.aarnet.edu.au/ and the correct one is http://mirror.aarnet.edu.au/pub/CRAN/. Regards, Andrew C. Ward CAPE Centre, The University of Queensland Brisbane Qld 4072 Australia andreww at cheque.uq.edu.au

Dear David, If you have sufficient RAM you can increase the amount R uses by adding the switch "--max-mem-size=1024M" to the command you use to start R (change 1024 to whatever amount of RAM you have). If you don't have lots of RAM, you may have to draw just a subset of data. If you're developing or testing a model this will be okay, as long as you choose enough samples.

Hello R users, Thank you very much to those who responded so quickly and helpfully to my question about update.packages() and Windows 2000. I am grateful to Thomas Lumley, Jason Turner, Prof Brian D Ripley, Uwe Ligges and Andy Liaw (at least). The solution was simply to use internet2.dll rather than internet.dll. I'm sorry I didn't spot this in the documentation. Thank you again. As

Dear R-Help, I'm hoping to find a Multinomial Nested Logit package in R. It would be great to find something analogous to "PROC MDC" in SAS: > The MDC (Multinomial Discrete Choice) procedure analyzes models > where the > choice set consists of multiple alternatives. This procedure > supports conditional logit, > mixed logit, heteroscedastic extreme value,

Dear R users, I am trying to convert a dataset from one format to several rectangular datasets. A consultant helped design the data entry program for our survey using Delphi/Pascal and for each household the information is stored in a file called "EA-HM-HH.TXT" where EA is the enumeration area number, HM is the homestead number and HH is the household number. Within this file the

Is there any library for fitting multinomial probit using either likelihood or the "method of simulated moments" or both. I presume it would be possible to write a family function in VGAM for the multinomial probit, but was hoping that someone has done it already. Thanking you as always. Vumani Dlamini CSO-Swaziland

Hello R users, I'm using xyplot (lattice) to plot some data. Unfortunately, the text that goes into each strip is often too big to fit. For instance, one of the factors may be labelled "Satisfaction with waiting time" which is too wide to fit. If I use split.string() to wrap it onto two lines it is then too high to fit in the strip. I'd prefer to apply split.string on the

Dear R users, Sorry if this has been discussed before. I was recently looking at the MySQL site and saw that they offer various levels of paid support. The support options range from basic email support (USD200) through to extensive telephone support from the development team (USD10000). These are annual prices. This is something I would be happy to see available with R (maybe it is?). My