search for: vumani

...d so on ##### #fit model to data # } When I use one year there is no problem, but I would like some data sets to span over years and I am not sure how to do this without having to change the body of my fitting program. I tried searching the R-list but to no avail. Thanking you in advance. Vumani _________________________________________________________________

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

...; is doing the right thing (I may be wrong). The strange thing here is that even the intercepts are the same. Lastly, these coefficients are the same as using clogit before adding the constant, but the intercept are different. I have attached my code for all this, maybe I messed it somewhere. Vumani ########################## xvar<-rnorm(100,0,1) data.mult<-rbinom(100,size=3,prob=(exp(0.5+0.4*xvar)/(1+exp(0.5+0.4*xvar)))) library(Nnet) mult.fit<-multinom(data.mult~xvar) coef(mult.fit) library(survival) choice<-c(ifelse(data.mult==0,1,0),ifelse(data.mult==1,1,0),ifelse(data.mult==2...

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

...), : C/Fortran function name not in load table When I use "transcan" it is; >f <- transcan(~y + x1 + x2, n.impute=10, shrink=FALSE) Error in transcan(~y + x1 + x2, n.impute = 10, shrink = FALSE) : Must specify data= when using R I am not sure what I am missing. Vumani

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

...t; , ", beta==scaleU,") Mean = ",meanU," days"), list(shapeU = round(mle.r$estimate[1],3), scaleU = round(mle.r$estimate[2],3), meanU = round(mle.r$estimate[3),4)))) ), lty=c(1,1),col=c(1,2),lwd=c(1,1),merge = TRUE) (Windows XP + R2.10 user) Thanks, Vumani

...connections inlcude a circle with the information on the indices ("daughter index", "mother index"): i.e. something similar to graphs produced by graphviz. I am clueless. I have looked at the libraries SEM and SNA but it looks like they can't help my cause. Hope you help. Vumani

...using glm. I use poistab <- data.frame(table(poisdata)) colnames(poistab) <- c("width","freq"); poistab[,"width"] <- as.numeric(poistab[,"width"]) glm(freq ~ 1-width, family = poisson, data = poistab) but the results are always different. cheers, vumani

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 */#include <R.h>#include <Rmath.h>#include <math.h>#include <Rdefines.h>SEXP binomial(SEXP r, SEXP n, SEXP p){ int i; SEXP out; PROTECT(out = allocVector(REALSXP,INTEGER(r)[0]));...

...e was a memory leak in "vmmin". I am not an expert in "C" and am only a windows user. Have tried some of the debugging tools available to us windows users: gdb, mpatrol, duma without success. Am not sure whether this is a bug. Version 2.3.1 (2006-06-01) Windows XP Regards, Vumani Rgui.exe caused an Access Violation at location 100dec2f in module R.dll Reading from location 3ff7c27d. Registers: eax=000063ac ebx=01e773a0 ecx=3ff7c27a edx=00000002 esi=01e254b4 edi=00000002 eip=100dec2f esp=00e0ecb0 ebp=00e0ed78 iopl=0 nv up ei pl nz na pe nc cs=001b ss=0023 ds...

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

...d/or not 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? Ar...

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

...001 0050 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

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

..., such 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 hav...

...the month resulting in proportions above 1 sometimes, for instance at the beginning where no individuals have the ailment but some are censored before the midpoint time. I have reason to believe there is a way within the survival packages but haven't figured out how. Thanking you in advance, Vumani ps: sorry for cross posting. wasnt sure whether "r-sig-epi@stat.math.ethz.ch" or "r-help@R-project.org" was more suitable for this question _________________________________________________________________ [[alternative HTML version deleted]]