the data is :
table.8.3-data.frame(expand.grid(
marijuana=factor(c(Yes,No),levels=c(No,Yes)),
cigarette=factor(c(Yes,No),levels=c(No,Yes)),
alcohol=factor(c(Yes,No),levels=c(No,Yes))),
count=c(911,538,44,456,3,43,2,279))
fit3-glm(count~.^3,poisson,table.8.3)
sumary(fit3)
...
Residual
I'm trying to fit a Gompertz sigmoid as follows:
x - c(15, 16, 17, 18, 19) # arbitrary example data here;
y - c(0.1, 1.8, 2.2, 2.6, 2.9) # actual data is similar
gm - nls(y ~ a+b*exp(-exp(-c*(x-d))), start=c(a=?, b=?, c=?, d=?))
I have been unable to properly set the starting value '?'s. All
Dear R-users,
is there anyone else with problems to get R-2.1.0 compiled on
FreeBSD6-CURRENT?
After typing '.configure' and then 'make' I get the following output:
-
[...snip...]
gcc -export-dynamic -L/usr/local/lib -o R.bin Rmain.o CConverters.o
Hello list,
I'd like to know if there is a function that tells in which package is a
given function.
Something like :
which.package(lda)
[1] MASS
Thank you.
Romain
--
~
~~ Romain FRANCOIS - http://addictedtor.free.fr
Ronggui asks:
: -Original Message-
: From: [EMAIL PROTECTED]
: [mailto:[EMAIL PROTECTED] On Behalf Of ronggui
: Sent: Saturday, 23 April 2005 5:41 PM
: To: r help
: Subject: [R] question about about the drop1
:
:
: the data is :
: table.8.3-data.frame(expand.grid(
:
Dear Lister,
I am working with Windows XP and R 2.1.0 and can check and build
home-made packages easily (just *.r, *.rmd, *.rda files, no compiled
code). However for some reasons, I cannot check or build the package
'foreign' from the source (I took it as an exercise...). After some
lines of
(Ted Harding) wrote:
On 21-Apr-05 [EMAIL PROTECTED] wrote:
Here's a bit of a refinement on Ted's first suggestion.
[ corrected from runif(M*k), N, k) to runif(N*k), N, k) ]
N - 1
graphics.off()
par(mfrow = c(1,2), pty = s)
for(k in 1:20) {
m - (rowMeans(matrix(runif(N*k), N, k)) -
On 4/23/05, Matthias Kohl [EMAIL PROTECTED] wrote:
here is another idea to illustrate the Central Limit Theorem which is
based on our package distr.
Thanks to all for your numerous suggestions. Since I am just beginning
with R, it will take a while before I can fully digest your help.
Have a
On 23-Apr-05 [EMAIL PROTECTED] wrote:
: -Original Message-
: From: [EMAIL PROTECTED]
: [mailto:[EMAIL PROTECTED] On Behalf Of
: michael watson (IAH-C)
: [...]
: I have a highly significant interaction term. In the context
: of the experiment, this makes sense. I can visualise the data
:
James Salsman wrote:
I'm trying to fit a Gompertz sigmoid as follows:
x - c(15, 16, 17, 18, 19) # arbitrary example data here;
y - c(0.1, 1.8, 2.2, 2.6, 2.9) # actual data is similar
It is a good practice to plot the data before trying to fit complicated
nonlinear models. If you are going
Kittler, Richard wrote:
I am using 'read.csv' in V2.0.1 to read in a CSV file with the
colClasses option and am getting an error from 'scan' when it encounters
a non-numeric value for a 'numeric' column, i.e.
ds - read.csv(in_file, nrows=irow, row.names=NULL,
colClasses=zclass,
Hui Han wrote:
Hi,
I used the density function in the R package, and got the following
results. I just wonder how to explain them.
What is Min, 1st Qu, Median, and so on? I could not find an explanation
from help(density). The plot doesn't seem to match
the x and y value either.
Thanks in
Romain Francois wrote:
Hello list,
I'd like to know if there is a function that tells in which package is a
given function.
Something like :
which.package(lda)
[1] MASS
Thank you.
Romain
getAnywhere(lda)
A single object matching 'lda' was found
It was found in the following places
Romain Francois wrote:
Hello list,
I'd like to know if there is a function that tells in which package is a
given function.
Something like :
which.package(lda)
[1] MASS
Thank you.
Romain
Perhaps easiest is
help.search(lda)
If you know that the name is recognized in your current session (which
Jessica Higgs wrote:
Hi all--
A quick question about the dr () function. I am using this function to
reduce the dimensions of a data set I have that involves 14 predictor
variables and one predictant or response. The goal is to discover which
variables play the most important role in
Thanks for your work on this topic, Marc, but the bug has already been
fixed by Brian 4 days ago (see below).
Uwe
r33976 | ripley | 2005-04-19 05:38:51 -0400 (Tue, 19 Apr 2005) | 1 line
Changed paths:
M /trunk/NEWS
M /trunk/src/library/utils/R/menu.R
missing braces
Dear All
Is there some way of restarting R, without quitting R? In other words,
I am looking for something like the commands reset (MuPAD) or restart
(Maple).
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
Patrick Giraudoux wrote:
Dear Lister,
I am working with Windows XP and R 2.1.0 and can check and build
home-made packages easily (just *.r, *.rmd, *.rda files, no compiled
code). However for some reasons, I cannot check or build the package
'foreign' from the source (I took it as an
Paul Smith wrote:
Dear All
Is there some way of restarting R, without quitting R? In other words,
I am looking for something like the commands reset (MuPAD) or restart
(Maple).
No, but of course you could write a function that fires up a new
instance of R and quits the current one after that.
A few moments ago I used Mozilla Thunderbird to send an email message
that contained a forwarded message sent from an R-project.org address.
The spell checker flagged this phrase as a potential spelling error and
offered the replacement R-project.Borg. If it wasn't that Mozilla is
another
Hi,
I'm reading some data from a file. The code below plots the all records in a
barplot. How can I plot only records 5 to 10?
chldata - read.table(test3.txt,header=TRUE)
barplot(rev(chldata$chlorophyll),horiz=TRUE,names.arg=rev(chldata$depth))
Thanks in advance,
Pieter
[[alternative
On Sat, 2005-04-23 at 16:14 +0200, Uwe Ligges wrote:
Thanks for your work on this topic, Marc, but the bug has already been
fixed by Brian 4 days ago (see below).
Uwe
r33976 | ripley | 2005-04-19 05:38:51 -0400 (Tue, 19 Apr 2005) | 1 line
Changed paths:
M /trunk/NEWS
On Sat, 2005-04-23 at 09:49 -0500, Douglas Bates wrote:
A few moments ago I used Mozilla Thunderbird to send an email message
that contained a forwarded message sent from an R-project.org address.
The spell checker flagged this phrase as a potential spelling error and
offered the
Pieter Provoost wrote:
Hi,
I'm reading some data from a file. The code below plots the all records in a barplot. How can I plot only records 5 to 10?
By indexing?
Please read the psoting guide and the docs.
chldata - read.table(test3.txt,header=TRUE)
Dear Uwe,
Is this a full installation of R (I guess you have not compiled it
yourself)? If yes, does gcc, perl and friends work?
In the first place, libR.a must be compiled at this point which should
happen automatically.
But it is really hard to say which of the tools fails exactly.
Uwe Ligges
Patrick Giraudoux wrote:
Dear Uwe,
Is this a full installation of R (I guess you have not compiled it
yourself)? If yes, does gcc, perl and friends work?
In the first place, libR.a must be compiled at this point which should
happen automatically.
But it is really hard to say which of the tools
Uwe Ligges ligges at statistik.uni-dortmund.de writes:
Bernd Weiss wrote:
Error in sprintf(gettext(unable to move temp installation '%d' to
'%s'), :
use format %s for character objects
Maybe your disc is full or a package is already in use and one of its
files is locked?
From: Romain Francois
Hello list,
I'd like to know if there is a function that tells in which
package is a
given function.
Something like :
which.package(lda)
[1] MASS
getAnywhere(), as Uwe suggested, will find objects that are in attached
packages. help.search(), as Doug
Ah! See R-2.1.0's version of the manual R Installation and
Administration, Sections Installing R under Windows and in
particular Appendix E The Windows toolset.
Uwe Ligges
OK. Got it... When I started packaging home made functions one year ago,
I did not use compiled code, and just install
From: [EMAIL PROTECTED]
On 23-Apr-05 [EMAIL PROTECTED] wrote:
[snip]
[technical setup snipped]
contrasts(dat$vac_inf) - ginv(m)
gm - aov(y ~ vac_inf, dat)
summary(gm)
Df Sum Sq Mean Sq F value Pr(F)
vac_inf 3 12.1294 4.0431 7.348 0.04190
Residuals4
Dear all,
I am using the optim() function which it stops with the following error
messagge:
error in optim(...) non-finite finite-difference value
I was wondering if somebody might suggest me a way to fix it please.
Thanks in advance to all of you.
Kind regards,
Tom
Just an FYI here as I think that this approach sometimes gets
overlooked.
When one has a very large dataset to read into R from a text file (as
opposed to an external database that can be queried) and in actuality
only needs a subset of sequential records for subsequent analysis, there
are
I just encountered that myself using optim in S-Plus 6.1. For a
certain set of parameter values, the function I wanted to minimize
returned NaN (not a number). To fix this problem, I set a trace with
something like options(trace=99) using, I think, method = CG. This
printed out all the
Hi,
I was trying to extract certain rows from a table. For example, if the value in
the first column is a, d, k, or z, then I would like to extract these
rows. In other words, if the value in the column is also an element in my list
(a,d,k,z), then I would like to extract the row.
I've
Hi,
I was trying to extract certain rows from a table. For example, if the value in
the first column is a, d, k, or z, then I would like to extract these
rows. In other words, if the value in the column is also an element in my list
(a,d,k,z), then I would like to extract the row.
I've
On 4/23/05, Uwe Ligges [EMAIL PROTECTED] wrote:
Is there some way of restarting R, without quitting R? In other words,
I am looking for something like the commands reset (MuPAD) or restart
(Maple).
No, but of course you could write a function that fires up a new
instance of R and quits
[No need to post the same message twice, please.]
Try something like:
dat - data.frame(a=sample(c(a, b, d, k, z), 30, replace=TRUE),
b=runif(30))
subset(dat, a %in% c(a, d, k, z))
Andy
From: Tim Smith
Hi,
I was trying to extract certain rows from a table. For
Hi r-help,
I have some troubles reading fortran binary file(from mm5) in R. Here
is what I have done.
1. Use a fortran subroutine to read this file in R. The subroutine is
as the following.
subroutine freadmm5(filenamet,out2d)
integer iflag,var1,miy,mjx,mkz,mt,z,t
character*4
[EMAIL PROTECTED]
This is a system-generated notification from Computer Centre Ngee Ann
Polytechnic.
Our system has received an email sent with your email address. The email
has been
blocked as it contained a file attachment disallowed under our Email
Security Policy.
As Internet email
Dear R experts,
I would like to explore if and to what extent bootstrapping and
permutation statistics can help me for my research (functional brain
imaging). I am looking for an introductory textbook, rather legible. I
have statistical knowledge, but I am definitely no statistical or
Hello,
may I ask a further question?
I have realized that data -
matrix(scan(file-name), ncol=29) will read the data differently than I
thought, i.e., (4,1) is the first column, (17,1) is the second column, and
(1,1) is the third and so on by this code - please see the data below.
Therefore,
You just need to try harder in reading the documentation. Try:
data - matrix(scan(file-name), ncol=29, byrow=TRUE)
Andy
From: Yoko Nakajima
Hello,
may I ask a further question?
I have realized that data -
matrix(scan(file-name), ncol=29) will read the data
differently than I
If I understand your problem, you are computing the difference between
your data and the quantiles of a standard gaussian variable -- in
other words, the difference between the data and the red line, in the
following picture.
N - 100 # Sample size
m - 1# Mean
s - 2# dispersion
x
The command
help.search(density)
which you should have tried if you had read the posting guide, returns,
among others, kde2d, in package MASS.
-- Vincent
On 4/22/05, Bernard Palagos [EMAIL PROTECTED] wrote:
hello
sorry for my english
I would like estimate density for multivariate
Hi - Just for the record -
Another option for the less command line orientated is
http://downthemall.mozdev.org/
Which is an extension to firefox that allows you to download all or only
selected links on a webpage. You can tick or untick whichever links you
want. Quite convenient.
Cheers
--
Hi,
I ran the following model using nlme:
model2-lme(log(malrat1)~I(year-1982),random=~1|Continent/Country,data=wbmal10)
I'm trying to run a Poisson GlMM to avoid the above transformation but I
don't know how to specify the model using lmer in the lme4 library:
Luis Fernando Chaves asks:
: -Original Message-
: From: [EMAIL PROTECTED]
: [mailto:[EMAIL PROTECTED] On Behalf Of Luis
: Fernando Chaves
: Sent: Sunday, 24 April 2005 11:29 AM
: To: [EMAIL PROTECTED]
: Subject: [R] A question on the library lme4
:
:
: Hi,
:
: I ran the following
47 matches
Mail list logo