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=
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
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
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
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
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é
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
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
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
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
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
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
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
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
) 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
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
?
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
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
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
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
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
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
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
)]-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
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
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
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
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
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
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
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
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
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
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
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
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
) + :
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
51 matches
Mail list logo