Displaying 20 results from an estimated 54 matches for "scherber".
2007 Nov 08
6
Extract correlations from a matrix
...0)
w=cor(m,use="pairwise.complete.obs")
How do I extract only those correlations that are >0.6?
w[w>0.6] #obviously doesn?t work,
and I can?t find a way around it.
I would very much appreciate any help!
Best wishes
Christoph
(using R 2.5.1 on Windows XP)
--
Dr. Christoph Scherber
DNPW, Agroecology
University of Goettingen
Waldweg 26
D-37073 Goettingen
Germany
phone +49(0)551 39 8807
fax +49(0)551 39 8806
homepage www.gwdg.de/~cscherb1
2005 Jan 07
6
coercing columns
Dear all,
I have a data frame that looks like this:
c1 c2 c3
A B C
B C A
A A B
and so on;
I?d like to produce one single vector consisting of the columns c1,c2,
c3, such that
vector=("A","B","A","B","C","A","C","A","B")
I guess it?s easy to do but I don?t know how...Can anyone
2008 Nov 06
2
replacing characters in formulae / models
...y[1]
Suppose I would like to replace the "x" by "w" in the formula object "y".
How can this be done? Somehow, the methods that can be used in character objects do not work 1:1 in
formula objects...
Many thanks and best wishes
Christoph
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2008 Feb 22
3
Simultaneously summarizing many models
...arks instead of m1 m2
Do you know a way around this? I guess I would need to use expression() or substitute() but I don?t
know how to construct m1, m2... without the quotation marks.
Thanks very much for your help!
Best wishes,
Christoph
(using R 2.5.1 on Windows XP SP2)
--
Dr. Christoph Scherber
DNPW, Agroecology
University of Goettingen
Waldweg 26
D-37073 Goettingen
Germany
phone +49(0)551 39 8807
fax +49(0)551 39 8806
homepage www.gwdg.de/~cscherb1
2008 Nov 25
4
glm or transformation of the response?
...ow would you deal with this kind of analysis? What would be your recommendation to the students,
given the fact that most of the standard glm models obviously don?t seem to produce good fits here?
Many thanks and best wishes
Christoph
(using R 2.8.0 on Windows XP)
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2005 May 02
4
"apply" question
Dear R users,
I??ve got a simple question but somehow I can??t find the solution:
I have a data frame with columns 1-5 containing one set of integer
values, and columns 6-10 containing another set of integer values.
Columns 6-10 contain NA??s at some places.
I now want to calculate
(1) the number of values in each row of columns 6-10 that were NA??s
(2) the sum of all values on columns 1-5
2004 May 04
1
RE: more on lm(y~x) question: removing NA´s
1. If your code actually runs, you should upgrade R, and quit using `_' for
assignment... 8-)
2. You seem to have an extraneous `]' after the na.exclude. Could that be
the problem?
Andy
> From: Christoph Scherber
>
> actually, the situation is much more complicated. I am producing
> multiple graphs within a "for" loop. For some strange reason, the
> plotting routine always stops once lm(y~x) encounters more than one
> missing value (I have marked the important bit with "**...
2008 Nov 03
2
standard errors for predict.nls?
Dear all,
Is there a way to retrieve standard errors from nls models? The help page tells me that arguments
such as se.fit are ignored...
Many thanks and best wishes
Christoph
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2005 Jan 24
4
lme and varFunc()
Dear R users,
I am currently analyzing a dataset using lme(). The model I use has the
following structure:
model<-lme(response~Covariate+TreatmentA+TreatmentB,random=~1|Block/Plot,method="ML")
When I plot the residuals against the fitted values, I see a clear
positive trend (meaning that the variance increases with the mean).
I tried to solve this issue using weights=varPower(),
2005 Jan 27
4
self-written function
Dear all,
I?ve got a simple self-written function to calculate the mean + s.e.
from arcsine-transformed data:
backsin<-function(x,y,...){
backtransf<-list()
backtransf$back<-((sin(x[x!="NA"]))^2)*100
backtransf$mback<-tapply(backtransf$back,y[x!="NA"],mean)
2009 Aug 20
1
nested, repeated measure lme
...ndom=~time|unique.ID.of.every.individual
I have read that (2) is the only approach that works. But how could I then still include the nesting information from (1)?
Many thanks for your help!
Best wishes
Christoph
(using R 2.9.0 and the nlme library on Windows XP)
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2008 Jul 24
2
What is wrong with this contrast matrix?
...on-zero rows in the matrix, P.
Is that correct? And how can I assess if the matrix is orthogonal (given that it is non-symmetrical,
hence det(P) and other matrix operations won?t work)
Many thanks for your help!
Best wishes
Christoph.
(using R 2.7.1 on Windows XP)
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2008 May 27
5
Rotated text on a regression line
Dear all,
I stumbled over a problem recently when trying to use srt with text() on a
windows device.
What I intended to do was to plot a simple regression line, and to rotate
a piece of text such that the text has the same angle as the regression
line.
However, the text is always plotted in a slightly wrong angle:
####
x=1:10 #create arbitrary x and y values
y=x*2-rnorm(1:10)
2008 Jun 07
1
Multivariate LM: calculating F-values after calling linear.hypothesis
...79)
##
The last two lines of code are what I am stuck with; I thought using
diag() on the SSPH and SSPE matrices should give me the sums of squares
for hypothesis and error; so division should yield the F values?
I would be most grateful for any help!
Best wishes
Christoph
--
Dr. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen, Germany
www.gwdg.de/~cscherb1
2007 May 07
4
Mardia's multivariate normality test
Dear all,
I got this error message
> library(dprep)
> mardia(Savg)
Error in cov(data) : 'x' is empty
But with the same data, I got
> library(mvnormtest)
> mshapiro.test(Savg)
Shapiro-Wilk normality test
data: Z
W = 0.9411, p-value = 0.6739
What does the error message "Error in cov(data) : 'x' is empty" mean? Thanks a lot!
Jiao
2005 May 04
2
JGR and R 2.0.X
Dear all,
I am running Windows XP with several parallel installations of R (2.0.1;
2.1 and so on). How can I install JGR for the 2.0.1 version? I keep on
getting error messages when trying to install it.
Best wishes
Christoph
2005 Nov 09
1
trellis: positioning of key
Dear R users,
Using xyplot(), how can I position the key in the *margin* outside the
plotting area ?
My problem is that the key always overlaps with the x axis labels, no
matter how I try to specify any of the par() arguments (e.g. oma()).
Many thanks for any suggestions!
Christoph
###
for information, here??s the code I use
par(oma=c(0,0,3,0)) ###this, I think, is what should be changed
2008 Oct 02
1
An AIC model selection question
...red, but it "wastes" 5 df compared to the other models.
Would it be allowed to select model3 instead, simply because it uses up less df and the delta-AIC
between model2 and model3 is just 1?
Many thanks for any suggestions/comments!
Best wishes
Christoph
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2009 Jan 07
1
Extracting degrees of freedom from a gnls object
Dear all,
How can I extract the total and residual d.f. from a gnls object?
I have tried str(summary(gnls.model)) and str(gnls.model) as well as gnls(), but couldn?t find the
entry in the resulting lists.
Many thanks!
Best wishes
Christoph
--
Dr. rer.nat. Christoph Scherber
University of Goettingen
DNPW, Agroecology
Waldweg 26
D-37073 Goettingen
Germany
phone +49 (0)551 39 8807
fax +49 (0)551 39 8806
Homepage http://www.gwdg.de/~cscherb1
2012 Jul 20
1
Creating a pdf with layers?
Dear all,
Is it possible to create a pdf file with layers using the pdf() device in R?
Many thanks for your help!
Christoph
(using R 2.15.1 on Windows 7 64-Bit)
--
PD Dr Christoph Scherber
Georg-August University Goettingen
Department of Crop Science
Agroecology
Grisebachstrasse 6
D-37077 Goettingen
Germany
phone 0049 (0)551 39 8807
fax 0049 (0)551 39 8806
http://www.gwdg.de/~cscherb1