[R] Hmisc::latex and the use of \phantom

Dear all, Is it possible to avoid the use of \phantom with latex function ? when I run latex(latex=FALSE,…) I get an error message. When I use Format(big.mark= ‘’ ‘’) the result appear correct in R console but not in Latex code. Is it a way to combine numprint with latex ? col.just=

[R] two phases sampling

is y=y1+y2. So the variance of y is v(y)=v(y1)+v(y2)+cov(y1,y2). Sampling theory indicate how to compute v(y1) and v(y2) but how can I compute cov(y1,y2) ? Can the survey package help me for that ?   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted

[R] Summary.formula question

to add a marginal column ?  Sincerly   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] Re : P values for vglm(zibinomial) function in VGAM

Estimation is realized by MLE, estimators are asymptotically normal Try this reg-vglm(...) p.value-1-pnorm(abs(coef(reg)/sqrt(diag(vcov(reg)   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : suuz suuz_b...@hotmail.com À : r-help@r-project.org Envoyé le

[R] Obtain beta regression estimation for betareg and VGAM package

and all coefficients are doubled in the second estimation. Is it possible to obtain the same estimation with the two packages ?   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] Re : Obtain beta regression estimation for betareg and VGAM package

I have found the solution betareg(formula, data) is equivalent to vglm(formula,betaff(zero=2),data) Sorry for the previous post.   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 - Mail original - De : justin bem justin_...@yahoo.fr À : R Maillist r-h...@stat.math.ethz.ch Cc : Envoyé

[R] vglm regression with weibull distribution

by 1.819e-12 9: In eval(expr, envir, enclos) :   MLE regularity conditions are violated(shape = 2) at the final iteration MLE regularity condition are still violated but a solution is find ! Can some one help me ? Justin BEM BP 1917 Yaoundé Tél (237) 76043774

[R] Install.package error

When try to use install.package I have this error message : Error in m[, 1L] : incorrect dimensions number What the matter with my installation ? I run R2.13 on Window 7 32bits.   Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted

[R] SAS Mixed

it's stupid to loose parcimony by enlarging explanatory variables space) ? Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Re : descriptive statistics

A nice way to obtain summary for data is to use summary.formula in Hmisc package. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : Jim Lemon j...@bitwrit.com.au À : effeesse scarpin...@gmail.com Cc : r-help@r-project.org Envoyé le : Lun 13 décembre 2010

[R] Re : Re : descriptive statistics

variable together in all table my.summary-function(x) c(mean(x[,1]),mean(x[,2]),mean(x[,3]),mean(x[,4]),mean(x[,5])) summary(cbind(v1,v2,v3,v4,v5)~v,data=df) Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : effeesse scarpin...@gmail.com À : r-help@r-project.org

[R] Fitting GLM with BFGS algorithm

Dear all, By default the glm function in the stats package use IWLS. How can I fit a glm model using BFGS algorithm ? Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing

[R] Error with read.ssd

Dear all, I'm trying to read and sas file hen using read.ssd, I got this error message: symbolic links are not supported I have debug read.ssd function and the function file.symlink is used. I run windows XP on virtual box, hosted by Vista Home  [[alternative HTML version

[R] Re : SAS and R on multiple operating systems

There is a possibility of integration of SAS with and other statistic package like Stata, etc. with Latex by using statweave software. http://www.cs.uiowa.edu/~rlenth/StatWeave/ or http://stat.uiowa.edu/~rlenth/StatWeave/OLD/SRC-talk.pdf Sincerly. Justin BEM BP 1917 Yaoundé Cameroon Tél

[R] Re : calling combinations of variable names

) results[i]-testfun(row) } You just have to replace testfun by your own function in this case ICC. Sincerly. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 De : Helter Two helter...@care2.com À : r-help@r-project.org Envoyé le : Lundi, 7 Septembre 2009, 18h17mn 22s

[R] Heckman probit ?

Is there a function to fit heckman probit model in R ? Sincerly.. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] IIA test

? Sincerly. Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] Re : Support Vector Machines

there is a SVM function in svmpath by Trevor Hastie. Before use it take time to read this http://www-stat.stanford.edu/~hastie/Papers/svmpath.pdf If you install then svmpath library ! run the attach source file. Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246

[R] Re : Re : Support Vector Machines

Of course SVM is for supervise learning method (classification or regression). You cannot use a boat to fly man !   Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 Envoyé le : Mardi, 5 Mai 2009, 16h07mn 10s Objet : Re : [R] Support Vector

[R] Sub-saharan Africa R user forum

A new forum on R for french speaker in sub-saharan africa  http://aurass.forum-free.org Justin BEM BP 1917 Yaoundé Tél (237) 76043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] Re : PCA and automatic determination of the number of components

See ade4 or mva package.  Justin BEM BP 1917 Yaoundé De : nikolay12 nikola...@gmail.com À : r-help@r-project.org Envoyé le : Lundi, 20 Avril 2009, 4h37mn 41s Objet : [R] PCA and automatic determination of the number of components Hi all, I have relatively

[R] Re : function question

Yes  because in OOP the constructor have the same name with class  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Bert Gunter gunter.ber...@gene.com À : Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no; r-help@r-project.org Envoyé le

[R] Re : function question

Functions and vectors are both objects in R, but a functions are not vectors try this foo-function() print(I'm the foo function) is.vector(foo) as.vector(foo)  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Edna Bell edna.bel...@gmail.com

[R] Re : Randomly splitting a data frame in half

)]-1 data$selected-selected rm(selected) or data$selected-rbinom(39622,1,.5) select case have the value 1, non-selected have value 0. In the second case, you will not get exactly .5 ! Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De

[R] Re : BOOTSTRAP_CROSS VALIDATION

See the boot package  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : per243 jose.perezsua...@csiro.au À : r-help@r-project.org Envoyé le : Lundi, 16 Mars 2009, 0h12mn 03s Objet : [R] BOOTSTRAP_CROSS VALIDATION I need a script that works

[R] nested logit estimation

Dear all, Is there a function to fit nested logit available somewhere ? Sincerly  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Re : Converting a dataframe to a matrix

Something like this will work m-matrix(df1$likes, nr=3,nc=3,byrow=T) colnames(m)-unique(df1$color) rowlnames(m)-unique(df1$name)  Sincerly. Justin BEM BP 1917 Yaoundé Cameroun Tél (237) 76043774 De : Jennifer Brea b...@fas.harvard.edu À : r-help@r

[R] Re : t-test, survey data

For survey data, the survey package provides methods for variance estimation.  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Marie Vandresse v...@plan.be À : r-help@r-project.org Envoyé le : Mercredi, 11 Mars 2009, 9h25mn 11s Objet : [R] t

[R] Re : (no subject)

see ?unique Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : arnaud_mosn...@uqar.qc.ca arnaud_mosn...@uqar.qc.ca À : r-help@r-project.org Envoyé le : Mardi, 10 Mars 2009, 17h15mn 57s Objet : [R] (no subject) Dear R users, I have a table

[R] Re : How to optimize a matrix

See ?optim See ?mle (stats4 package) or help.search(MLE) Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : b88207...@ntu.edu.tw b88207...@ntu.edu.tw À : r-help@r-project.org Envoyé le : Lundi, 9 Mars 2009, 6h55mn 20s Objet : [R] How

[R] Re : PCA and categorical data

See homals package in R. But also look documents for ade4 package.    Justin BEM BP 1917 Yaoundé Tél (237) 76043774   De : Galanidis Alexandros a...@env.aegean.gr À : r-help@r-project.org r-help@r-project.org Envoyé le : Vendredi, 6 Mars 2009, 10h09mn 18s Objet

[R] Re : SVM regression code

there is svmpath package by Trevor Hastie  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Alex Roy alexroy2...@gmail.com À : r-help@r-project.org Envoyé le : Jeudi, 19 Février 2009, 9h19mn 18s Objet : [R] SVM regression code Dear R

[R] Re : standard error of logit parameters

Run outfit-nlm(..., hessian=T) and then standards error are se-diag(solve(outfit$hessian))   Justin BEM BP 1917 Yaoundé Tél (237) 76043774   De : Bomee Park bom...@stanford.edu À : r-help@r-project.org Envoyé le : Jeudi, 29 Janvier 2009, 4h01mn 56s

[R] Re : Need help on running Heckman Correction Estimation using R

See the micEcon package. there is and heckit function  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Kishore gladikish...@gmail.com À : r-help@r-project.org; r-h...@stat.math.ethz.ch Envoyé le : Mardi, 27 Janvier 2009, 11h54mn 00s

[R] Specify nonlinear equation with nlme

Dear all, I want to estimate a nonlinear mixed effect model using nlme package or lme4. How can I specify this y=exp(B'X)/(1+exp(B'X)). the same model can be obtained with logit(y)=B'X but I have many value of 1 in y logit(1)=Inf. Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237

[R] specify formula for nonliear mxed model equation

How can I specify formula using nlme to estimate y=exp(X'B)/(1+exp(X'B)) ? Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] problem with Surv

) + : Response must be a survival object Can you help me ? Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] Re : A package to set up a questionnaire enter data

I don't believe there exist one ! but for that purpose I use CSPRO who is free !  http://www.census.gov/ipc/www/cspro/  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 À : r-help@r-project.org Envoyé le : Jeudi, 11 Décembre 2008, 1h35mn

[R] Re : growing a list sequentially -- memory management

that 100 loops for ( i in 1: length(sh)){     index-st(i):st(i+1)     myres(tapply(index, factor(index),gen     sqlwrite(myres,connexion,append=(i1)) } It will be more faster .  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Patrick

[R] A R forum for sub-saharan African R user

in the describe in fortune(83) Justin BEM BP 1917 Yaoundé Tél (237) 0023776043774 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] Re : Regarding lme functions

in nlme package you the intervals function. myfit-lme(...) intervals(myfit) Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : ratna ghosal [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Lundi, 10 Novembre 2008, 17h51mn 51s Objet

[R] sample selection mixed-effects model with nlme

Dear all, I want to estimate this equation with repeated data. log(y)=b0+b1*log(x1)+b2*log(x1)+B'Z+u if y=0 then log(y) is undefined.  It seem to be a tobit model with panel data, can I do this in R with nlme, lme or plm ? If not is it valid to use a two step procedure ? Justin BEM BP

[R] Re : panel data maximum likelihood function

see plm, or lmer packages Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Migle Purzelyte [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Lundi, 3 Novembre 2008, 0h26mn 21s Objet : [R] panel data maximum likelihood function Hi

[R] Re : how to restrict a parameter in optim()

Where is the problem ? then you have to write the objective function with the restriced value before call optim() Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Luis SAGAON TEYSSIER [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le

[R] Re : (a) Credit Scoring models and (b) aceesing earlier emails

Go to this link http://defaultrisk.com/ you will find a lot about credit scoring Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : saggak [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Mardi, 14 Octobre 2008, 9h11mn 53s Objet : [R

[R] Multilevel SEM

Dear all, Is it possible to estimate a SEM with panel data in R ? Sincerly  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] Tool for simple two class Fisher linear discriminant analysis?

Use lda in MASS package Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : Yip [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Jeudi, 18 Septembre 2008, 8h04mn 23s Objet : [R] Tool for simple two class Fisher linear discriminant analysis

[R] Re : How to import package into R script

Read documentation carefully ! see ?library (e.g library(tseries)) Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : Edward Wijaya [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Vendredi, 23 Mai 2008, 4h23mn 06s Objet : [R] How

[R] Re : Unable to save plot as .pdf , ps , eps

I simply use pdf(test.pdf) plot(rnorm(10)) dev.off() postscript(test.ps) plot(rnorm(10)) dev.off() Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : Luis Ridao Cruz [EMAIL

[R] Re : writing a categorical var. with condition

see also ?cut Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : [EMAIL PROTECTED] [EMAIL PROTECTED] À : sigalit mangut-leiba [EMAIL PROTECTED] Cc : r-help [EMAIL PROTECTED] Envoyé le : Vendredi, 2 Novembre 2007, 22h33mn 52s Objet : Re: [R

[R] Re : how to create this graph?

Try this plot.new() plot.window(xlim=c(0,10), ylim=c(0,15)) rect(1,0,2,10) rect(2,0,4,4.1,col=grey) rect(2,4.1,4,10,col=grey) segments(0, 0, 5, 0) segments(0, 10, 5, 10) text(3,3.5,41) text(3,6, 51) text(4.5,3..5,G/V) text(4.5,6, U) Justin BEM BP