Displaying 20 results from an estimated 10000 matches similar to: "Need to fit a regression line using orthogonal residuals"
2007 Jan 29
1
[Fwd: Need to fit a regression line using orthogonal residuals]
[Originally sent this to r-help at lists.R-projects.org, but in case that's
the wrong list I'm re-posting. Apologies if this becomes a re-post]
-------------- next part --------------
An embedded message was scrubbed...
From: Jonathon Kopecky <jkopecky at umich.edu>
Subject: Need to fit a regression line using orthogonal residuals
Date: Mon, 29 Jan 2007 14:52:24 -0500
Size: 1138
2005 Jul 07
1
Orthogonal Distance Regressions
Hi,
I work with environmental data and want to determine correlations between
variables that either have no "dependent/independent" relationship or the
relationship is unknown. Therefore I prefer to use orthogonal distance
regression (orthogonal linear regression, perpendicular sum of squares,
etc.). I am trying to get set up to do this in R, but the various
terminologies are
2008 Sep 01
1
intercept of 3D line? (Orthogonal regression)
I posted before recently about fitting 3D data x, y, z where all have
error attached.
I want to predict z from x and y; something like
z = b0 + b1*x + b2*y
But multiple regression is not suitable because all of x, y, and z have errors.
I have plotted a 3D scatterplot of some data using rgl. I see that the
data form a cigar-shaped cloud. I think multiple regression is only
suitable when the points
2005 Jun 14
2
ordinary polynomial coefficients from orthogonal polynomials?
How can ordinary polynomial coefficients be calculated
from an orthogonal polynomial fit?
I'm trying to do something like find a,b,c,d from
lm(billions ~ a+b*decade+c*decade^2+d*decade^3)
but that gives: "Error in eval(expr, envir, enclos) :
Object "a" not found"
> decade <- c(1950, 1960, 1970, 1980, 1990)
> billions <- c(3.5, 5, 7.5, 13, 40)
> #
2003 Aug 13
3
A question on orthogonal basis vectors
Hey, R-listers,
I have a question about determining the orthogonal
basis vectors.
In the d-dimensinonal space, if I already know
the first r orthogonal basis vectors, should I be
able to determine the remaining d-r orthognal basis
vectors automatically?
Or the answer is not unique?
Thanks for your attention.
Fred
2011 Feb 02
2
unequally spaced factor levels orthogonal polynomial contrasts coefficients trend analysis
Hello [R]-help
I am trying to find
> a package where you can do ANOVA based trend analysis on grouped data
> using orthogonal polynomial contrasts coefficients, for unequally
> spaced factor levels. The closest hit I've had is from this web site:
>(http://webcache.googleusercontent.com/search?q=cache:xN4K_KGuYGcJ:www.datavis.ca/sasmac/orpoly.html+Orthogonal+polynomial
>l
but I
2010 Dec 03
3
Checking for orthogonal contrasts
A common point made in discussion of contrasts, type I, II, III SS etc
is that for sensible comparisons one should use contrasts that are
'orthogonal in the row-basis of the model matrix' (to quote from
http://finzi.psych.upenn.edu/R/Rhelp02/archive/111550.html)
Question: How would one check, in R, that this is so for a particular
fitted linear model object?
Steve Ellison
2005 Apr 06
2
two methods for regression, two different results
Please forgive a straight stats question, and the informal notation.
let us say we wish to perform a liner regression:
y=b0 + b1*x + b2*z
There are two ways this can be done, the usual way, as a single
regression,
fit1<-lm(y~x+z)
or by doing two regressions. In the first regression we could have y as
the dependent variable and x as the independent variable
fit2<-lm(y~x).
The second
2004 Nov 02
2
Matrix decomposition: orthogonal complement
Hello,
How I can compute in R the orthogonal complement of one matrix?
If A (n x m ) matrix of full column rank (n>m), its orthogonal
complement is denoted by A_ .
A_ is n X (n-m) matrix of full column rank and such that A'A_=0.
I need to compute A_. How I can compute A_ in R?
Best Regards,
/Florin
--
Florin G. Maican
Ph.D. candidate, Department of Economics
School of
2006 Apr 20
1
Mutually Orthogonal Latin Squares
Hi all,
The package crossdes could contruct a complete sets of mutually orthogonal latin squares.
The construction works for prime powers only.
I hope to know whether there is a way to construct a mutually orthogonal Lation square for
10 or other numbers that could not be prime powers.
Thanks for any suggestions.
Best wishes,
Jinsong Zhao
2008 Jul 01
1
Orthogonal polynomials and poly
Dear All,
I have found in the poly help this sentence:
The orthogonal polynomial is summarized by the coefficients, which can be
used to evaluate it via the three-term recursion given in Kennedy & Gentle
(1980, pp. 343–4), and used in the predict part of the code.
My question: which type of orthogonal polynomials are used by this function?
Hrmite, legendre..
TIA
Giovanni
[[alternative HTML
2002 Oct 08
2
Orthogonal Polynomials
Looking to the wonderful statistical advice that this group can offer.
In behavioral science applications of stats, we are often introduced to
coefficients for orthogonal polynomials that are nice integers. For
instance, Kirk's experimental design book presents the following
coefficients for p=4:
Linear -3 -1 1 3
Quadratic 1 -1 -1 1
Cubic -1 3 -3 1
In R orthogonal
2004 Oct 30
2
How to plot PDF which is in the form of orthogonal polynomial
Dear all
using the orthogonal polymial on a set of data, I get an approximate
density which basically is in the form: exp(-polynomial),
as you know, the
parameters are the converged coeeficients.
obviously, It is hard, if not impossible, to use the inverse CDF method to
get a
sample and then plot density. then how can I plot the approximated density
in order to have a graphical comparision
2010 Mar 22
2
problem with abline and lines
Dear R users,
I need to plot to perpendicular straight lines. However, although I set the coefficients
so that the lines are perpendicular, they do not look to be so in the plot. Here is a minimal working example:
plot(x=c(-1, 1), y=c(-1, 1)); abline(a=0, b=1/sqrt(2)); abline(a=0, b=-1/sqrt(2))
Please tell me if the same problem is valid by you. I am running R-2.10.1 on Linux.
Is there a way
2010 Aug 23
4
Draw a perpendicular line?
Hi,
I am trying to draw a perpendicular line from a point to two points.
Mathematically I know how to do it, but to program it, I encounter some
problem and hope can get help. Thanks.
I have points, A, B and C. I calculate the slope and intercept for line
drawn between A and B.
I am trying to check whether I can draw a perpendicular line from C to line
AB and get the x,y value for the point D
2004 May 06
5
Orthogonal Polynomial Regression Parameter Estimation
Dear all,
Can any one tell me how can i perform Orthogonal
Polynomial Regression parameter estimation in R?
--------------------------------------------
Here is an "Orthogonal Polynomial" Regression problem
collected from Draper, Smith(1981), page 269. Note
that only value of alpha0 (intercept term) and signs
of each estimate match with the result obtained from
coef(orth.fit). What
2005 Jul 08
1
Orthogonal regression
Dear R-Users,
is there any statement to fit a orthogonal regression
in R environment?
Many thanks in advance.
Best regards,
Vito
Diventare costruttori di soluzioni
Became solutions' constructors
"The business of the statistician is to catalyze
the scientific learning process."
George E. P. Box
"Statistical thinking will one day be as necessary for efficient
1999 Dec 10
1
orthogonal and nested model
I'm working with a orthogonal and nested model (mixed).
I have four factors, A,B,C,D;
A and B are fixed and orthogonal
C is nested in AB interaction
and finally, D is nested in C.
I would like to model the following
Y_ijklm=Mu+A_i+B_j+AB_ij+C_k(ij)+D_l(k(ij))+Error_m(...)
I used the next command
>summary(aov(abund~A*B + C % in % A:B + D % in % C % in % A:B ,datos))
Is it the correct
2005 Apr 28
1
standard errors for orthogonal linear regression
Could someone please help me by giving me a reference to how one computes standard errors for the coefficients in an orthogonal linear regression, or perhaps someone has some R code? (I would accept a derivation or formula, but as a former teacher, I know how that can rankle.) I tried to imitate what's done in the code for lm() but went astray somewhere and got nonsense.
(This type of
2003 Jan 21
1
Orders of terms in formulae
Hi,
Given that R reports Type I sums of squares, isn't it a bit anachronistic
that it re-orders terms in formulae?
> d <- expand.grid(y=rnorm(8),
+ A=factor(c(1,2)),
+ B=factor(c(1,2)),
+ C=factor(c(1,2)))
> summary(aov(y ~ A+B+A:B+C,data=d))
Df Sum Sq Mean Sq F value Pr(>F)
A 1 8.294e-34 8.294e-34 1.027e-33 1