search for: dlamini

...llowing 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 like c(3,0,0) using "lapply" but I have no splitting cases of the form c(1,1,0) and thats where I need your help. Thanking you as always. Vumani Dlamini

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 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

...r2[indices$var2]<-mean(xvar) temp$xvar3[indices$var3]<-mean(xvar) totalmodel<-clogit(status~intercept1+intercept2+intercept3+pspline(xvar1)+pspline(xvar2)+pspline(xvar3)+strata(stratum),method="exact",temp) coef(totalmodel) >From: ripley at stats.ox.ac.uk >To: Vumani Dlamini <dvumani at hotmail.com> >CC: R-help at stat.math.ethz.ch >Subject: Re: [R] R analogue >Date: Mon, 20 Jan 2003 13:17:02 +0000 (GMT) > >On Mon, 20 Jan 2003, Vumani Dlamini wrote: > > > Is there any R analogue for the S+ function "na.gam.replace". > >No,...

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

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"] <-

...t very promising: 1. It is my understanding that the multinom function in the nnet package does NOT do nested models. 2. The MNP (multinomial probit) package is close, but I want a logit, not probit, model. 3. http://www.r-project.org/nocvs/mail/r-help/2002/4394.html On Wed, 29 May 2002, Vumani Dlamini wrote: > > Has anyone implemented the conditional logit and/or the nested > logit model > > in R? > > Conditional logistic regression is clogit() in the survival package. But the survival package didn't seem to have what I needed. Am I overlooking anything? Are there...

Dear useRs, Am trying to use xtable on the following data.frame and I don't get what I expect: example.table <- data.frame(rbind( c("Gender"," "," "," "), cbind(rep(" ",2),c("Male","Female"),c(3.0,4.0),c(3/7,4/7)) )) colnames(example.table) <- c(" "," ","number of

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

...50 There are several similar files where the EA number, HM number or HH number changes, thus I would also like to know whether it is possible to create a script where for all household there is only one dataset for demographics, and one for agricultural inputs. Thanking you as always. Vumani Dlamini

...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 _________________________________________________________________ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subsc...

...as xyplot()? In that case, the legend (key in lattice) can be placed pretty much anywhere. Regards, Andrew C. Ward CAPE Centre Department of Chemical Engineering The University of Queensland Brisbane Qld 4072 Australia andreww at cheque.uq.edu.au On Friday, January 17, 2003 10:22 PM, Vumani Dlamini [SMTP:dvumani at hotmail.com] wrote: > 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 le...

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, 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 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

Dear Listers, I am trying to use "rbinom" in my C code, but i always get zeros as output no matter the probability. Am not sure what I am doing wrong because the function has worked before. Attached in an example. Noticed that "rbinom" expects 'n' to be REAL. Regards, Vumani R 2.3.1 (2006-06-01) Windows XP Gcc /* Called this file binom.c and then ran rcmd shlib on it

Dear listers, Am currently using MCMC approaches to estimate some parameters of my model. One parameter has to be updated using a tuned gamma distribution. So at each iteration I estimate the mean and variance of the density of the gamma approximation using "vmmin" (i also supply the gradient argument). For moderate replications the procedure works, but if I increase them R crashes.

Dear useRs, I'm trying to estimate the proportion of individuals with a without a certain recurring ailment at several times points. The data are of the survival type, with "start"-"stop" dates and whether the individual had the ailment in that interval. Some cases are observed until database closure and some died or are lost to followup. The interest is not on death. I