[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= ‘n{#}{

[R] two phases sampling

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

[R] Summary.formula question

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

[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 À : r-help@r-project.org Envoyé le : Mardi 2

[R] vglm regression with weibull distribution

epsilon) :   1525 elements replaced by 1.819e-12 8: In checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) :   1525 elements replaced by 1.819e-12 9: In eval(expr, envir, enclos) :   MLE regularity conditions are violated(shape <= 2) at the final iteration MLE regularity co

[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 À : R Maillist Cc : Envoyé le : Mercredi 17 Août 2011 20h30 Objet

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

hapes 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/listi

[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

No" for me 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/mai

[R] Re : Re : descriptive statistics

all 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 À : r-help@r-project.org Envoyé le :

[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 À : effeesse Cc : r-help@r-project.org Envoyé le : Lun 13 décembre 2010, 11h 23min 15s Objet : Re: [R

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

[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

;", n) for (i in 1:n) { row<-unique(unlist(eval[i,])) if (length(row)>=3) 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 : H

[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

this case ? 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.o

[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] 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] 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 À : 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 small dataset on which

[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 À : Wacek Kusnierczyk ; r-help@r-project.org Envoyé le : Jeudi, 19 Mars 2009, 21h16mn 13s Objet : Re: [R

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

:39622,39622/2)]<-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 __

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

[R] Re : BOOTSTRAP_CROSS VALIDATION

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

[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 : 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 À : r-help@r-project.org Envoyé le : Mercredi, 11 Mars 2009, 9h25mn 11s Objet : [R] t-test, survey

[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 À : r-help@r-project.org

[R] Re : (no subject)

see ?unique Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : "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 with the foll

[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" À : r-help@r-project.org Envoyé le : Lundi, 9 Mars 2009, 6h55mn 20s Objet : [R] How to

[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 À : "r-help@r-project.org" Envoyé le : Vendredi, 6 Mars 2009, 10h09mn 18s Objet : [R] PCA and categorica

[R] Re : Have a function like the "_n_" in R ? (Automatic count function )

[1]),nr=1))     if (i==4) break } colnames(dat)<-c("counter", "x","seed") dat  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 De : Nash À : r-help Envoyé le : Mercredi, 25 Février 2009, 13h25mn 18s Objet : [R] Have a

[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 À : r-help@r-project.org Envoyé le : Jeudi, 19 Février 2009, 9h19mn 18s Objet : [R] SVM regression code Dear R user,           Â

[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 À : r-help@r-project.org Envoyé le : Jeudi, 29 Janvier 2009, 4h01mn 56s Objet

[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 À : r-help@r-project.org; r-h...@stat.math.ethz.ch Envoyé le : Mardi, 27 Janvier 2009, 11h54mn 00s Objet : [R] Need help on

[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

[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://st

[R] problem with Surv

d <= 1) ~ SOLVA + log(AF088) + : 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 http

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

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

ve around 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=(i>1)) } It will be more faster .  Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246

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

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 gu

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

[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 ? Ju

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

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

[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

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

[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] About Bootstrap

more accurate. Is it an explanation to this ? Is bootstrap appropriate for all situations or there are drawback? Sincerly. Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246

[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

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

[R] Re : how to create this graph?

text(4.5,3..5,"G/V") text(4.5,6, "U") Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : squall44 <[EMAIL PROTECTED]> À : r-help@r-project.org Envoyé le : Lundi, 1 Octobre 2007, 8h41mn 44s Objet : [R] how to creat

[R] Re : 财务心声。

So you expect answer to your post ! hum Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : 宝华工业 <[EMAIL PROTECTED]> À : r-help <[EMAIL PROTECTED]> Envoyé le : Jeudi, 13 Septembre 2007, 5h26mn 23s Objet : [R] 财务心