similar to: Variance-covariance by factor

Dear all; Searching the R site for answers to my problem, but found none. Here is the run and error: > bwplot(dev ~ Dbhcl | Period, data = DbhValid2, font = 2, + main=list(" "), axis.font =2, + ylab = list(label = "Residual (cm)", font = 2), + xlab = list(label = "Dbh class (cm)", font = 2), + par.strip.text=list(cex=0.8, font=2), + panel=function(x,

Dear All; I tried to use fitdistr() in the MASS library to fit a mixture distribution of the 3-parameter Weibull, but the optimization failed. Looking at the source code, it seems to indicate the error occurs at if (res$convergence > 0) stop("optimization failed"). The procedures I tested are as following: >w3den <- function(x, a,b,c)

Dear All; Searching on the achieve, many questions on optim() have been asked, but I haven't seen the following. The question began with my original inquiry on "Optimization failed in fitting mixture 3-parameter Weibul l distribution using fitdistr()" which I posted on Jul. 28, Prof. Ripley kindly advised me to look into options of optim() for the answer. Following his advice and

Hello, I installed RGui 2.5.1 and package DBI on Windows XP and tried to connect to Oracle database which is on a Linux server. When I tried to use dbDriver("Oracle"), I got an error as follows: > drv <- dbDriver("Oracle") Error in do.call(as.character(drvName), list(...)) : could not find function "Oracle" > Could anyone tell me how

Dear forumities, It's seem that there is no straight forward way to calculate R2 of a cluster solution in R. So, I would like to know if I'm right when calculating a R2-like statistic for a given clustering solution. In fact, I have different cluster solution for a given set of data. I would like to know which cluster solution gives the highest R2. My data (5 variables) are scaled to a 0

Hi all, I want to plot a facet plot with column names as x and column values as y. One plot for each row. here is part of my dataset: Gene T0h T0.25h T0.5h T1h T2h T3h T6h T12h T24h T48h NM_001001130 68 95 56 43 66 62 68 90 63 89 NM_001001144 0 1 4 0 1 1 1 4 1 2 NM_001001152 79 129 52 50 24 45 130 154 112 147 NM_001001160 1 1 2 0 0 1 0 0 0 1 NM_001001176 0 0 0 0 0 0 0 0 0 0 NM_001001177

hello folks, i am learning R and microarray analysis from scratch using different sites. today i am doing an exercise from http://manuals.bioinformatics.ucr.edu/home/R_BioCondManual#R_functions the section i am at is 2. Affymetrix data analysis. I understand the syntax given in this section up until: design <- model.matrix(~ -1+factor(c(1,1,2,2,3,3))) # Creates appropriate design

I have a data frame like this: T0h T0.25h T0.5h T1h T2h T3h T6h T12h T24h T48h C0h C0.25h C0.5h C1h C2h C3h C6h C12h C24h C48h NM_001001130 68 95 56 43 66 62 68 90 63 89 65 85 58 49 81 49 76 73 48 77 NM_001001152 79 129 52 50 24 45 130 154 112 147 56 85 67 33 52 31 93 77 65 109 ...

Dear R community, I'm trying to remove a loop from my code but I'm stock and I can't find a good way to do it. Hopefully one of you will have something clever to propose. Here is a simplified example: I have a squared matrix: > nom.plac2 <- c("102", "103", "301", "303","304", "403") > poids2 <- matrix(NA, 6,6,

Hello Arnaud, You posted this question a long long time ago, however I found your answer so I decided to post it anyway in case somebody else have the same problem as you and me. You were actually very close in finding your solution. The function DoWritedbf is an internal function from the foreign package. To access it outside of the package just do: foreign:::DoWritedbf so in your line:

Hello, First of all, thanks very much for this piece of software, it's working like a charm. I'm using it to boot a Linux distro I'm working on. What I want to deal with is the display of accented characters in ISOLINUX. First, I set EXTLINUX up to boot an ext2-formatted USB stick. No problem so far. I wrote a config file ordering to load a font (lat-9w) to display french

Hello R users, I have some troubles to modify the "per default" legend in the scatterplot.matrix using the car package. Here is my code: -------------------------------------------------------------------- scatterplot.matrix(~ a + b + c, groups=treatment, by.groups=F, data=newfile3, smooth=T, pch=c(16,4), cex=1.3,

Dear R users, Is someone know how to get a simple analysis of variance table using other random distribution than normal (ex: Poisson or Binomial)? When I try, I recieved this message: (See my R code below) --------------------- R response ---------------------------- Erreur dans anova(fit.poisson) : single argument anova for GLMMs not yet implemented OR Erreur dans anova(fit.binomial) :

In my var/log.smb file I am getting the error: String overflow by 18 in safe_strcopy [adaigle ahenry ajleblanc ] a long list My share is setup with the list of valid users i.e. valid users = user1 user2 \ user3 user4 writeable =yes Is there some limit to the number of users listed this way. A bad character in there? Some users are not getting their drives mapped.

Hi, I try to compute circular means for a matrix with NAs, but the function "circ.mean" return only means for lines with complete values and do not accept "na.omit=T" or "na.rm=T", or "na.action=na.omit", or "na.fail=T". Also, I try to use "polar.plot" of the package poltrix with the same matrix, but angles are not display because of

Dear all, I am trying to merge two data frames based on a common column but for this common column both data frame do not have the same length and associated information. I checked previous exemples in the list but was not able to apply them in my case... Is someone know how to do that? Below is my code with the expected result: # data frame 1 Id1 <- c(1,1,1,2,2,2,3,3,3) Habit1 <-

A monsieur Andrew Je suis ?tudiant en fin de cycle de formation en infomatique ? l'institut Africain d'informatique de Libreville.J'?tais sur vos sites et je suis persuad?s que vous ?tes la personne la bien plac?e pour m'orienter et m'aider sur un sujet concernant "l'authentification linux des sessions windows".Ainsi on me demande de concevoir et de mettre en

Hello everybody, I'm trying to create my own color palette on R, in order to interpolate some different temperature data on different maps (daily means, seasonal means,...). I would like to create a color palette which works for each map, so I need a color palette between -40 and +40?C. Sometimes my data for one map range from -10 to +20, sometimes from 10 to 30, etc... but always between

I presently have a Samba server (3.0.21b) set up as a member server in an NT4 domain (with a real Windows NT4 PDC). We are migrating to an Active Directory domain (with a real Windows 2003 domain controller). We have set up a two-way trust between the old NT4 domain "CLUNKY" and the new ADS domain "SLEEK" (aka sleek.local). The Samba server is a member of the CLUNKY domain

Hi, I am having troubles trying to install package "sp" in R (2.13.0) on mac OSX. I have tried installing the package using GUi or function install.packages but it didn't work. Here is the error message I get: also installing the dependency ‘rgdal’ trying URL 'http://cran.univ-lyon1.fr/src/contrib/rgdal_0.6-33.tar.gz' Content type 'application/x-gzip' length