similar to: plot() and lines() multivariate problem

Dear list, I have a problem with plotting multivariate data and models. I suppose the answer is obvious for experienced users but I didn't find it looking in usual documentation or in the mail archive... Let's say I use a dataset like : a <- sort(round(rnorm(90),3)) niv <- factor(c(rep('B',30), rep('E', 30),

dear list, I'm a newbie with R but it seems very interesting to me. First, to train myself and students, I'd like to use it for teaching labs next year. one of my labs deals with the barn owl predation in my region. tudents have to work on a contingency table tytogrp : station1 station2 station3 ... n1,1 n1,2 n1,3 ... n2,1 n2,2 n2,3 ...

Hi useRs, I'm trying for the first time to use a sem. The model finally runs, but gives a warning saying : "In sem.default(ram = ram, S = S, N = N, param.names = pars, var.names = vars, : Could not compute QR decomposition of Hessian. Optimization probably did not converge. " I found in R-help some posts on this warning, but my attemps to modify the code didn't change

Dear list member, My question is related to input file format to an Anova from car package. Here is an example of what I did: My file format is like this (and I dislike the idea that I will need to recode it): Hormone day Block Treatment Plant Diameter High N.Leaves SH 23 1 1 1 3.19 25.3 2 SH 23 1 1 2 3.42 5.5 1 SH 23 1 2 1 2.19 5.2 2 SH 23 1 2 2 2.17 7.6 2 CH 23 1 1 1 3.64 6.5 2 CH 23 1 1 2

I am trying to define groupings from levels of factor variables and this the warning message that R give "& not meaningful for factors". The nature of my task is this. I have a variable stage which has the levels (1B, 2A, 2B) - these are the AJCC TNM stages of cancer, and another variable diameter with factor levels ("=< 4", "4 - 6.5, > 6.5; limit values are

Hi, I have a function that generates a set of data but I am having problems determining the parameters using the nls fitting procedure. #### "MH"<-function(field,diameter,mu=10e-7,sig=0.1,Ms=100,chi=0){ #variables mu, sig, chi, Ms #input: field and diameter #all in CGS rho <- 5 kb <- 1.38e-16 t <- 300 length.d<-length(diameter) length.H<-length(field)

Regards, I have lately experienced a rather peculiar problem, details follows: Trying to minimise the number of samba servers, I have compiled in support for automounting in a Solaris 8 version of samba 3.0.14a. Amongst others, this server mounts a filesystem from a Solaris 6 box by request. The samba server mounts the Solaris 6 filesystem fine, but when the user tries to access the directory

Regards, Sent this one before, here's another go: My sambaserver, let's call it sserv, is running Solaris 8 (latest Recommended set) or 9 (tried both) compiled with automount support. When PC client A tries to access a file in a automounted mapping from sserv, the client is freezing and Windows Explorer has to be restarted. The NFS server which the samba server automounts from is a

Dear all, The 'make check' step fails for the pacakge mva on IBM AIX. The tail of the Rout log file looks like: > for(factors in 2:4) print(update(Harman23.FA, factors = factors)) Call: factanal(factors = factors, covmat = Harman23.cor) Uniquenesses: height arm.span forearm lower.leg weight 0.170 0.107 0.166

Regards, Samba 3.0.14a, running on Solaris 8 or Sgi IRIX 6.5.27. Scenario: 1. User has a symlink on his homeshare to a directory, let's call it dir1, which is automounted to the samba server (tried both Solaris samba server and Sgi, same result). 2. In a subdirectory under dir1 there's another directory (dir2) which contains a file (filename.wrl). 3. When the user directs Windows

Hello, I want to fit a Gompertz model for tree diameter growth that depends on a 4 levels edaphic factor (?Drain?) and I don?t manage to introduce the factor variable in the formula. Dinc is the annual diameter increment and D is the Diameter. >treestab > Dinc D Drain [1,] 0.03 26.10 2 [2,] 0.04 13.05 1 [3,] 0.00 24.83 1 [4,] 0.00 15.92 4

Hi, Maybe this is more a programming questions than a specific R-project question, but maybe there is someone who can point me in the right direction. I have a picture of cirkels which I took with a digital camera. Now I want to use the diameter of the cirkels on the picture for analysis in R. I can use pixmap to import the picture, but how do I find the outside cirkels and calculate the

I'm experiencing some confusion with the ellipsis argument (...). In my .Rprofile, I have the following functions: stderr <- function(x, ...) { sqrt( var(x, ...) / length(x) ) } se <- stderr I can use tapply to calculate some means: > tapply( Diameter, factor(Region), mean, na.rm=TRUE ) I II III IV V 0.02896429

I have a data set with the following structure (with many more obs.): var1 expt day diameter 1 1 2 0.5 1 1 3 0.9 1 1 4 1.3 1 1 5 1.7 1 2 2 0.3 1 2 3 0.5 1 2 4 0.9 1 2 5 1.6 2 1 2 0.7 2 1 3 1.2 2 1 4 1.6 2 1 5 2.3 2 2 2 0.4 2 2 3 0.8 2 2 4 1.6 2 2 5 3.2 I can get separate regression analysis for each level of var1 and expt with the command: by(data3, data3$var1:data3$expt, function(x)

Hi all, I would like to test the homoegeneity of variances between several linear model for some analysis of covariance. It seems that the Bartlett's test is a good test to use but I am having problem using with linear model and I cannot find any examples on the internet. There are some examples for comparisons of variances but not linear models. If I take the hellung data set, which is

Dear All, I wish to calculate the Gini index (ineq from same package) and some other indices for the diameter distribution of each plot (df dgtot). dgtot: IDPlot Diameter(cm) 1 4 34.0 2 4 23.0 3 4 38.0 ... 51 5 16.0 52 5 8.0 53 5 9.0 ... 5301 140 25.0 5302 140 12.0 5303 140 7.0 I use: >

Hi there, Am try to find tree diameter increment in a hierarchical data with tree diameter measurements, date of measurements, in plots that are located in different region. How do I do it? am very new to R. Thanks for you help. Maina Kariuki Level 3, 8 Nicholson Street East Melbourne VIC 3002 Ph: (03) 9637 9809 Email: maina.kariuki@dse.vic.gov.au Notice: This email and any

Hi everyone, I intend to do a discriminant analyse for 2 measures(eye diameter and body length) and for different areas to show differences between those areas if there are any. The raw data (eye diameter, body length) make one cloud of points so it seems there aren't any differences between those areas. If I take the mean values (eye diameter) for length classes, I nearly get linear functions

How do I create a legend with ggplot? I think should be getting the FuelTypeNum in the legend. Plot: http://r.789695.n4.nabble.com/file/n4632471/Rplot.jpeg My code is: ggplot(data=tempTable, aes(x=Bands8, y=SubPercent, fill=FuelTypeNum)) + geom_bar(position="stack", stat="identity") + scale_colour_hue('my legend', breaks = levels(tempTable$FuelTypeNum),

Hi there, I have been having trouble running negative binomial regression (glm.nb) using library MASS in R v2.15.0 on Mac OSX. I am running multiple models on the variables influencing the group size of damselfish in coral reefs (count data). For total group size and two of my species, glm.nb is working great to deal with overdispersion in my count data. For two of my species, I am getting a