similar to: Change default parameters of panel.smooth

How do you install a package from CRAN ? I want to install pls.pcr, so I have downloaded pls.pcr_0.1.1.tar.gz but when I try to install it using the install package(s) from local zip file(s) option it says : > install.packages(choose.files('',filters=Filters[c('zip','All'),]), .libPaths()[1], CRAN = NULL) Error in file(file, "r") : unable to open connection

Does anybody know how to make pairs graphics with dots of different colours depending on the value of a categorical variable ? Thanks, Arnaud ************************* Arnaud DOWKIW Department of Primary Industries J. Bjelke-Petersen Research Station KINGAROY, QLD 4610 Australia T : + 61 7 41 600 700 T : + 61 7 41 600 728 (direct) F : + 61 7 41 600 760 **************************

Hi all, I have received the following error from summary.manova: Error in summary.manova(manova.test, test = "Pillai") : residuals have rank 36 < 64 The data is simulated data for 64 variables. The design is a 2*2 factorial with 10 replicates per treatment. Looking at the code for summary.manova, the error involves a problem with qr(). Does anyone have a suggestion as to how to

Could anyone please explain how to set up contrasts between means in R. I want to know if "before I conduct an experiment and believe the mean for 1 and 2 will be different from means 3 and 4, Is this true?" That is what I have to prove or disprove, I thought that contrasts would be the way to go. Thanks for the help. Igor [[alternative HTML version deleted]]

Hi all, I am trying to get a handle on gls (package nlme). I have a toy problem: 3 fixed factors (A, B, C), two levels each, 5 replicates per treatment. The response variable is continuous, normal. I have a correlation matrix of the form: > mat A B C A 1.00 0.75 0 B 0.75 1.00 0 C 0.00 0.00 1 which is common to all observations. How do I construct the call to gls? I think I need to

Hi, everyone, I am a new user of R. Recently, I tried the following commands, but couldn't make them work. If any one of you has some ideas about it, please help me. The commands are >std<-1000 >mean<-8000 >prior<-function(n){1/(sqrt(2*pi)*std)*exp(-1.0*(n-mean)^2/(2*std^2))} >plot(prior,1,15000) >post<-function(n){ + if(n < 9543) + 0 + else + prior(n)/n +

Hello everyone, Does anyone know how to implement the argument "unit" in R's plclust function ? I used to use Splus where this argument exists but it has not been implemented in R's plclust. The reason why I switched from Splus to R is that Ward's method is not implemented for S's hclust whereas it is implemented for R's hclust. What I would need is S's plclust

Dear all, Does anyone know how to run Principal Coordinates Analysis (PCoA) from a squared euclidean dissimilarity matrix with R ? Thanks, ************************* Arnaud DOWKIW Department of Primary Industries J. Bjelke-Petersen Research Station KINGAROY, QLD 4610 Australia T : + 61 7 41 600 700 T : + 61 7 41 600 728 (direct) F : + 61 7 41 600 760 **************************

Dear R-helpers, I have performed a PLS regression with the mvr function from the pls.pcr package an I have 2 questions : 1- do you know if mvr automatically centers the data ? It seems to me that it does so... 2- why in the situation below does the output say that the optimal number of latent variables is 4 ? In my humble opinion, it is 2 because the RMS increases and the R2 decreases when 3 LVs

I am looping over many data files and reading in the data with F <- read.table(filename) to read in a 22000 by 15 matrix. Works fine on the first matrix F, but I get the following error when the second file is read into F: Error in "row.names<-.data.frame"(*tmp*, value = row.names) : duplicate row.names are not allowed I have tried picking a column of the matrix and

Hello, My name is Rodrigo, I am using R program and I have a trouble. I am trying to do a dendrogram with genetics information. Let me explain... The Similarity Matrix was already did, and with this matrix I want to construct a dendrogram. So, the distance is done. I need to transform this matrix (that I have) in a dendrogram, I woud be very grateful if someone could help me. PS: I am sending

Forgive me please if this has been answered previously, but the r-help archive doesn't have a "search" capability. It's a pity, as I'm sure that there's a wealth of informative help in there. I am unable to get the "cloud" function to produce a fundamental 3D scatterplot. I have a 3-column numeric matrix: >is.matrix(sanity.MIF); is.numeric(sanity.MIF);

Dera R helpers, I have a (x,y,z) data file where x and y are spatial coordinates and z a variable. I have some missing values in the z column and I would like to replace them with an optimized estimation from the neighbour cells. I could not find any function in R to do that. Is anybody aware of such function ? Maybe someone knows how to implement the Papadakis method in R... Thanks for your

dear group members, I am looking for a function that assess the stability of cluster. The result of hclust function is an hclust object which can be plot as a dendrogram. However to have confidence in the tree topology usualy bootstap is applied. I understand that I can apply bootstarp on the original data and then run hclust(dist() ) as much as I resampled but how to comapre the topologies the I

Hi, I'm probably being dense, but could somebody tell me how I can get tick marks on the axes of my wireframe plot to show up? This gets me a beautiful looking plot, but no ticks: wireframe(z ~ x * y, data, drape=T, colorkey=F) Thanks in advance, -James [[alternative HTML version deleted]]

I would appreciate some advice on the following task. I have some data that currently looks like this: t1 <- data.frame(id=c(1,1,2,2), aspect=c("A","B","A","B"), score=c(10,9,11,12)) I'd like it to look like this: id A B 1 10 9 2 11 12 reshape() looks like a good candidate for this job but I'm not really sure about the roles of

Hello I was wondering if it was possible to enter R commands from an external text file. If it is possible, it will be easy for repetetive tasks. Does anyone have an idea ? thanks in advance Ahmet Temiz TURKEY ______________________________________ ______________________________________ The views and opinions expressed in this e-mail message are ...{{dropped}}

Hi, Can someone set me straight as to how to write formulas in R to indicate: complete independence [A][B][C] joint independence [AB][C] conditional independence [AC][BC] nway interaction [AB][AC][BC] ? For example, if I have 4 factors: hair colour, eye colour, age, sex does > mosaicplot( frequency ~ hair + eye + age + sex) mean that the model fitted is of complete independence of

Dear all, So far i could do (in an informal way) to draw a Standardized Resisual plot in the following way- --------------------- >x <- c(104.1, 106.6, 105.5, 107.5, 109.6, 113.3, 115.5, 117.7, 119.9, 122.1, 124.3, 126.5, 128.2) >y <- c(53732, 52912, 57005, 61354, 67682, 71602, 71961, 75309, 82931, 93310, 102161, 103068, 108927) >

Hi R helpers, I was wondering if anybody knows if is possible to generate bootstrap values for a cluster analysis in R. What I am trying to do is obtain some confidence on the clusters formed by resampling the data set. A similar type of analysis is used in molecular taxonomy and the confidence values of each cluster are placed in the nodes of the dendogram. Any ideas on how to do this