Displaying 20 results from an estimated 10000 matches similar to: "error in color-labeling in plot.lm"
2004 Feb 24
1
rstandard does not produce standardized residuals
Dear all,
the application of the function rstandard() in the base package
to a glm object does not produce residuals standardized to
have variance one:
the reason is that the deviance residuals are divided
by the dispersion estimate and not by the
square root of the estimate for the dispersion.
Should the function not be changed to produce residuals
with a variance about 1?
R 1.8.1 on
2003 Nov 17
1
confint: which method attached?
the function
confint
uses the profiling method of the function of the package MASS
confint.glm
even after the package has been detached!
1: might this be the intenden behavior?
2. How does the function remember its 'MASS' functionality after detaching the package?
R: 1.8.0; Windows 2000
Here is a sample program
> set.seed(7882)
> x<-rep(c(0,1),c(20,20))
>
2006 May 03
1
Problem in using confint method on polr model object
I fit a proportional odds model
with the polr-function of the MASS package from
Venables and Ripley
Applying the confint method to calculate confidence intervals for the
parameters I get
the following error message
Waiting for profiling to be done...
Re-fitting to get Hessian
Error in X[, -i, drop = FALSE] : incorrect number of dimensions
Can someone explain the error-message?
(The
2011 Mar 04
2
glht: Problem with symbolic contrast for factors with number-levels
Using a factor with 'number' levels the straightforward
symbolic formulation of a contrast in 'glht' of
the 'multcomp' package fails.
How can this problem be resolved without having to redefine the factor levels?
Example:
#A is a factor with 'number' levels
#B similar factor with 'letter' levels
dat<-data.frame(y=1:4,A=factor(c(1,1,2,2)),
2011 Sep 05
0
glht (multcomp): NA's for confidence intervals using univariate_calpha (fwd)
fixed @ R-forge. New version should appear on CRAN soon.
Thanks for the report!
Torsten
>
> ---------- Forwarded message ----------
> Date: Sat, 3 Sep 2011 23:56:35 +0200
> From: Ulrich Halekoh <Ulrich.Halekoh at agrsci.dk>
> To: "r-help at r-project.org" <r-help at r-project.org>
> Subject: [R] glht (multcomp): NA's for confidence intervals using
2004 May 04
2
Seeing the definition of a function
Dear all,
I was trying to see how the function 'confint' is defined. Doing
> confint
function (object, parm, level = 0.95, ...)
UseMethod("confint")
<environment: namespace:stats>
does not really enlighten me. How can I get to see the implementation (I guess it should be possible according to the general philosophy of the R project)?
Thanks in advance
S??ren
2011 Sep 03
0
glht (multcomp): NA's for confidence intervals using univariate_calpha
Hej,
Calculation of confidence intervals for means
based on a model fitted with lmer
using the package multcomp
- yields results for calpha=adjusted_calpha
- NA's for calpha=univariate_calpha
Example:
library(lme4)
library(multcomp)
### Generate data
set.seed(8)
d<-expand.grid(treat=1:2,block=1:3)
e<-rnorm(3)
names(e)<-1:3
d$y<-rnorm(nrow(d)) + e[d$block]
2002 Dec 09
2
R as a COM client - is it possible?
Dear all,
In S+, there are functions like
create.ole.object
call.ole.method
release.ole.object
for communicating with other programs which work as a COM server (on
Windows).
Is it possible to do something similar in R (I've studied the 'connections'
facilities, but they do not seem to work).
==========================================
S?ren H?jsgaard, PhD, Senior Scientist
2005 Jun 13
1
Error in load(zfile, envir = envir) : input has been corrupted, with LF replaced by CR
I am trying to build a package binary, and get the message below. Can anyone point me to a solution to that problem.
Thanks in advance
S?ren
....
installing data files
installing man source files
installing indices
Error in load(zfile, envir = envir) : input has been corrupted, with LF replaced
by CR
Execution halted
make[2]: *** [indices] Error 1
make[1]: *** [all] Error 2
2006 Jan 31
0
lattice: combining panel.xyplot with panel.abline - is thispossible?
Well, one way to do it is
xyplot(y~time|id, data=dat,type='l',
panel=function(x,y,subscripts,...){
panel.xyplot(x,y,subscripts,...)
panel.abline(v=dat[subscripts,"cm1"])
panel.abline(v=dat[subscripts,"cm2"])
}
)
Since I don't know what the dataframe 'mergeData' is I have used 'dat'.
Best regards
Frede Aakmann
2007 Jul 18
2
remove columns having a partial match name
Dear all,
I would like to know how can I retrieve a data.frame without the columns that have a partial match name. Let´s say that I have a data.frame with 200 columns and 100 of them have the name "StartX", with X being the unique part for each column name. I want to delete all columns that have the name starting with "Start". I´ve tried to do this but it doesn´t work:
>
2007 Jul 20
3
binned column in a data.frame
Dear all,
I would like to know how can I create a binned column in a data.frame. The output that I would like is something like this:
Start Binned_Start
1 0-5
2 0-5
6 5-10
8 5-10
13 10-15
...
Best regards
João Fadista
Ph.d. student
UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of Genetics and Biotechnology
Blichers
2011 Jun 06
1
Log file of building vignette in RCMD check
Hello,
I am trying to run "RCMD check" on a package. It performs OK, with the exception of a single warning:
* checking package vignettes in 'inst/doc' ... WARNING
Package vignette(s) without corresponding PDF:
AnnotationFuncsUserguide.Rnw
As I understand the problem, the vignette is not built. However, I have no problems in doing so manually.
How do I fix this problem? I
2007 Mar 23
2
concatenate 2 data.frames
Dear all,
I would like to know how can I concatenate 2 data.frames into a single one. Both data frames have the same number of columns and the same class type in each correspondent column. So what I want is to have a new data.frame where I have first the values from one data.frame and then the values from a second data.frame would came after in this new data.frame.
Thanks in advance.
Med
2008 Apr 08
4
permutation test assumption?
Dear all,
Can I do a permutation test if the number of individuals in one group is much bigger than in the other group? I searched the literature but I didin´t find any assumption that refers to this subject for permutation tests.
Best regards
João Fadista
Ph.d. student
UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of Genetics and Biotechnology
Blichers Allé 20, P.O.
2007 Mar 30
1
Model comparison
Dear all,
I would like to know if I can compare by a significance test 2 models with different kind of parameters. Perhaps I am wrong but I think that we can only compare 2 models if one is a sub model of the other.
Med venlig hilsen / Regards
João Fadista
Ph.d. studerende / Ph.d. student
AARHUS UNIVERSITET / UNIVERSITY OF AARHUS
Det Jordbrugsvidenskabelige Fakultet / Faculty of
2007 May 15
0
sliding window approach
Dear all,
I would like to know if there is any R package that uses a sliding window approach to assess statistical significance out of my data. My data is composed of DNA sequences (of variable length) that are mapped to a genome with a determined score of alignment.
So, I want to see if I can find more tags in a given region of the genome as opposed to finding them by chance. In this sliding
2007 May 31
0
distribution of peaks in random data results
Dear all,
I have the positions of N points spread through some sequence of length L (L>N), and I would like to know how can do the following:
1- Permute the positions of the N points along the whole sequence.
Assuming a uniform distribution I did: position1 <- runif(N, 1, L)
2- Apply a kernel convolution method to the resulting permuted points profile.
For this I applied the
2007 Jul 19
0
test about distribution of data in a single population
Dear all,
I would like to know how can I test which are the intervals of my data that have significant less or more counts than the other intervals.
Example:
Interval [1:200] [200:400] [400:600] [600:800] ... more 900 hundred columns
Count 12 28 7 5
Thanks in advance,
Best regards
João Fadista
Ph.d.
2007 Sep 05
1
Running geeglm unstructured corstr
? stato filtrato un testo allegato il cui set di caratteri non era
indicato...
Nome: non disponibile
Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070905/6d1002c1/attachment.pl