similar to: multi-comparison of means

Dear all, downloaded alpha13, did: ./autogen.sh This result in: [root at node1 source4]# ./autogen.sh Setting up for waf build done. Now run ./configure or ./configure.developer then make. ./configure.devloper ?ok Then: Make Error: data.mk:1881: *** Befehle beginnen vor dem ersten Ziel. Schluss. Make does nothing?! Any idea1? Daniel ----------------------------------------------- EDV Daniel

Hi, A new try: What I did. Took the old source directory I made for samba4 12: Cd /samba4 Rsync the new samba4 13 in it by: rsync -avz samba.org::ftp/unpacked/samba_4_0_test/ . cd samba4/source made: make clean then ./autogen.sh then ./configure.developer then make idl_full The same error nothing: data.mk:1389: *** Befehle beginnen vor dem ersten Ziel. Schluss. The file data.mk is the newest

Hello, I've got a matrix consisting of one column with n rows. Each field in the matrix is filled with a character vector. I would like to convert this matrix into a character vector containing the B03_MAH-type entries from the beginning of each row and a data.frame that contains the numeric data, but I am stuck. I have tried to use textConnection with sep = " " but since

Dear all, It looks like the parameters var-cov matrix returned by gee() is not adjusted for the scale parameter: > fm1 <- gee(nbtrp ~ strate * saison + offset(log(surf)), family = poisson, data = Eff2001, + id = loc, tol = 1e-10, corstr = "exchangeable") [1] "Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27" [1] "running glm to get initial

This command cdmoutcome<- glm(log(value)~factor(year) > +log(gdppcpppconst)+log(gdppcpppconstAII) > +log(co2eemisspc)+log(co2eemisspcAII) > +log(dist) > +fdiboth > +odapartnertohost > +corrupt > +log(infraindex) > +litrate > +africa >

Dear WizaRds, A run of nls produces the following concise summary: > summary(cs.wt) Formula: 0 ~ wt.MM(conc, time, A1, a1, A2, a2) Parameters: Estimate Std. Error t value Pr(>|t|) A1 4.814e+02 2.240e+01 21.495 0.0296 * a1 7.401e-01 7.435e-02 9.956 0.0637 . A2 1.613e+02 1.738e+01 9.280 0.0683 . a2 1.770e-02 7.324e-03 2.417 0.2497

Dear Dr.Goslee and anyone may intrested in matrix manipulate, I am using your ecodist to do mantel and partial mantel test, I have locality data and shape variation data, and the two distance matrixs are given as belowings. When I run the analysis, it is always report that the matrix is not square, but I didn't know what's wrong with my data. Would you please help me on this. I am quite

I have a data frame with 1 factor, one date, and 37 numeric values: str(waterchem) 'data.frame': 3525 obs. of 39 variables: site : Factor w/ 64 levels "D-1","D-2","D-3",..: 1 1 1 1 1 ... $ sampdate : Date, format: "2007-12-12" "2008-03-15" ... $ CO3 : num 1 1 6.7 1 1 1 1 1 1 1 ... $ HCO3 : num 231 228 118 246

I have a fairly large set of data with the following attributes: >str(raw.data) `data.frame': 1429 obs. of 16 variables: $ TStamp :`POSIXlt', format: chr "2001-11-25 02:00:00" "2001-11-25 01:55:00" "2001-11-25 01:50:00" "2001-11-25 01:45:00" ... $ iPDT.AHU14.14: num 0.0122 0.0125 0.0120 0.0120 0.0122 ... $ iPDT.AHU14.15: num 0.0121

Hello!!, for linear models fit I use Gretl, but now I'm starting to use R, I would like to know if is there some function to obtain a extended summary like in Gretl. I will write a example in Gretl Modelo 1: MCO, usando las observaciones 1968-1982 (T = 15) Variable dependiente: Invest Coeficient St error t-ratio p-value const 377,631 35,0955 10,7601 <0,00001 *** GNP

Hola!! A ver si alguien puede ayudarme!! Para ajuste de modelos lineales normalmente uso Gretl. Ahora estoy empezando a hacerlo en R. Me gustaría saber si existe alguna función que haga un summary extendido como el de Gretl. Os pongo un ejemplo del summary de Gretl. Modelo 1: MCO, usando las observaciones 1968-1982 (T = 15) Variable dependiente: Invest              Coeficiente   Desv. Típica

Hello All, I have searched many help forums, message boards, etc. and I just can't apply the comments to what I need my program to do. I am running R 2.5.1 on an XP system, and my desire is to produce replicate datasets for a simulation study I am running. Essentially, I have sets of parameters (a's, b's, and c's) that define a function which produces a decimal value. This

It seems that the following should work: > x<-read.table('1A.tab',header=T,skip=1); > mode(x[1]) [1] "list" > x[1] p01 R1-00 0.0295 R1-01 0.0283 R1-02 0.0145 R1-03 0.0235 R1-04 0.0339 R1-05 0.0239 R1-06 0.0414 R1-07 0.0259 > hist(x[1]) Error: hist: x must be numeric > as.numeric(x[1]) Error: (list) object cannot be coerced to vector type 14 >

I am playing with the a 1-way anova with and without the "-1" option. I have a simple cooked up example below but it behaves the same on a more complex real example. From what I can tell: 1) the estimated means of the different levels are correctly estimated either way (although reported as means with the -1 and as contrasts without the -1 as expected) 2) the residuals are

dear listers, I use gam (from mgcv) for evaluation of shape and strength of relationships between a response variable and several predictors. How can I interpret the 'F' values viven in the GAM summary? Is it appropriate to treat them in a similar manner as the T-statistics in a linear model, i.e. larger values mean that this variable has a stronger impact than a variable with smaller F?

I'm using the fixed argument in arima. Shouldn't ar4, ar5, and ar6 display as zero in the output? Call: arima(x = window(log(hhprice), start = c(1990, 1), end = c(2003, 3)), order = c(7, 1, 0), xreg = window(ts.union(exa1 = lag(exa, -1), exa12 = lag(exa, -12), exb1 = lag(exb, -1), exc1 = lag(exc, -1), exc12 = lag(exc, -12)), start = c(1990, 1), end = c(2003, 3)),

Hola Ruben, Sí precisamente es lo que comentas, en matemáticas no se suele llamar bucketización (este término se emplea más en informática) sino datos agrupados. Pero la idea es la que tu mismo dices. Respecto a las gráficas que has puesto, me han aclarado mucho sobre el tema, gracias. Si realizo lo mismo, por ejemplo con nbucket=1000 sigo obteniendo un p-valor de 1. Es decir, que casi le

On Fri, Sep 14, 2007 at 11:42:18PM -0700, Tanya Lattner wrote: > The 2.1 pre-release (version 1) is available for testing: > http://llvm.org/prereleases/2.1/version1/ > > [...] > > 2) Download llvm-2.1, llvm-test-2.1, and the llvm-gcc4.0 source. > Compile everything. Run "make check" and the full llvm-test suite > (make TEST=nightly report). > > Send

I am struggling to generate p values for comparisons of levels (post-hoc tests) in a glm with a negative binomial distribution I am trying to compare cell counts on different days as grown on different media (e.g. types of cryogel) so I have 2 explanatory variables (Day and Cryogel), which are both factors, and an over-dispersed count variable (number of cells) as the response. I know that both

Hi all, I'm handling massive data.frames and matrices in R (30000 x 400). In the 1st column, say, I have 0s and 1s indicating rows that matter; other columns have probability values. One simple task I would like to do would be to get the column mean values for signaled rows (the ones with 1) As a very fresh "programmer" I have build a simple function in R which should not be very