Displaying 20 results from an estimated 6000 matches similar to: "Solution to non-linear equation problem"
2005 Nov 30
8
Solving Systems of Non-linear equations
I am trying to write a function that will solve a simple system of
nonlinear equations for the parameters that describe the beta
distribution (a,b) given the mean and variance.
mean = a/(a+b)
variance = (a*b)/(((a+b)^2) * (a+b+1))
Any help as to where to start would be welcome.
--
Scott Story
Graduate Student
MSU Ecology Department
319 Lewis Hall
Bozeman, Mt 59717
406.994.2670
sstory at
2005 Nov 30
1
Loop within nlme
I am trying to mimic the SAS code below in R. The trick is that each
row in the dataset has variable "t" which controls how many times the
do-loop below will be iterated (that is, the model is fit to the
response, ifate, 0 to t-1 times for each row of data). Is it possible to
incorporate a loop like this into nlme by writing a function? Can
anybody provide some hints to get me on my
2003 Nov 17
2
Symbolic math?
Hi Folks,
I am using Windows 2000 and was wondering what (Open Source) software R
users use or might recommend for symbolic computations (aside from the ol'
noggin, e.g., Maxima, Mathomatic) .
Thanks,
Hank
Dr. Martin Henry H. Stevens, Assistant Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056
Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
2003 Dec 24
0
Solution to "Can anyone help me reproduce this SAS Mixed output??"
To those who might be interested -- following is the solution to my
previous post regarding reproducing output from SAS Proc Mixed for a
two-factor crossed random effects ANOVA model.
I am graciously endebted to the kind replys from two statisticians for
this solution whose names I will refrain from mentioning for the sake of
privacy.
I hope this helps someone?!
-- Phil Turk
> hw7 <-
2008 Jun 13
2
Quartile regression question
I have data that looks like
lake,loglength,logweight
1,2.369215857,1.929418926
1,2.426511261,2.230448921
1,2.434568904,2.298853076
1,2.437750563,2.298853076
1,2.442479769,2.230448921
1,2.445604203,2.356025857
...
102,2.722633923,3.310268367
102,2.781755375,3.502153893
102,2.836324116,3.683407299
102,2.802773725,3.583312152
102,2.790285164,3.546419267
102,2.806179974,3.599118565
2008 Aug 19
1
Histogram binning
I am trying to produce frequencies in defined intervals however I can't seem
to figure out how to get R to bin my data the way I want it to.
I have several thousand lengths of fish that I want to be binned as follows:
Ex.
Length Bin
209 200
219 210
431 430
727 720
That is, bins with any length equal to or greater than the lower
2007 Jun 11
1
Rearranging Capture History Data in R
What code can i use to convert a table like this:
Tag# Date
1 1
2 1
3 1
4 1
2 2
4 2
1 3
2 3
4 4
Into one like this:
Tag 1 2 3 4 #Date header
1 1 0 0 1
2 1 1 1 0
3 1 0 0 0
4 1 1 0 1
Thanks,
Ben Cox
Research Assistant (M.S.)
Montana Cooperative
2005 Jul 19
1
deriv - accessing numeric output listed under gradient attribute
Hi,
I am interested in using the numeric output from the "gradient" attribute of
deriv's output in subsequent analyses.
But, I have so far been unable to determine how to do so.
I will use the example from the deriv help to illustrate.
> ## function with defaulted arguments:
> (fx <- deriv(y ~ b0 + b1 * 2^(-x/th), c("b0", "b1", "th"),
2010 Mar 26
1
Trouble loading package
I am trying to load a package called Rmark, but when I run
library(Rmark)
I get the following:
> library(RMark)
Error in !character.only : invalid argument type
Error in library(RMark) : .First.lib failed for 'RMark'
When I try to load Rmark from the packages menu, I get:
> local({pkg <- select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg,
2010 Aug 23
2
Quantile Regression and Goodness of Fit
All -
Does anyone know if there is a method to calculate a goodness-of-fit
statistic for quantile regressions with package quantreg?
Specifically, I'm wondering if anyone has implemented the
goodness-of-fit process developed by Koenker and Machado (1999) for R?
Though I have used package quantreg in the past, I may have overlooked
this function, if it is included.
Citation:
Koenker, R. and
2009 Mar 14
1
dispcrepancy between aov F test and tukey contrasts results with mixed effects model
Hello,
I have some conflicting output from an aov summary and tukey contrasts
with a mixed effects model I was hoping someone could clarify. I am
comparing the abundance of a species across three willow stand types.
Since I have 2 or 3 sites within a habitat I have included site as a
random effect in the lme model. My confusion is that the F test given by
aov(model) indicates there is no
2008 Jul 03
1
lm() question
I have data that looks like
YC Age Num
82 11 2
83 10 0
84 9 8
85 8 21
86 7 49
87 6 18
88 5 79
89 4 28
90 3 273
91 2 175
with a program
mod1=lm(log(Num+1)~YC, data=box44)
plot(log(Num+1)~YC, data=box44, pch=19, xlab="Year Class",
ylab="Loge Number at age", ylim=c(0,6), xlim=c(91,82))
abline(lm(log(Num+1)~YC), col="blue", lwd=2)
summary(mod1)
I need to
2008 Jun 12
2
Predicting from an nls model
I keep running up against the same error when I try to plot a line from a nls model. The data is fisheries length/weight data. Code follows:
require(graphics)
pow = nls(Weight~alpha*Length^beta, data=wae,
start=list(alpha=0.0000001, beta=3.0), trace=TRUE)
predict(pow)
plot(Weight~Length,
2003 Feb 07
2
Data manipulation
I am interested in building a model with a subset of data from a column.
The first 6 lines of my data look like this:
QUAD YEAR SITE TREAT HERB TILL PLANT SEED Kweed
1 A4 2002 s 1 N N N N 55.00
2 A10 2002 s 1 N N N N 60.00
3 B2 2002 s 1 N N N N 35.00
4 C2 2002 s 1 N N N N 23.00
5 C9
2008 Dec 18
2
Contextstack overlow
All -
I have a number of rows that I am assigning length classes to via
l.class<-with(wae,
ifelse((Length>=120)&(Length<130),"125",
ifelse((Length>=130)&(Length<140),"135",
ifelse((Length>=140)&(Length<150),"145",
ifelse((Length>=150)&(Length<160),"155",
2011 May 22
2
Finding solution set of system of linear equations.
I have a simple system of linear equations to solve for X, aX=b:
> a
[,1] [,2] [,3] [,4]
[1,] 1 2 1 1
[2,] 3 0 0 4
[3,] 1 -4 -2 -2
[4,] 0 0 0 0
> b
[,1]
[1,] 0
[2,] 2
[3,] 2
[4,] 0
(This is ex Ch1, 2.2 of Artin, Algebra).
So, 3 eqs in 4 unknowns. One can easily use row-reductions to find a
homogeneous solution(b=0) of:
X_1
2008 Aug 22
2
Newbie programming help
All -
Not sure if this is a real programming question, but here goes:
I have data that looks like
Lake Length Weight
1 158 45
1 179 70
1 200 125
1 202 150
1 206 145
1 209 165
1 210 140
1 215 175
1 216 152
1 220 150
1 221 165
...
where lake goes from 1 - 84 and the number of rows for each lake is variable (but > ~20).
I'm trying to do two things: 1) build a simple linear model of the
2008 Sep 02
2
nls.control()
All -
I have data:
TL age
388 4
418 4
438 4
428 5
539 10
432 4
444 7
421 4
438 4
419 4
463 6
423 4
...
[truncated]
and I'm trying to fit a simple Von Bertalanffy growth curve with program:
#Creates a Von Bertalanffy growth model
VonB=nls(TL~Linf*(1-exp(-k*(age-t0))), data=box5.4,
start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE)
#Scatterplot of the data
plot(TL~age, data=box5.4,
2008 Jun 10
3
newbie nls question
I'm tyring to fit a relatively simple nls model to some data, but keep coming up against the same error (code follows):
Oto=nls(Otolith ~ Linf*(1-exp(-k(AGE-to))),
data = ages,
start = list(Linf=1000, k=0.1, to=0.1),
trace = TRUE)
The error message I keep getting is "Error in eval(expr, envir, enclos) : could not find function "k"". I've used this
2006 Apr 06
0
algebra
hello all. i have been trying to develop a representation (in the S4 sense) for a floating cash object, which would store a cash amount as a function of an arbitrary number of variables (with unknown values). for example, an interest rate swap may call for a payment in one year that can be represented as a function of a 3-month libor rate to be determined in nine months. this floating cash