Displaying 20 results from an estimated 34 matches for "soybeans".
Did you mean:
soybean
2011 Feb 08
4
Interactions in a nls model
I am interested in testing two similar nls models to determine if the
lines are statistically different when fitted with two different data
sets; one corn, another soybean. I know I can do this in linear models
by testing for interactions. See Introductory Statistics with R by
Dallgaard p212-218 for an example. I have two different data sets I am
comparing to lai. ci.re should have very
1999 Nov 25
1
gnls
Doug,
I have been attempting to learn a little bit about nlme without too
much documentation except the online help. The Latex file in the nlme
directory looks interesting but uses packages that I do not have so
that I have not been able to read it.
I have run the example from gnls to compare it with the results I
get from my libraries (code below - I have not included output as it
is rather
2003 Apr 02
4
randomForests predict problem
Hello everybody,
I'm testing the randomForest package in order to do some simulations and I
get some trouble with the prediction of new values. The random forest
computation is fine but each time I try to predict values with the newly
created object, I get an error message. I thought I was because NA values
in the dataframe, but I cleaned them and still got the same error. What am
I
2010 May 11
1
merge two data frames
...ul/10",
"COPPER May/10", "CORN Jul/10", "CORN May/10", "COTTON NO.2 Jul/10",
"CRUDE OIL miNY May/10", "GOLD Jun/10", "HENRY HUB NATURAL GAS May/10",
"ROBUSTA COFFEE (10) Jul/10", "SILVER May/10", "SOYBEANS Jul/10",
"SPCL HIGH GRADE ZINC USD", "STANDARD LEAD USD", "SUGAR NO.11 Jul/10",
"SUGAR NO.11 May/10", "WHEAT Jul/10", "WHEAT May/10"), class = "factor"),
pl = c(4.75000000000003, -14, -1.5, 7.75, 3.74999999999999,...
2001 Jan 17
1
Pinheiro/Bates Soybean nlme failure
Dear Mixed Effect Friends,
Somehow, R(1021, Windows) seem to run differently from S Plus:
The soybean example from Pinheiro/Bates on page 290 fails
in R. (Soybean1 is Soybean with the NA and "critical" case
removed. Same procedure with full Soybean).
> fm1Soy.lis<-nlsList(weight~SSlogis(Time,Asym,xmid,scal),data=Soybean1)
> fm1Soy.nlme<-nlme(fm1Soy.lis)
Error: Singularity
2003 Apr 03
0
debugging question
Dear R-devel,
A user reported a strange problem with predict.randomForest in the
randomForest package yesterday, and I'm baffled by it. The code at the end
of the message produces the error. The problem is that, in
predict.randomForest, there's a .Fortran call to "runforest". One of the
arguments passed in is "countts", which is a vector of doubles. The error
2010 Apr 23
2
read.csv data frame thousands separator
Dear group,
Here is my df, trades1 :
trades1 <-
structure(list(Instrument.Long.Name = c("CORN", "CORN", "CORN",
"CORN", "CORN", "SOYBEANS", "SOYBEANS", "SOYBEANS", "SOYBEANS",
"SOYBEANS", "SOYBEANS", "STANDARD LEAD USD", "STANDARD LEAD USD",
"SPCL HIGH GRADE ZINC USD", "SPCL HIGH GRADE ZINC USD", "SPCL HIGH GRADE
ZINC USD",...
2004 Feb 20
1
nlme and multiple comparisons
This is only partly a question about R, as I am not quite sure about the
underlying statistical theory either.
I have fitted a non-linear mixed-effects model with nlme. In the fixed
part of the model I have a factor with three levels as explanatory
variable. I would like to use Tukey HSD or a similar test to test for
differences between these three levels.
I have two grouping factors:
2010 Oct 20
0
bootstrapping nonlinear mixed effects models
In order to bootstrap nonlinear regression, the following code works.
library(nlme)
data(Soybean)
fm1.nls <- nls(weight ~ SSlogis(Time, a, b, c), data=Soybean)
summary(fm1.nls)
bstat <- function(A, indices) {
mboot <- nls(weight ~ SSlogis(Time, a, b, c), data=Soybean[indices, ])
return(coef(mboot))
}
require(boot)
boot1 <- boot(Soybean, bstat, 2000)
boot1
For nonlinear mixed
2010 May 11
1
merging data frame
...merge.
Here they are:
pose16 <-
structure(list(DESCRIPTION = structure(c(1L, 2L, 3L, 4L, 5L,
6L, 7L, 8L, 9L, 12L), .Label = c("COPPER May/10", "COTTON NO.2 Jul/10",
"CRUDE OIL miNY May/10", "GOLD Jun/10", "ROBUSTA COFFEE (10) Jul/10",
"SOYBEANS Jul/10", "SUGAR NO.11 Jul/10", "SUGAR NO.11 May/10",
"WHEAT Jul/10", "PRIMARY NICKEL USD", "PRM HGH GD ALUMINIUM USD",
"SPCL HIGH GRADE ZINC USD", "STANDARD LEAD USD"), class = "factor"),
POSITION = c(2, -3,...
2010 Apr 28
1
gsub,regex and replacing
...gsub and the
regex functionalities to do this.
Here is the data frame (df):
DESCRIPTION prix quantity
1 CORN Jul/10 -1.5 0
2 CORN May/10 -1082.0 -3
3 ROBUSTA COFFEE (10) Jul/10 11084.0 8
4 SOYBEANS Jul/10 1983.5 2
5 SPCL HIGH GRADE ZINC USD Jul/10 -2464.0 -1
6 STANDARD LEAD USD Aug/10 -118.0 0
For each df$DESCRIPTION element containing "USD" I want to remove the last
part of it(i.e Jul/10, or Aug/10...).
I was thinking of something like that:
>Df$D...
2002 May 08
3
R 1.5.0 terminated while plotting groupedData
Dear R-help;
Following the installation of R 1.5.0 and library packages to my Win 2K
machine, I tested the new R system:
> library(nlme)
Loading required package: nls
> library(lattice)
Loading required package: grid
>data(Soybean)
>plot(Soybean)
A program error appears: "Rgui.exe has generated errors and will be closed
by Windows. You will need to restart the program. An
2004 Oct 01
4
gnls or nlme : how to obtain confidence intervals of fitted values
Hi
I use gnls to fit non linear models of the form y = alpha * x**beta
(alpha and beta being linear functions of a 2nd regressor z i.e.
alpha=a1+a2*z and beta=b1+b2*z) with variance function
varPower(fitted(.)) which sounds correct for the data set I use.
My purpose is to use the fitted models for predictions with other sets
of regressors x, z than those used in fitting. I therefore need to
2010 Aug 04
3
merge two data frames
...data frames:
av <-
structure(list(DESCRIPTION = c("COFFEE C Sep/10", "COPPER Sep/10",
"CORN Dec/10", "CRUDE OIL miNY Sep/10", "GOLD Aug/10", "HENRY HUB
NATURAL GAS Sep/10",
"PALLADIUM Sep/10", "SILVER Sep/10", "SOYBEANS Nov/10", "SPCL HIGH
GRADE ZINC USD",
"SUGAR NO.11 Oct/10", "WHEAT Sep/10"), prix = c(-168.3, -1.60000000000002,
-773.75, -78.75, -1168.3, -0.0919999999999996, -470.75, 1758.5,
-975.25, 1964, -19.09, -605.75), pos = c(-1, 0, -2, -1, -1, 0,
-1, 1, -1, 1, -1, -1), PL...
2000 Feb 11
1
R CMD check [nlme|MASS] fails (PR#431)
Mmmh, seems as if I really should change my options as I seem to keep
sending off empty bug-reports ;-/ Sorry guys. Here is the content
that should have been in the last e-mail:
`R CMD check nlme' fails on my machine. The final output in
nlme-Ex.Rout is:
> library(nlme)
> data(Soybean)
> fm1 <- nlme(weight ~ SSlogis(Time, Asym, xmid, scal), data = Soybean,
+ fixed =
2001 May 31
1
nlme and memory
I am trying to follow an example in Pineiro & Bates's book
library(nlme)
data(Soybean)
soy.lis<-nlsList(weight~SSlogis(Time,Asym,xmid,scal),data=Soybean)
soy.nlme<-nlme(soy.lis)
Error: Calloc could not allocate (151974496 of 8) memory
This is not a large problem-- only 412 observations. I am using R-1.2.3
with automatic memory allocation. Does this message mean that I need to
run
2009 May 04
1
how to change nlme() contrast parametrization?
How to set the nlme() function to return the answer without the intercept parametrization?
#=========================================================================================
library(nlme)
Soybean[1:3, ]
(fm1Soy.lis <- nlsList(weight ~ SSlogis(Time, Asym, xmid, scal),
data = Soybean))
(fm1Soy.nlme <- nlme(fm1Soy.lis))
fm2Soy.nlme <- update(fm1Soy.nlme,
2005 Mar 02
1
Using varPower in gnls, an answer of sorts.
Back on January 16, a message on R-help from Ravi Varadhan described a
problem with gnls using weights=varPower(). The problem was that the
fit failed with error
Error in eval(expr, envir, enclos) : Object "." not found
I can reliably get this error in version 2.0.1-patched 2004-12-09 on
Windows XP and 2.0.1-Patched 2005-01-26 on Linux.
The key feature of that example is that the
2010 Apr 28
3
data frame move columns
Dear group,
Here is my df:
avprix <-
structure(list(DESCRIPTION = c("CORN Jul/10", "CORN May/10",
"ROBUSTA COFFEE (10) Jul/10", "SOYBEANS Jul/10", "SPCL HIGH GRADE ZINC USD",
"STANDARD LEAD USD"), prix = c(-1.5, -1082, 11084, 1983.5, -2464,
-118), quantity = c(0, -3, 8, 2, -1, 0), fees = c(-8.64, -30.24,
-26.4, -25.92, -37.5, -25)), .Names = c("DESCRIPTION", "prix",
"quantity"...
2009 May 01
2
combine column names
Dear R-list,
Is there any way I can extract part or whole column name from two
variables and combine them? Not values, but column names.
For example,
I have data set with 70 variables. The column names are "Alfalfa",
"Corn","Soybean",..."Sunflower". After I combine the two variable (say
Alfalfa and Corn), I want to rename it as "AlfCor" or