Hello,
I have difficulties combining boxes and lines in plot legend. I
searched previous R-posts and found this (with no solution):
http://tolstoy.newcastle.edu.au/R/help/06/07/30248.html. Is there a
way to avoid boxes behind the line legends?
x1 - rnorm(100)
x2 - rnorm(100, 2)
hist(x1, main = ,
Check out:
http://tolstoy.newcastle.edu.au/R/e2/help/07/05/16777.html
On 9/10/07, Lauri Nikkinen [EMAIL PROTECTED] wrote:
Hello,
I have difficulties combining boxes and lines in plot legend. I
searched previous R-posts and found this (with no solution):
Hello,
I would like to know what can I do if I use strplit with a string and I want
to use the middle left,I mean I have this:
strsplit(bA531F16-rep,\\-)
[[1]]
[1] bA531F16 rep
I would like to work just with bA531F16 in another variable, what could I
do?, Thank you
Why not use the more simple
xyplot(total.fat~x|variable,groups=Group,
data=tmp1,type=c(p,r))
???
See ?panel.xyplot and especially the type argument of that panel function.
Best regards
Frede Aakmann Tøgersen
Scientist
UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of
Thanks Gabor, I got it!
For example:
x1 - rnorm(100)
x2 - rnorm(100, 2)
hist(x1, main = , col = orange,ylab = density, xlab = x, freq
= F, density = 55, xlim = c(-2, 5), ylim = c(0, 0.5))
par(new = T)
hist(x2, main = , col = green, ylab = , xlab = ,axes = F, xlim
= c(-2, 5), ylim = c(0, 0.5),
Dear Friends.
I found something very puzzling with constOptim(). When I change the
parameters for ConstrOptim, the error messages do not seem to be
consistent with each other:
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
Error in constrOptim(c(0.5, 0.3, 0.5), f =
Dear Rosa,
please be more specific. Statistical tests for which hypothesis?
For example, some tests can be made robust using Heteroskedasticity-
*and Autocorrelation-* Consistent (HAC) covariance matrices in package
'sandwich': see
- waldtest{lmtest} for a redundant variables test much like
unlist(strsplit(bA531F16-rep,\\-))[1]
[1] bA531F16
--- Carlos Morales [EMAIL PROTECTED]
wrote:
Hello,
I would like to know what can I do if I use
strplit with a string and I want to use the middle
left,I mean I have this:
strsplit(bA531F16-rep,\\-)
[[1]]
[1] bA531F16 rep
Yogesh Tiwari wrote:
Hi,
How to plot a variance over a data point, something like error bar.
Hi Yogesh,
The generic method is to use the arrows function, and there are quite
a few variations on this theme (in alpha order):
brkdn.plot(plotrix)
dispbars(plotrix)
errbar(Hmisc and sfsmisc)
Thanks Frede
I didn't know about the r type.
Ross Darnell
-Original Message-
From: Frede Aakmann Tøgersen [mailto:[EMAIL PROTECTED]
Sent: Mon 10-Sep-07 4:45 PM
To: Ross Darnell; r-help@stat.math.ethz.ch
Subject: SV: [R] lattice panel.lmline problem
Why not use the more simple
prcomp() in stats handles matrices with n p well, IMO.
--
Bjørn-Helge Mevik
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide
PS == Paul Smith [EMAIL PROTECTED]
on Sun, 9 Sep 2007 12:17:32 +0100 writes:
PS On 9/9/07, kevinchang [EMAIL PROTECTED] wrote:
I tried to made the matrix with this size by either matrix() or array().
However, there seems to be default limit of number for rows made. I got
sort
Hi,
I have been suggested to use Information theoretic approach instead of a
Geneal linear model for a multiple regression. Is is possible to perform
this type of anlaysis in R?
Thanks,
Adela
--
Adela González Megías
Depto. Biología Animal
Fac. Ciencias
Universidad de Granada
18071
How do I install and load the BradleyTerry add on package in R 2.5.1 in
MSWindowsXP environment?
Kalyan Roy
Indian Market Research Bureau (IMRB) International
New Delhi, India
[[alternative HTML version deleted]]
__
However you may also need to install the package brlr, since the
BradleyTerry package depends on this.
For Windows users, it's usually easiest to install packages using the
Packages menu in the RGui - any dependencies are then automatically
installed.
Heather
-Original Message-
From:
On 9/10/07, Martin Maechler [EMAIL PROTECTED] wrote:
PS On 9/9/07, kevinchang [EMAIL PROTECTED] wrote:
I tried to made the matrix with this size by either matrix() or
array().
However, there seems to be default limit of number for rows made. I
got sort
of error message
Kalyan Roy (DEL/MSG) wrote:
How do I install and load the BradleyTerry add on package in R 2.5.1 in
MSWindowsXP environment?
Hi Kalyan,
If
R CMD INSTALL
doesn't work, you can use WinZip or Zip Reader to unzip the package to:
C:\Program Files\R-2.5.1\library
or whatever your path to the
Hello,
Some weeks ago, thanks to you, I managed to install R, to connect to a local
MySQL Database and to launch some queries with a script written with Tinn-R.
My script is now ok and would like to test it with the real database.
I did the same installation of R, DBI package and RMySQL package
Hi,
I'm looking for an implementation of the corrected resampled t-test
[1] to compare to different machine learning classifiers. Any idea,
whether this is implemented in a R package?
Thanks in advance,
Christian
[1] Nadeau Bengio, 2003, Inference for the Generalization Error,
Machine
This is not quite what you want but you can try this:
legend(3, 0.45, legend = c(x1, x2, mean(x1), mean(x2)), col =
c(orange, green),pch = c(15,15,-1,-1), lty=c(-1,-1,2,2))
Although pch=22 should draw a filled square with a border - but it draws only
the border instead
Monica
I'm trying to build an extension that uses the \u characters, but
building the library fails. The Writing R Extensions guide says
that these characters should not be used, but I really can't find any
way around them. Is there any way to configure the building of a
library so that it will
Hi,
sapply(formals(readBin), mode)
con what n sizesignedendian
namename numeric logical logicalcall
returns for the mode of size logical. But in the documentation is said
that size should be integer. Does anyone know why the mode is logical?
Thanks
The following is one of the examples in the help page for histogram:
histogram( ~ height | voice.part, data = singer,
xlab = Height (inches), type = density,
panel = function(x, ...) {
panel.histogram(x, ...)
On 9/10/2007 10:26 AM, Sigbert Klinke wrote:
Hi,
sapply(formals(readBin), mode)
con what n sizesignedendian
namename numeric logical logicalcall
returns for the mode of size logical. But in the documentation is said
that size should be
You could directly use
strsplit(bA531F16-rep,\\-)[[1]][1]
and
strsplit(bA531F16-rep,\\-)[[1]][2]
Regards
Carlos Morales-2 wrote:
Hello,
I would like to know what can I do if I use strplit with a string and I
want to use the middle left,I mean I have this:
I have recently been introduced to the ggplot package by Hadley Wickham
and must say I am quite impressed so far at how easy it is to make
attractive plots, but one thing I am struggling over is how to
consolidate legends.
It's not currently possible to consolidate them (although in the
Hi,
I was trying to install the packageRWinEdt in my computer with Vista OS.
Once I finished the installation, R just cannot load this library at all,
reporting some error information. I reinstalled R and Winedit and reloaded
this package. RWinEdt seemed to register in the computer but never work
Hi, there:
just a little bit off-topic:
any algorithm in classification is good for ordered features, like all
variables used are 1, 2, 3,...; not really continuous but ordered. i
tried random forest and dlda already.
thanks.
--
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.
Did you
Hi,
I was trying to install the packageRWinEdt in my computer with Vista OS.
O
Use Tinn-R it works with Vista although has some minor issues (which are
probably system specific).
https://sourceforge.net/projects/tinn-r
Stefan
-=-=-
... Time is an illusion, lunchtime doubly so. (Ford
Yuchen Luo wrote:
Dear Friends.
I found something very puzzling with constOptim(). When I change the
parameters for ConstrOptim, the error messages do not seem to be
consistent with each other:
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
Error in
On Thu, 6 Sep 2007, eugen pircalabelu wrote:
Good afternoon!
I'm trying to use the Survey package for a stratified sample which has
4 criteria on which the stratification is based. I would like to get the
corrected weights and for every element i get a weight of 1
E.g: tipping
On Thu, 6 Sep 2007, Tobias Verbeke wrote:
eugen pircalabelu wrote:
I'm trying to use the survey package to get a better point of view
for my data, but i need some piece of advice:
i have some data from a survey which has been stratified using 2
criteria: region(7 values), size of
Steve,
This example works:
x-TukeyHSD(fm1, tension, ordered = TRUE)
as.table(x$tension)
Steve Powers wrote:
So I've come across a few cases where complex outputs from functions
will not write to tables. The most recent case involves the TukeyHSD
function in the stats package. If I save
The error message about the feasible region comes from constrOptim(),
before your function is called. The error message about missing lambda1
comes from calling your function.
-thomas
On Sun, 9 Sep 2007, Yuchen Luo wrote:
Dear Friends.
I found something very puzzling with
On Sun, 9 Sep 2007, eugen pircalabelu wrote:
A short example:
stratum id weight nh Nh y sex
1 1 3 5 15 23 1
1 2 3 5 15 25 1
1 3 3 5 15 27 2
1 4 3 5 15 21 2
1 5 3 5 15 22 1
2 6
This example works fine:
test-matrix(c(1,2,'VOICIUNPETITTES',3),ncol=2,nrow=2)
write.csv(test,file='C:/xavier/test.csv')
Could you provide the same small example when it doesn't work?
kwaj wrote:
Hello,
I have a peculiar problem which I am hoping I can get help
On Mon, 10 Sep 2007, Sigbert Klinke wrote:
Hi,
sapply(formals(readBin), mode)
con what n sizesignedendian
namename numeric logical logicalcall
returns for the mode of size logical. But in the documentation is said
that size should be integer. Does
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is what is the best OS on PC (laptop) for statistical
computing and why.
I really appreciate your insight.
Have a nice day.
__
Hi Hadley,
I just tried out your suggestion, but it does not look like the
get_legends function is working correctly. Instead of returning a grob
back to me it returns NULL.
Here is my modified code and the results of running it.
Any help would be appreciated. I believe that once I can get the
Dear R users,
Below I have written 4 functions CIT1, CIT2a and CIT2b and CIT3 which
recode a variable CLD_ISCH into 3 new variables(T1 T2 T3), I wish to use
T1, T2 and T3 based on the values of tf1 and tf2.
(NOTE:- T2a is used to create T2 in a long winded manner due to my lack
of programming
Hello,
I am new to R exploratory data analysis and plotting. Is anyone aware of a way
to overlay a set of conditional histograms with conditional PDFs? Below, I
generate a lattice plot of precipitation histograms based on different months
and stations, given a subset of the dataset:
You want whatever all the people you are working with are using
to make it as easy as possible to work together with them.
On 9/10/07, Wensui Liu [EMAIL PROTECTED] wrote:
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is
Are there any packages in R that reproduce the package resample of S-Plus?
The sample() function in R doesn't provide equivalent flexibility of
bootstrap() and bootstrap2().
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL
I have 316 files storing a number of cycles of the respiratory signal
from 316 different patients. Each file record is made up of the
following data:
amplitude,phase,timestamp,validflag,ttlin,mark,ttlout
amplitude: is the position in cm relative to an arbitrary reference
(signal amplitude)
useRs,
I am looking to find the minimum positive value of some data I have.
Currently, I am able to find the minimum of data after I apply some other
functions to it:
x
[1] 1 0 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10 10
sort(x)
[1] 0 1 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10
Try this:
min(diff(sort(x))[diff(sort(x))0])
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O
On 10/09/2007, dxc13 [EMAIL PROTECTED] wrote:
useRs,
I am looking to find the minimum positive value of some data I have.
Currently, I am able to find the minimum of
On Mon, 2007-09-10 at 11:20 -0700, dxc13 wrote:
useRs,
I am looking to find the minimum positive value of some data I have.
Currently, I am able to find the minimum of data after I apply some other
functions to it:
x
[1] 1 0 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10 10
Hello friends,
I loaded R 2.4.1 onto a Fedora Core 6 Linux box (taking all defaults). Then
I ran these commands from within R:
options(CRAN=http://cran.stat.ucla.edu;)
install.packages(CRAN.packages()[,1])
As a new user of R, I was shocked when I finished loading R and discovered
the following
Sorry, I should have mentioned that get_legend won't work on the plots
that you are actually plotting - you have turned their legends off!
You'll need a plot which isn't plotted, but is used to produce the
legends.
Hadley
On 9/10/07, Te, Kaom [EMAIL PROTECTED] wrote:
Hi Hadley,
I just tried
http://www.burns-stat.com/pages/Tutor/bootstrap_resampling.html
includes a synopsis of R packages that do bootstrapping.
It is brief and incomplete, but hopefully useful.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling
I am having some trouble getting the persp() package to change the x
and y axis on a 3d plot. It defaults to the [0,1] interval and when I
try to change it I get errors.
Example:
This works:
D - c(1,2,3,4,5,6,7,8,9,10)
M - c(11,12,13,14,15,16,17,18,19,20)
DM - cbind(D,M)
persp(DM,
A Lenzo wrote:
Hello friends,
I loaded R 2.4.1 onto a Fedora Core 6 Linux box (taking all defaults). Then
I ran these commands from within R:
options(CRAN=http://cran.stat.ucla.edu;)
install.packages(CRAN.packages()[,1])
As a new user of R, I was shocked when I finished loading R and
On 9/10/07, Ross Darnell [EMAIL PROTECTED] wrote:
Thanks Frede
I didn't know about the r type.
For the record, this is probably what you wanted:
xyplot(...
panel = panel.superpose,
panel.groups = function(x, y, ...) {
panel.xyplot(x, y, ...)
if (length(x)
I am just trying to teach myself how to use the mle function in R because it is
much better than what is provided in MATLAB. I am following tutorial material
from the internet, however, it gives the following errors, does anybody know
what is happening to cause such errors, or does anybody know
On 11/09/2007, at 4:22 AM, Wensui Liu wrote:
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is what is the best OS on PC (laptop) for statistical
computing and why.
I really appreciate your insight.
Have a nice day.
Linux! Mac OS is ok to me2.
On 9/10/07, Wensui Liu [EMAIL PROTECTED] wrote:
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is what is the best OS on PC (laptop) for statistical
computing and why.
I really appreciate
Dear Professor Murdoch.
Thank you for your help!
1. I believe c(0.5,0.3,0.5) satisfies the constrain because I did the
following experiment
ui=-1*ui
ci=-1*ci
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
The same error message pops up. Any theta ( in this case,
Terence Broderick wrote:
I am just trying to teach myself how to use the mle function in R because it
is much better than what is provided in MATLAB. I am following tutorial
material from the internet, however, it gives the following errors, does
anybody know what is happening to cause such
Hi
I was wondering if someone could help me with an rpart problem. I can see
that cross-validation is the default for tree selection in rpart -- has a
bootstrap method been implemented anywhere? I think this is a different
thing to 'bagging' or 'boosting' -- I still want 'one' tree at the end,
From ?trackObjs:
Overview of trackObjs package
Description:
The trackObjs package sets up a link between R objects in memory
and files on disk so that objects are automatically resaved to
files when they are changed. R objects in files are read in on
demand and do not
Hello All,
I have searched many help forums, message boards, etc. and I just can't
apply the comments to what I need my program to do. I am running R 2.5.1 on
an XP system, and my desire is to produce replicate datasets for a
simulation study I am running. Essentially, I have sets of
See the article Cluster in R Task Views
http://cran.at.r-project.org/src/contrib/Views/Cluster.html
It lists names of packages for clustering analysis, which you can install.
Or, go to choose another CRAN mirror closest to you from the R web site,
then click on 'Task views' in the left frame.
My sense is that R users are even split between UNIX and Windows
users so either will do in terms of the larger community.
Some R packages may not be avaliable on every platform or will
be available on one platform before another or there will be
certain platform-specific issues. So in the end
Hi Jonathan,
What exactly do you mean by replication?
Do you want to keep a1,b1,c1,... unchanged but have 30
different sets of random numbers?
Regards,
Moshe.
--- VTLT1999 [EMAIL PROTECTED] wrote:
Hello All,
I have searched many help forums, message boards,
etc. and I just can't
Here are some solutions each of which
1. has only one line,
2. x only occurs once so you can just plug in a complex expression
3. no temporary variables are left
min(sapply(x, function(z) if (z 0) z else Inf))
(function(z) min(ifelse(z 0, z, Inf))) (x)
with(list(z = x), min(z[z 0]))
local({
x-seq(-1,1,length=10)
y-seq(-1,1,length=10)
a-matrix(c(1,2,2,1),2,2)
b-matrix(c(2,1,1,2),2,2)
fv-function(x,y) {
m-x*a+y*b
t-m[1,1]+m[2,2]; d-m[1,1]*m[2,2]-m[1,2]^2
return((t-sqrt(t^2-4*d))/2)
}
gv-function(x,y) {
t-x*(a[1,1]+a[2,2])+y*(b[1,1]+b[2,2])
Its a FAQ:
http://hermes.sdu.dk/Rdoc/faq.html#Why%20does%20outer()%20behave%20strangely%20with%20my%20function%3f
On 9/10/07, Jan de Leeuw [EMAIL PROTECTED] wrote:
x-seq(-1,1,length=10)
y-seq(-1,1,length=10)
a-matrix(c(1,2,2,1),2,2)
b-matrix(c(2,1,1,2),2,2)
fv-function(x,y) {
Dear Listers,
I am a little tired of installing all packages I want every time when
I instill a new version of R.
Say, if I have a list of packages I need to use, is it possible to
tell R to install them all for me automatically rather than I install
them one by one?
Thx.
Try
?update.packages
Ross Darnell
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Wensui Liu
Sent: Tuesday, 11 September 2007 11:45 AM
To: r-help@stat.math.ethz.ch
Subject: [R] install packages automatically
Dear Listers,
I am a little tired of
Version 1.1-8 of package 'scuba' has been uploaded to CRAN.
'scuba' is a package for scuba diving calculations and decompression
models. It supports dive profiles (tables, plotting etc), analysis of dive
profiles using decompression models, gas toxicity calculations, and gas
usage calculations.
Mange tak!
FYI, this is the way it is able to run (I was going to attach
station.precip.R, but I read that attaching files is not recommended
- let me
know if you would like it)
x - dget(file=C://Documents and Settings/Bradley/My
Documents/Arizona/CourseResources/ATMO529/station.precip.R)
I am using 'xyplot' in lattice to plot some data where the x-axis is a
POSIXct date. I have data which spans a 6 month period, but when I
plot it, only the last month is printed on the right hand side of the
axis. I would have expected that at least I would have a beginning
and an ending point
Either
min(diff(sort(x))[diff(sort(x))0])
or
min(diff(sort(unique(x
--- dxc13 [EMAIL PROTECTED] wrote:
useRs,
I am looking to find the minimum positive value of
some data I have.
Currently, I am able to find the minimum of data
after I apply some other
functions to it:
x
Fiona Callaghan wrote:
I was wondering if someone could help me with an rpart problem. I can see
that cross-validation is the default for tree selection in rpart -- has a
bootstrap method been implemented anywhere? I think this is a different
thing to 'bagging' or 'boosting' -- I still want
Dear friends.
I use ConstrOptim( ) and got error message initial value not available.
My understanding of initial value not available is that one of the
following 3 cases happens:
1.The objective function is not well defined at the point of the initial
value.
2. The differentiation of the
Hi guys,
I tried to made the matrix with this size by either matrix() or array().
However, there seems to be default limit of number for rows made. I got sort
of error message from R .To be specific,
m--matrix(ncol=3,nrow=10)
error message:[ reached getOption(max.print) -- omitted 7
On 9/9/07, kevinchang [EMAIL PROTECTED] wrote:
I tried to made the matrix with this size by either matrix() or array().
However, there seems to be default limit of number for rows made. I got sort
of error message from R .To be specific,
m--matrix(ncol=3,nrow=10)
error message:[ reached
Hello,
using Sweave, is there any option to preserve the
original aspect ratio of plots generated from R code?
Consider this Sweave chunk:
test,echo=F,fig=T,width=2,height=2=
x - 1:10
y - sin(x)
par(mar=c(4,4,0,4))
plot(x,y,
xlab=x label,
ylab=y label
)
@
In Latex, I want to
Your problem is not related to missing file permissions but to a problem
with the linker (ld). As you can see from the error message, ld cannot
find a compatible libR.so so compilation fails. You did not specify
which architecture you're on but in case it is a 64-bit platform, the
following thread
On 09/09/2007 7:01 AM, Yuchen Luo wrote:
Dear friends.
I use ConstrOptim( ) and got error message initial value not available.
My understanding of initial value not available is that one of the
following 3 cases happens:
1.The objective function is not well defined at the point of the
On 09/09/2007 7:24 AM, Werner Wernersen wrote:
Hello,
using Sweave, is there any option to preserve the
original aspect ratio of plots generated from R code?
Consider this Sweave chunk:
test,echo=F,fig=T,width=2,height=2=
x - 1:10
y - sin(x)
par(mar=c(4,4,0,4))
plot(x,y,
xlab=x
Have a look at rbind.fill in the reshape package.
Hadley
On 9/9/07, Muenchen, Robert A (Bob) [EMAIL PROTECTED] wrote:
Hi All,
If I need to stack two data frames, I can use rbind, but it requires
that all variables exist in both sets. I can make that happen, but other
stat packages would
Perfect. Thanks Hadley!
-Original Message-
From: hadley wickham [mailto:[EMAIL PROTECTED]
Sent: Sunday, September 09, 2007 10:11 AM
To: Muenchen, Robert A (Bob)
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] stacking data frames with different variables
Have a look at rbind.fill
Wolfgang,
This looks great and certainly puts to shame my code .
Thanks,
Monica
Date: Fri, 7 Sep 2007 21:04:32 +0100 From: [EMAIL PROTECTED] To: [EMAIL
PROTECTED] CC: r-help@stat.math.ethz.ch Subject: Re: [R] confusion matrix -
better code? Dear Monica, try this: cm =
Hi,
This may be a stupid question, but why when creating an AR model, why does R
set the first n components of the ar object's resid vector to NA (where n is
the order of the model)?
Regards,
NC
[[alternative HTML version deleted]]
__
This is bizzare. The for loop I set seems to go infinitly because the
hourglass cursor never disappers and I never get the result. But when I
click on X of the window try to exit R, the execution finishes right after
I abort exiting when RGui asks me to save the workplace or not. Please help
On 09/09/2007 1:36 PM, kevinchang wrote:
This is bizzare. The for loop I set seems to go infinitly because the
hourglass cursor never disappers and I never get the result. But when I
click on X of the window try to exit R, the execution finishes right after
I abort exiting when RGui asks me
I am trying to fit the chi-squared distribution to a set of data using the
fitdistr function found in the MASS4 library, the data set is called ONES3, I
have loaded it using the command
ONES3-read.table(ONES3.pdf,header=TRUE,na=NA)
I print out the dataset ONES3 to the screen to make
I assume that you want to do the fitdistr on one of the columns of the
dataframe that you have read in. What does 'str(ONES3)' show? If the
data is in the first column, try:
fitdistr(ONES3[[1]],chi-squared)
On 9/9/07, Terence Broderick [EMAIL PROTECTED] wrote:
I am trying to fit the
Hi R-users!
I have a problem with the survey package and i would be very grateful if you
can help me.
A short example:
stratum id weight nh Nh y sex
1 1 3 5 15 23 1
1 2 3 5 15 25 1
1 3 3 5 15 27 2
1 4 3 5
Dear Professor Mordoch.
Thank you very much for your help! Your time is highly appreciated!
I do intend to optimize over 3 parameters and the way I did it is
constrOptim(c(0.5,0.3,0.5), fit.error, fit.error.grr, ui=ui,ci=ci).
Also, There is a missing value for both sigmae and ss. (The last 3rd
Many thanks for the quick reply, Duncan.
Now I see that something is going wrong in R already.
I overlooked this because the R window was so big. I
have to persuade plotViewport() not to create a square
area but one which takes the different margins into
account.
Best regards,
Werner
---
QCA implements the Qualitative Comparative Analysis using a boolean
minimization algorithm for data coded with presence/absence of the causal
conditions that affects a phenomenon of interest.
This new release has an experimental function that obtains the same exact
solutions as the main
Thanks all!
On 09/09/2007, Nick Chorley [EMAIL PROTECTED] wrote:
Hi,
This may be a stupid question, but why when creating an AR model, why does
R set the first n components of the ar object's resid vector to NA (where n
is the order of the model)?
Regards,
NC
[[alternative
I am wanting to generate panels showing scatterplots with the linear
fitted line for two groups within each panel superimposed.
I have two conditioning factors, variable and Group and I want
separate panels for each level of variable
with different symbols and lmlines for each level of Group.
Hi Gabor,
Wow, this is awesome although I eventually should learn MySQL for
integrating it on web-based DB management using PHP or Perl, this is a
very helpful tool for me to start with!
Thank you very much
Gabor Grothendieck wrote:
Others have already pointed out %in% but regarding
Great. Regarding the web, note that there are actually quite a few R
web projects as well:
http://www.lmbe.seu.edu.cn/CRAN/doc/FAQ/R-FAQ.html#R-Web-Interfaces
I have used rpad (www.rpad.org) which has an integrated web server right
in the R package making setup a non-issue.
On 9/8/07,
François Pinard wrote:
[Roland Rau]
[François Pinard]
I wonder what happened, for R to hide the underlying Scheme so fully,
at least at the level of the surface language (despite there are
hints).
To further foster portability, we chose to write R in ANSI
Hi. I need to use a few different clustering functions. I managed to run the
kmeans() one which is in my stats library, but I can't use any function,
such as agnes(), that is in my cluster library. Any idea how to access
other libraries?
Thanks!
--
View this message in context:
I am looking at some extremes, and I want to implement covariates in my
analysis. My work is on daily datasets, and I use the cluster of exceedences
approach in fpot (evd-package) to estimate the parameters in the General
Pareto Distribution. How to include covariates like e.g. time is my
1 - 100 of 88815 matches
Mail list logo