similar to: nested design; intercept

Dear R-help, I posted this on 4 Oct but got no response (I wasn't even told to go away and do some more background reading ;) ). I am reposting it in the, perhaps, vain hope that someone with knowledge of the subject will reply, if only to point me in a different direction to which I am now facing. Earlier Posting:--- I am sorry if this is more of a stats question than an R-question, but I

Dear List, I would greatly appreciate help on the following matter: The RuleFit program of Professor Friedman uses partial dependence plots to explore the effect of an explanatory variable on the response variable, after accounting for the average effects of the other variables. The plot method [plot(summary(rq(y ~ x1 + x2, t=seq(.1,.9,.05))))] of Professor Koenker's quantreg program

Dear List, I would appreciate help on the following matter: I am aware that higher dimensional contingency tables can be analysed using either log-linear models or as a poisson regression using a generalized linear model: log-linear: loglm(~Age+Site, data=xtabs(~Age+Site, data=SSites.Rev, drop.unused.levels=T)) GLM: glm.table <- as.data.frame(xtabs(~Age+Site, data=SSites.Rev,

Dear R-help, Has anyone implemented sparklines in the strips of a lattice plot? What I have in mind is, say, highlighting that part of a time series that one is examining in more detail in a set of lattice plots. Regads,. Mark Difford. PS: (Andreas Loffler has implemented a simple but functional version for TeX/LaTeX: http://www.tug.org/tex-archive/help/Catalogue/entries/sparklines.html)

Dear R-core/development-team, The problem noted in the subject-line has been a problem in the last three development versions of R for Windows that I have downloaded and tested, the most recent of them being a version I downloaded this morning. Update.packages() using the graphical option, i.e. called as update.packages(ask='graphics', checkBuilt=TRUE) does not work as it should, but

On May 24, 2011; 11:06pm Belle wrote: > Does anyone know how to run Multinomial logistical Model in R in order to > get predicted probability? Yes. I could stop there but you shouldn't. The author of the package provides plenty of examples (and two good vignettes) showing you how to do this. Suggest you do some work in that area. Look especially at how model formulas are

On Nov 12, 2011 at 8:29pm Alex wrote: > Has anyone used SIAR package add on? I posted a reply to an earlier question from you on this subject. See http://r.789695.n4.nabble.com/Errors-in-SIAR-td4029804.html. In it I note that there are problems with the function from siar (not SIAR) you are using, but that this may not be "your" problem, that the function calls for matrices (you

Dear R-users, I am a newbie to this site and a relative new-comer to S/R, so please tread lightly, for you tread... There have been several posting relating to problems with augPred() from the nlme library. Here is a "fix" for one of these problems which may lie at the root of others. In my case the problem with augPred() lay in gsummary(), which augPred() uses, causing it to fail.

Thanks Dimitris!!! That's much clearer now. Still have a lot of work to do this weekend to understand every bit but your code will prove very useful. Cheers, Aziz -----Original Message----- From: Dimitrios Rizopoulos [mailto:Dimitris.Rizopoulos at med.kuleuven.be] Sent: May 12, 2006 4:35 PM To: Chaouch, Aziz Subject: RE: [R] Maximum likelihood estimate of bivariate

HI Every one Could some one provide me definitions of following bivariate distributions gamma, exponencial, Weibull, half-normal , Rayleigh, Erlang,chi-square thanks A.S. Qureshi

Hi All, I have been trying to get to plot my PCA legend outside of the PCA plot, but success still alludes me. Can you guys please advise how I can achieve this. I used locater() to obtain coordinates for below the Comp.1 axis. Using these coordinates the legend disappears. Below is the code for the PCA and legend. Thanks in advance for the help. Regards Tinus r.cols <-

Hi, I'm dealing with wind data and I'd like to model their distribution in order to simulate data to fill-in missing values. Wind direction are typically following a vonmises distribution and wind speeds follow a weibull distribution. I'd like to build a joint distribution of directions and speeds as a VonMises-Weibull bivariate distribution. First is this a stupid question? I'm

Hi: I have a large dataset mydata, of 1000 rows and 1000 columns. The rows have gene names and columns have condition names (cond1, cond2, cond3, etc). mydata<- read.table(file="c:/file1.mtx", header=TRUE, sep="") I applied PCA as follows: data_after_pca<- prcomp(mydata, retx=TRUE, center=TRUE, scale.=TRUE); Now i get 1000 PCs and i choose first three PCs and make a

Dear List-Members, Is the application of multiple comparison procedures (using the multcomp package) to the output from a rank-based ANOVA straightforward, or do I need to take heed ? That is, is it as simple as: glht( aov(rank(NH4) ~ Site, data=mydat), linfct=mcp(Site="Tukey") ) Thanks in advance for your help. Regards, Mark Difford. -- View this message in context:

Hi, I have a right-censored (positive) random variable (e.g. failure times subject to right censoring) that is observed for N subjects: Y_i, I = 1, 2, ..., N. Note that Y_i = min(T_i, C_i), where T_i is the true failure time and C_i is the censored time. Let us assume that C_i is independent of T_i. Now, I would like to generate another random variable U_i, I = 1, 2, ..., N, which is

Greetings, I'm interested in generating data from various bivariate or mulitivariate distributions (e.g. gamma, t, etc), where I can specify the parameter values, including the correlations among the variables. I haven't been able to dig anything up on the faq, but I probably missed something. A nudge in the right direction would be appreciated. David --

Hello Anyone that now of any function in R that can calculate copulas? Or if anyone have any code avaible I would be more than interested. Thank you in advance /Thomas ______________________________________________ R-help at stat.math.ethz.ch mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Hi All, I'm using the following code to produce some stacked bar graphs. *setwd("C:\\Users\\Tinus\\Documents\\NMMU\\R\\Seamounts")* *SChla <- read.csv("SM_Chla_data.csv")* * * *#Extract mean values from data file* * * *Coral <- SChla[185:223,] #Reduce SChla to Coral only* *coral <- with(Coral , aggregate(cbind(Pico, Nano, Micro), list(Depth),FUN=mean))*

Good day, I read some data from a PostgreSQL database by a following script: library(Rdbi) library(RdbiPgSQL) # conn becomes an object which contains the DB connection: conn <- dbConnect(PgSQL(), host="localhost", dbname="BVS", user="postgres", password = "*******") query_duj_kal <- dbSendQuery(conn, "select zdroj as well, cas as date, fe2,

Hi all, I'm trying to do model reduction for logistic regression. I have 13 predictor (4 continuous variables and 9 binary variables). Using subject matter knowledge, I selected 4 important variables. Regarding the rest 9 variables, I tried to perform data reduction by principal component analysis (PCA). However, 8 of 9 variables were binary and only one continuous. I transformed the data by