Re: [R] nls and factor

2006-04-25 Thread Manuel Gutierrez
: On Thu, 20 Apr 2006, Manuel Gutierrez wrote: Is it possible to include a factor in an nls formula? Yes. What do you intend by it? If you mean what it would mean for a lm formula, you need A[a] and starting values for A. There's an example on p.219 of MASS4. I've searched the help

[R] nls and factor

2006-04-20 Thread Manuel Gutierrez
Is it possible to include a factor in an nls formula? I've searched the help pages without any luck so I guess it is not feasible. I've given it a few attempts without luck getting the message: + not meaningful for factors in: Ops.factor(independ^EE, a) This is a toy example, my realworld case is

Re: [R] nls and factor

2006-04-20 Thread Manuel Gutierrez
, Manuel Gutierrez wrote: Is it possible to include a factor in an nls formula? I've searched the help pages without any luck so I guess it is not feasible. I've given it a few attempts without luck getting the message: + not meaningful for factors in: Ops.factor(independ^EE

Re: [R] se.fit in predict.nls

2006-01-19 Thread Manuel Gutierrez
: 0x88a076c But I can not find what is form. Any help, please. Manuel Manuel Gutierrez wrote: The option se.fit in predict.nls is currently ignored. Is there any other function available to calculate the error in the predictions? Thanks, Manuel

[R] se.fit in predict.nls

2006-01-18 Thread Manuel Gutierrez
The option se.fit in predict.nls is currently ignored. Is there any other function available to calculate the error in the predictions? Thanks, Manuel __ LLama Gratis a cualquier PC del Mundo. Llamadas a fijos y móviles desde 1

[R] error propagation

2005-12-29 Thread Manuel Gutierrez
Are there any functions to do error propagation in R? I have done a search with little success. Any pointers to read about this topic are greatly welcomed. My specific problem is: I use a linear model (lm) to predict the biomass of an individual in a population, then I add up the biomass of all

[R] addition of error terms

2005-12-29 Thread Manuel Gutierrez
Are there any functions to do error propagation in R? I have done a search with little success. Any pointers to read about this topic are greatly welcomed. My specific problem is: I use a linear model (lm) to predict the biomass of an individual in a population, then I add up the biomass of all

[R] dealing with multicollinearity

2005-04-11 Thread Manuel Gutierrez
I have a linear model y~x1+x2 of some data where the coefficient for x1 is higher than I would have expected from theory (0.7 vs 0.88) I wondered whether this would be an artifact due to x1 and x2 being correlated despite that the variance inflation factor is not too high (1.065): I used

Re: [R] Re: testing slopes different than a given value

2005-02-11 Thread Manuel Gutierrez
Thanks to all, Yes, I meant a single test for both coefficients. Peter's reply is what I wanted. I've tried with linear.hypothesis but I must confess that with my limited statistical experience and without the car book at hand, the nomenclature for the function was a bit difficult to understand

[R] testing slopes different than a given value

2005-02-10 Thread Manuel Gutierrez
In a multiple linear regression with two independent variables is there any function in R to test for the coefficients being different than some given values? Example: x1-rnorm(100) x2-rnorm(100) y-3+0.6*x1+0.3*x2 lm(y~x1+x2) Obtain a test for the coefficients for x1 being different than 0.6 and

[R] Rare Cases and SOM

2005-02-04 Thread Manuel Gutierrez
I am trying to understand how the SOM algorithm works using library(class) SOM function. I have a 1000*10 matrix and I want to be able to summarize the different types of 10-element vectors. In my real world case it is likely that most of the 1000 values are of one kind the rest of other (this is

[R] memory and swap space in ncdf

2005-01-21 Thread Manuel Gutierrez
I've a linux system with 2Gb of memory which is not enough for reading a 446Mb netcdf file using ncdf: library(ncdf) ncold - open.ncdf(gridone.grd) Error: cannot allocate vector of size 1822753 Kb When I look at the free memory in my system I can see that none of the Swap space is being used by

[R] unlist kills R

2005-01-05 Thread Manuel Gutierrez
When I try to unlist a very large list, R is killed without any other warning: A-as.list(as.data.frame(matrix(1:21639744,nrow=3578,ncol=6048))) with AA-unlist(A) or AA-c(A,recursive=TRUE) I get a R terminado (killed) and the end of the session. I think I'll need to get more RAM (now 1Gb, any

[R] test multiple objects for being equal length

2004-12-09 Thread Manuel Gutierrez
I could not find any help pages on How to test many objects for being of equal length Something like identical for more than two objects? x-1:6 y-1:10 z-3:5 ## For two objects I can do: identical(length(x),length(y)) ## For more than two I currently can do:

[R] substitute accents

2004-11-25 Thread Manuel Gutierrez
I have an openoffice spreadsheet with a column of character strings. Some of them contain accents. I want to read it in R so I have saved it as a csv file using Western Europe (ISO-8859-1) character set (the default, I've tried other sets but it doesn't help). R reads it fine with

Re: [R] substitute accents

2004-11-25 Thread Manuel Gutierrez
of a UTF-8 locale, which R does not currently support and which some Linux distros have made their default. However, R does issue a warning -- so did you get one? On Thu, 25 Nov 2004, Manuel Gutierrez wrote: I have an openoffice spreadsheet with a column of character strings. Some