Re: [R] Fitted Value Pareto Distribution

2007-06-13 Thread Markus Jäntti
livia wrote: I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the fitted turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha),

Re: [R] generalized moment method

2007-06-11 Thread Markus Jäntti
Sebastian Kruk wrote: Dear everyone: I have to finish my thesis to graduate as Bs. in Economics. I choose to estimate a New Keynesian Phillips Curve (NKPC) for Uruguay using Generalized Moment Method (GMM). I do not know programming or R but I would like to use it. Should I use gee,

Re: [R] help with gls

2007-01-06 Thread Markus Jäntti
On Sat, 2007-01-06 at 01:21 -0500, Zhenqiang Lu wrote: Hello R-users, I am using gls function in R to fit a model with certain correlation structure. The medol as: fit.a-gls(y~1,data=test.data,correlation=corAR1(form=~1|aa),method=ML) mu-summary(fit.a)$coefficient With the toy data I

Re: [R] using weighted.mean with tapply()

2005-08-03 Thread Markus Jäntti
roger bos wrote: I am trying to calculate the weighted mean for a of 10 deciles and I get an error: decile - tapply(X=mat$trt1m, INDEX=mat$Rank, FUN=weighted.mean, w=mat$mcap) Error in FUN(X[[1]], ...) : 'x' and 'w' must have the same length All three of my inputs have the same length,

Re: [R] poly() in lm() leads to wrong coefficients (but correct residuals)

2005-06-29 Thread Markus Jäntti
On Wed, 2005-06-29 at 18:19 +0200, Andreas Neumann wrote: Dear all, I am using poly() in lm() in the following form. 1 DelsDPWOS.lm3 - lm(DelsPDWOS[,1] ~ poly(DelsPDWOS[,4],3)) 2 DelsDPWOS.I.lm3 - lm(DelsPDWOS[,1] ~ poly(I(DelsPDWOS[,4]),3)) 3 DelsDPWOS.2.lm3 -

Re: [R] Problem loading library Design

2005-03-24 Thread Markus Jäntti
Emmanuelle Comets wrote: I have used library Design (Frank Harrell) in the past but when I tried this week, I get : library(Design) Error in eval(expr, envir, enclos) : Object formula.default not found I updated to the latest version of the library (2.0) with the same result. My version of R

Re: [R] Quantiles of data in a contingency table

2005-03-17 Thread Markus Jäntti
Matt Mohebbi wrote: Hello, I have data of the following form: data - data.frame(type=c(c,d,e), size=c(10,20,30), count=c(20,10,5)) data type size count 1c 1020 2d 2010 3e 30 5 I would like to compute the quantiles of size given the counts. For instance, in this

Re: [R] Std Err on Concentration measures

2005-02-07 Thread Markus Jäntti
Angelo Secchi wrote: Hi, I'm using the ineq package to calculate some concentration measures (Gini, Herfindal, ...) and I was wondering if there's around also a function to calculate standard error on these measures. If not, is anybody aware of where I can find a reference on this point?

Re: [R] weighted mean

2003-11-26 Thread Markus Jäntti
/mailman/listinfo/r-help -- Markus Jäntti [EMAIL PROTECTED] Abo Akademi University __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] rbind.data.frame: character comverted to factor

2003-02-12 Thread Markus Jäntti
- rbind(d1, d2) str(d1) # - id is character str(d2) # - id is character str(d3) # - id is factor Regards, Markus -- Markus Jäntti [EMAIL PROTECTED] Statistics Finland __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help