If the matrices are not all the same size, then the order of
computation will make a difference. simple example: A is 1xn, B is
nx1, C is 1xn.
A(BC) takes n^3 multiplies, while (AB)C requires 2n.
albyn
Quoting Todd Schneider todd.w.schnei...@gmail.com:
Hi all,
I am looking for a
On Mon, Oct 26, 2009 at 12:19 PM, Ordaz, Sarah ord...@upmc.edu wrote:
Hi,
I'd like to clarify my previous posting. See below - my updates are noted
with *s
Thanks,
Sarah Ordaz
ord...@upmc.edu
Hello,
I am a newbie to the lattice package in R, and I'm trying to make a plot
using the
On Fri, Oct 23, 2009 at 1:54 AM, Jim Lemon j...@bitwrit.com.au wrote:
On 10/23/2009 06:07 AM, Lasse Kliemann wrote:
I wish to save a scatter plot comprising approx. 2 million points
in order to include it in a LaTeX document.
Using 'pdf(...)' produces a file of size about 20 MB, which is
2009/10/27 Joel Fürstenberg-Hägg joel_furstenberg_h...@hotmail.com:
Hello everybody,
I'm using the lattice package and the xyplot to make several graphs like
below. However, I can just print the three grouped plots onto one page as I'm
putting them into a pdf-file, which gives me a huge
Windows XP 32 bit machine.
allocated 2gb as memory for R 2.9.2
data will be 20mb and was running calculations on the data .
but the r crashes with out any warning or so.
--
View this message in context:
http://old.nabble.com/R-crashes-tp26110355p26124740.html
Sent from the R help mailing list
On Thu, 29-Oct-2009 at 05:22PM -0400, zrl wrote:
| Hi,
|
| I am wondering how to control the label position on panels of trellise plot.
| I have 6 panels, I want to all x label positions at the bottom of each
| panel. How should I do that?
| Similar, how should I control y label positions of
On 10/30/2009 05:17 AM, premmad wrote:
Windows XP 32 bit machine.
allocated 2gb as memory for R 2.9.2
can you try again whatever you do with 2.10.0
data will be 20mb and was running calculations on the data .
but the r crashes with out any warning or so.
That is a bit better, but nowhere
hello all
I am using the R package Zelig for some tobit regression with robust
standard errors.
I have got R version 2.9.2 (2009-08-24)
and Zelig Version: 3.4-5
when i do demo(robust)
It ends like this way
data(coalition)
# Fit the model with robust standard error
user.prompt()
Press return
Hi R users:
I got this code to generate a graphic for each member of a lists.
list1-list(A=data.frame(x=c(1,2),y=c(5,6)),B=data.frame(x=c(8,9),y=c(12,6)))
names1-names(list1)
sapply(1:length(list1),function(i)
with(list1[[i]],plot(x,y,type=l,main=paste(Graphic of,names1[i]
Is there a more
The 'Gower' metric is one that is commonly used.
R
On Thu, Oct 29, 2009 at 10:22 AM, karuna m m_karuna2...@yahoo.com wrote:
I am doing hierarchical clustering with cluster package. I couldnot find
similarity measures like matching coefficient, Jaccard coefficient and sokal
and sneath.
The 'Gower' metric is one that is commonly used.
Rod
On Thu, Oct 29, 2009 at 10:22 AM, karuna m m_karuna2...@yahoo.com wrote:
I am doing hierarchical clustering with cluster package. I couldnot find
similarity measures like matching coefficient, Jaccard coefficient and sokal
and sneath.
2009/10/30 hadley wickham h.wick...@gmail.com:
I read anything that mentions
ggplot2 no matter where it is.
... one should hope this statement only applies to the Internet
though, does it? Please do share your regexp if it's not the case.
:)
baptiste
On 10/30/09, Luna Laurent luna.laurent@gmail.com wrote:
Could anybody tell me how to test for stationarity in time series?
http://www.rseek.org/?cx=010923144343702598753%3Aboaz1reyxd4q=stationarity+time-seriessa=Search+functions%2C+lists%2C+and+morecof=FORID%3A11
Liviu
Le 30 oct. 2009 à 02:21, R_help Help a écrit :
Ok. I have the following likelihood function.
L - p*dpois(x,a)*dpois(y,b+c)+(1-p)*dpois(x,a+c)*dpois(y,b)
where I have 100 points of (x,y) and parameters c(a,b,c,p) to
estimate. Constraints are:
0 p 1
a,b,c 0
c a
c b
I construct a
Hi;
Does anybody aware of R implementation for AIB (Agglomerative Information
Bottleneck) by Slomni and Tishby?
in which package should i find it?
Thanks
_
Hmm, so if read correctly you want to remove exactly duplicated rows. So
maybe try the following to begin with.
duplicated(newdf[ , c(id, loc, clm)])
[1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE
TRUE TRUE
Then you can remove the duplicated rows before proceeding
Dear R-users,
my present problem is related to lattice.
I would like to put some text in each panel, namely a p-value.
I therefore wrote a simple panel function as reported here below.
I'm able to write the value in each panel at the maximum value of y for
each panel,
but this obviously
Readers,
I have the following command:
expression(A[1]B~%)
but that causes an error so I changed to:
expression(A[1]~B~%)
but the result is too much whitespace between subscript and B. Is
there a way to reduce this whitespace?
Yours,
rhelp at conference.jabber.org
r 251
mandriva 2008
I'm having a problem where I have to apply a function to a subset of a
variable, where the subset is defined by the n nearest neighbours of a
second variable.
Here's an example applied to the 'iris' dataset:
$ head(iris)
Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1
Hi Karuna,
You can try vegdist() function in the vegan package. The function
computes dissimilarity indices that are useful for or popular with community
ecologists. All indices use quantitative data, although they would be named
by the corresponding binary index, but you can calculate the
Hi,
also interested...
If you are checking for non-normal behaviour, first let us define normal
behaviour: only small temperature changes and no steep ramps?
If so, maybe you can make an rolling average of the last x points, and check
if the following point deviates more than ... ?
Or is it
Thanks for all the suggestions.
Cut + plotmeans (gplots) solved my problem.
groups - cut(t[,1],breaks=5)
plotmeans(t[,2]~groups)
Thanks!
Federico
--- El jue, 29/10/09, Jim Lemon j...@bitwrit.com.au escribió:
De: Jim Lemon j...@bitwrit.com.au
Asunto: Re: [R] x-y plot as an histogram
Para:
Hi Luna,
tseries package provides various stationarity and unit root
tests including Augmented Dickey-Fuller, Phillips-Perron, and KPSS.
Alternative implementations of the ADF and KPSS tests are in the
urca package, which also includes further methods such as
Elliott-Rothenberg-Stock,
e-letter wrote:
Readers,
I have the following command:
expression(A[1]B~%)
but that causes an error so I changed to:
expression(A[1]~B~%)
but the result is too much whitespace between subscript and B. Is
there a way to reduce this whitespace?
Use * rather than ~
Uwe Ligges
Yours,
On Fri, 30 Oct 2009 10:28:49 +0100 Karl Ove Hufthammer k...@huftis.org
wrote:
$ (pos=which(order(abs(iris$Sepal.Length-x)) %in% 2:6))
This should of course be:
(pos=order(abs(iris$Sepal.Length-x))[2:6])
--
Karl Ove Hufthammer
__
Dear All,
I am struggling with extracting data from a data frame:
x=data.frame(a=1:11,b=100:110)
What I want is a list/vector in this sence: 1 100 2 101 3 102...
For single rows, this works fine:
as.matrix(x)[1,]
For, say 2 rows, this works fine:
z=c(as.matrix(x)[1,],as.matrix(x)[2,])
But
On Fri, 30 Oct 2009 10:42:12 +0100 Christian Lerch t@gmx.net
wrote:
I am struggling with extracting data from a data frame:
x=data.frame(a=1:11,b=100:110)
What I want is a list/vector in this sence: 1 100 2 101 3 102...
y=t(as.matrix(x))
as.vector(y)
Happy to help.
--
Karl Ove
On Fri, Oct 30, 2009 at 2:10 AM, Ottorino-Luca Pantani
ottorino-luca.pant...@unifi.it wrote:
Dear R-users,
my present problem is related to lattice.
I would like to put some text in each panel, namely a p-value.
I therefore wrote a simple panel function as reported here below.
I'm able to
Nice workaround :-). Thank you.
Best,
Christian
Original-Nachricht
Datum: Fri, 30 Oct 2009 10:54:49 +0100
Von: Karl Ove Hufthammer k...@huftis.org
An: r-h...@stat.math.ethz.ch
Betreff: Re: [R] data.frame extracting data row-wise
On Fri, 30 Oct 2009 10:42:12 +0100 Christian
Hi,
This is a newbie question.
I would to be able to convert annual time series of flow data into quarterly
data. I wonder if there is any existing R-function which permits to do it? In
what package ?
I the archive, i found that some poeple speak about tempDis package for
performing time
sapply(names(list1), function(.data.){
with(list1[[.data.]], plot(x, y, type='l', main=paste(Graphic of, .data)))
})
On Fri, Oct 30, 2009 at 3:41 AM, Kenneth Roy Cabrera Torres
krcab...@une.net.co wrote:
Hi R users:
I got this code to generate a graphic for each member of a lists.
Hello,
ascii package has been updated to version 0.3.2.
ascii is a R package for writing asciidoc, txt2tags, sphinx or org documents
with embeded R code. See http://eusebe.github.com/ascii/ for some examples.
News since version 0.2 are:
- new sphinx driver and output (see
You could lapply over the names:
junk - lapply(names(list1), function(nm)
with(list1[[nm]], plot(x, y, type = l, main = paste(Graphic
of, nm
which is a slight improvement though that still involves dealing with
a list and its names separately.
On Fri, Oct 30, 2009 at 3:41 AM, Kenneth
Hi,
I want to parallelize some computations when it's possible on multicore
machines.
Each computation produces a big objects that I don't want to store if
not necessary: in the end only the object that best fits my data have to
be returned. In non-parallel mode, a single gloabl object is
See:
https://stat.ethz.ch/pipermail/r-sig-finance/2009q3/004677.html
On Fri, Oct 30, 2009 at 6:09 AM, Axel Leroix axel.ler...@yahoo.fr wrote:
Hi,
This is a newbie question.
I would to be able to convert annual time series of flow data into quarterly
data. I wonder if there is any existing
Hi, everybody.
My name is Ricardo and I'm doing doctorate at USP (Brazil). I have used a
lot the package 'fractal', developed by William Constantine and Donald
Percival, in my research.
However, I have a doubt about the execution of RoverS function. I was
reading an article
Hi,
I am an [R] and lattice beginner.
I an doing an epidemiological study and want to present my data in a
condensed and clear manner as possible.
The example below is fake, but representative for the the data I want to show.
I have a data frame with different categorical data and one column of
Dear gbm users,
When running predict.gbm() on a large dataset (150,000 rows, 300 columns,
500 trees), I notice that the memory used by R grows beyond reasonable
limits. My 14GB of RAM are often not sufficient. I am interpreting this as a
memory leak since there should be no reason to expand memory
Dear Robert,
The differences have to do with diffent scaling defaults.
lrm by default standardizes the covariates to unit sd before applying
penalization. penalized by default does not do any standardization, but
if asked standardizes on unit second central moment. In your example:
x =
Hi there,
when conducting a multidimensional scaling analysis, it is possible to do
this on a correlation matrix?
I calculated a dissimilarity matrix by substracting the corr.matrix from 1.
Then I calculated the distance matrix with d - dist(dissimilarity matrix)
and ran the isoMDS function.
Is there anyway to make the lines in the dispersion command come forward in
a plot and allow the fill in the dispersion parameter be transparent so all
of the lines I am using to note confidence intervals are shown?
---
Joe King, M.A.
Ph.D. Student
Dear all,
Does anyone know if the PLM package (to run Panel Data Analysis) accepts
quarterly data?
The package vignette and documentation only use annual data -and the only time
index available seems to work for years.
José
Mr José Luis Iparraguirre
Senior Research Economist
Hi,
I could not find a function that returns a time in hours and minutes
given seconds as argument. I looked also in the chron package.
I would need something like this:
mytime(61) #61 are seconds
[1] 00:01:01
mytime(7385) #7385 are seconds
[1] 02:03:05
Am I missing something? or should I
Hello R users:
I'm trying now to apply the package strucchange to see whether there is
a structural change in linear regression. I have noted the following
problem that arises in my case with recursive-based CUSUM: generic
function recresid() in efp() generates an error, since (probably) it
Dear friends,
I will be very happy if anyone tell me the way to change work directory
permanently?
I mean not use the function setwd() which can only change temporary, when
you close the console, it will the old directory.
Sys.setenv(R_USER = '') also doesn't work.
[[alternative HTML
Try this:
library(chron)
times(100/(24*60*60)) # 100 seconds
[1] 00:01:40
On Fri, Oct 30, 2009 at 5:37 AM, Stefano Leonardi
stefano.leona...@unipr.it wrote:
Hi,
I could not find a function that returns a time in hours and minutes given
seconds as argument. I looked also in the chron
In my .Rprofile I have the following functions which display the
current directory in the main R window title bar,
and modify base::setwd() to keep this up to date. I like this because I
can always tell where I am in the file system.
cd - function(dir) {
base::setwd(dir)
Hi All:
I have a data frame called all_corn. This has 31 columns. The first column is a
character key. The next 15 columns (stat1,stat2,...,stat15) are the statistical
forecast. The last 15 columns (sls1,sls2,...,sls5) are actual sales.
I want to calculate textbook tracking signal and cuulative
On Oct 29, 2009, at 5:23 PM, PDXRugger wrote:
Terrific help thank you.
dupbuild-aggregate(DF$Acres, list(Bldgid), sum)
This line worked best.
Now im going to challenge everyone (i think?)
Consider the following:
Acres-c(100,101,100,130,156,.5,293,300,.09,100,12.5)
On Wed, Oct 28, 2009 at 8:19 PM, Bill Gillespie bi...@metrumrg.com wrote:
I currently use lattice functions to produce multiple pages of plots using
the layout argument to specify the number of rows and columns of panels,
e.g.,
xyplot(price ~ carat | clarity, diamonds, layout = c(2, 2))
Date: Fri, 30 Oct 2009 09:29:06 +0100
From: Christophe Dutang duta...@gmail.com
Subject: Re: [R] [R-SIG-Finance] Fast optimizer
To: R_help Help rhelp...@gmail.com
Cc: r-help@r-project.org
Ok. I have the following likelihood function.
L -
Hi
r-help-boun...@r-project.org napsal dne 30.10.2009 14:40:06:
Hi All:
I have a data frame called all_corn. This has 31 columns. The first
column is
a character key. The next 15 columns (stat1,stat2,...,stat15) are the
statistical forecast. The last 15 columns (sls1,sls2,...,sls5) are
In order to instrument functions for optimization or DE solving etc., I
want to set up a private ensemble of working data somewhere in my
workspace that won't collide with the problem-data. I can do this
already with named objects within .GlobalEnv, but I'd like to have
functions like
On Oct 30, 2009, at 9:11 AM, Erik Svensson wrote:
Hi,
I am an [R] and lattice beginner.
I an doing an epidemiological study and want to present my data in a
condensed and clear manner as possible.
The example below is fake, but representative for the the data I
want to show.
I have a data
On Oct 30, 2009, at 7:02 AM, Joe King wrote:
Is there anyway to make the lines in the dispersion command come
forward in
a plot and allow the fill in the dispersion parameter be transparent
so all
of the lines I am using to note confidence intervals are shown?
Got code?
Hi All,
Once again , would like to know about your suggestions regarding the
integration of PHP and R.I have read **
*4. R Web Interfaces* here:
http://cran.r-project.org/doc/FAQ/R-FAQ.txt
It says about many tools.Here my requirement is :--
I have a tool developed in PHP.I want to
Hi:
In general, how to I cast a character to the class that I am trying to change.
For example, if I have a data frame df1. df1 has a column x. suppose I want to
a substring of x (the first 3 chars). Then I want to do something like
df1$new = substring(of x)
Example
Data frame df1
x
abcd
efgh
On Wed, 28 Oct 2009, Kingsford Jones wrote:
I thought I'd share this with the list since it appears to provide a
quick fix to some memory problems, and I haven't see it discussed in
relation to R.
Seems you didn't look in the places mentioned in the posting guide, fo
this is an actual FAQ!
On Oct 30, 2009, at 11:24 AM, Vadlamani, Satish {FLNA} wrote:
Hi:
In general, how to I cast a character to the class that I am trying
to change.
For example, if I have a data frame df1. df1 has a column x. suppose
I want to a substring of x (the first 3 chars). Then I want to do
Dear R users,
this is a follow up of this message
http://tolstoy.newcastle.edu.au/R/e6/help/09/05/13897.html
I'm reproducing the core of it for convenience.
//
/ data(Oats, package = MEMSS) /
/ tp1.oats - xyplot(yield ~ nitro | Variety + Block, /
/ data = Oats, /
/
On 10/30/2009 9:30 AM, Michael Friendly wrote:
In my .Rprofile I have the following functions which display the
current directory in the main R window title bar,
and modify base::setwd() to keep this up to date. I like this because I
can always tell where I am in the file system.
cd -
How do you do a double loop through a matrix or data frame , and
within each iteration , do a calculation and append it to a new,
second data frame?
(So, if your original matrix or data frame is 4 x 5 , then 20
calculations are done, and the new data frame, which had 0 rows to
start with, now has
On Fri, Oct 23, 2009 at 6:05 PM, Ravi Varadhan rvarad...@jhmi.edu wrote:
Frank,
I have heard this (i.e. only head-to-head comparisons are valid) and various
other folklores about AIC and BIC based model selection before, including
one that these information criteria are only applicable for
OK, here's the first-cut code, which does execute properly. What I
would like to do, simply for 'cleanliness', is to assign nothing to
the variable l , rather than assigning seq(1,dim(data)[timedim]).
Then the result of all the get(unlist(dims)) operations below will be
func[j,k,]
Where
At least in the Unix version of the ``Rscript`` front-end to R (but
oddly enough not ``R CMD BATCH``) you could put something like this in
test.R:
cat(Save workspace image? [y/n]: )
answer - tolower( scan(stdin, , n=1) )
if (substr(answer, 1, 1) == y) {
q(yes)
} else {
q(no)
}
Which
On Oct 30, 2009, at 11:46 AM, Robert Wilkins wrote:
How do you do a double loop through a matrix or data frame , and
within each iteration , do a calculation and append it to a new,
second data frame?
(So, if your original matrix or data frame is 4 x 5 , then 20
calculations are done, and the
premmad wrote:
Windows XP 32 bit machine.
allocated 2gb as memory for R 2.9.2
data will be 20mb and was running calculations on the data .
but the r crashes with out any warning or so.
Ahh, well 20 MB is not huge data for R, so there should be no problem with
processing it. The
Thanks a lot for the hints. The sapply method may work for me.
But how can I extract just the .grad expression from the return value
of the deriv function? (and secondly store in a matrix of expressions?)
Thanks again.
__
R-help@r-project.org mailing
Is there any R package that implements the same capability of MatLab toolbox
called SimBiology ?
We are expecially interested in protein-protein interactions and network
analysis.
As far as I know SimBiology implements a system of ODEs reflecting the kinetic
chemical reactions.
We would be more
Hi,
Two questions:
1 - Say I have average speed and directions for tide and I would like to
plot them on a polar plot, but with different colors so I can indicate the
two directions. I'm using polar.plot from the plotrix library. How can I add
a second b and dir.b series to a polar.plot?
Hi,
I've been trying to work with the gnls() function in the nlme package. My
decision to use gnls() was so that I could fit varPower and such to some of
the data. However, in working with a small dataset, I've found that the
results given by gnls() don't seem to make any sense and they
Thank you for the response, Prof Ripley. Some comments below
On Fri, Oct 30, 2009 at 9:27 AM, Prof Brian Ripley
rip...@stats.ox.ac.uk wrote:
On Wed, 28 Oct 2009, Kingsford Jones wrote:
I thought I'd share this with the list since it appears to provide a
quick fix to some memory problems, and
Which OS and how do you start R?
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of Huidong TIAN
Sent: Friday,
Hi there,
when conducting a multidimensional scaling analysis, it is possible to do
this on a correlation matrix?
I calculated a dissimilarity matrix by substracting the corr.matrix from 1.
Then I calculated the distance matrix with d - dist(dissimilarity matrix)
and ran the isoMDS function.
Julia:
I'm trying now to apply the package strucchange to see whether there is a
structural change in linear regression. I have noted the following problem
that arises in my case with recursive-based CUSUM: generic function
recresid() in efp() generates an error, since (probably) it cannot
Hello,
We (two mac users) have been attempting to install rgdal from
http://www.kyngchaos.com/software:frameworks;, given that it is not available
as a binary on the CRAN (binaries) of the Package Installer.
I have also tried to solve this problem by looking on the net for an old
question,
Hello,
I have a data file originally as .mat file from matlab that I imported
in R using the R.matlab package (readMat). This loaded a list
containing amongst other things a 4x4 array of rainfall data that is
quite heavy (around 40 Mb, it took 10 minutes to load on my computer).
I tried
Dear all
Is it possible to programmatically detect which commands have
generated a graphic to the x11() graphics device? When in front of the
computer, it is easy to see that after a command---say, plot(1:10)---a
graphics window opens/activates and displays a graphic. But is there a
way to detect
On Oct 30, 2009, at 1:51 PM, Pablo Alvarez wrote:
Hello,
We (two mac users) have been attempting to install rgdal from http://www.kyngchaos.com/software:frameworks
, given that it is not available as a binary on the CRAN (binaries)
of the Package Installer.
I have also tried to solve this
Can you at least show us what you were doing. It sounds like you were
trying to 'load' a text file, not something that was saved with
'save'. Here is what I get:
load('trackV3.r')
Error: bad restore file magic number (file may be corrupted) -- no data loaded
In addition: Warning message:
file
Hi all,
I am trying to estimate a simple logit model.
By using MLE, I am maximizing the log likelihood, with optim().
The thing is, each observation has different set of choice options, so I need a
loop inside the objective function,
which I think slows down the optimization process.
The
On 10/30/2009 3:13 PM, Liviu Andronic wrote:
Dear all
Is it possible to programmatically detect which commands have
generated a graphic to the x11() graphics device? When in front of the
computer, it is easy to see that after a command---say, plot(1:10)---a
graphics window opens/activates and
Heres my code:
era1 - seq(1.5,5,.25)
plot(era1,yhyp1$pe, xlab=ERA, ylab=Probability of Winning Cy
Young,type=l, col=black, lwd=2)
dispersion(era1,yhyp1$pe,y1upper,y1lower,type=l,
fill=blue,arrow.cap=0.01,intervals=FALSE)#Requres plotrix
dispersion(era1,yhyp2$pe,y2upper,y2lower,type=l,
(R 2.9.2 on Windows XP)
Turns out that
my.data.frame[c(1,3)] - list(NULL, NULL)
and
my.data.frame[c(3,1)] - list(NULL, NULL)
are not the same. Is this expected?
Example:
d0 - structure(list(Year = c(2009L, 2009L, 2009L), Season = c(1L,
1L, 1L), Series = c(1L, 1L, 1L), Mst =
On Oct 30, 2009, at 3:49 PM, Joe King wrote:
Heres my code:
era1 - seq(1.5,5,.25)
plot(era1,yhyp1$pe, xlab=ERA, ylab=Probability of Winning Cy
Young,type=l, col=black, lwd=2)
This plots.
library(plotrix) # should go here.
dispersion(era1,yhyp1$pe,y1upper,y1lower,type=l,
Dear R users:
I need to compare character strings stored in 2 separate data frames. I need
an exact match, so finding a in animal is no good.
I've tried regexpr, match, and grepl, but to no avail.
Anybody know how to accomplish this?
Ben
--
View this message in context:
I would like to increase the precision/accuracy of R. That is, I'm
dealing with numbers exceedingly close to 1, and I would like to
increase the number of significant digits used in computation. In
matlab, you can use vpa() to accomplish this. I'm wondering if there's
an equivalent or a
I am fitting a model using zeroinfl() and it runs without errors,
returning results that are generally consistent with my hypotheses.
One of my variables is percent black (pblack). This variable was highly
significant in some of the other count models I ran on the way to my
current
I would like to increase the precision/accuracy of R. That is, I'm
dealing with numbers exceedingly close to 1, and I would like to
increase the number of significant digits used in computation. In
matlab, you can use vpa() to accomplish this. I'm wondering if there's
an equivalent or a
David,
You are correct. I think the frist two assumptions can be thrown out and
only the latter two (c,d) can be considered. So how would i combine Acres
for matching Bldgids based on assumptions c,d?
David Winsemius wrote:
On Oct 29, 2009, at 5:23 PM, PDXRugger wrote:
Terrific
Unless this is a homework problem, you would be much better off using
glm().
Giovanni
Date: Fri, 30 Oct 2009 12:23:45 -0700
From: parkbomee bbom...@hotmail.com
Sender: r-help-boun...@r-project.org
Importance: Normal
Precedence: list
--Boundary_(ID_/D+lL9iK1qLhrkPBeoxH+Q)
Dear list,
my problem seems to be primarily a statistical one, but maybe there is a
misspecification within R (and hopefully a solution).
I have two groups with two measured variables as training data. According to
the variables, the groups differ totally. I know that this is a very easy
Hi Greg,
Others on this list are much more knowledgeable than I am, and I'm
sure they will correct me if I'm wrong.
My understanding is that R uses double precision (see
http://en.wikipedia.org/wiki/Double_precision_floating-point_format),
which has the following implication:
x -
Thank you.
But I'd prefer using a written function which allows me more flexible model
specification.
Later on, I could have random parameters.
So I want to know if there is any more efficient way so that I can speed it up.
Date: Fri, 30 Oct 2009 16:10:29 -0600
To: bbom...@hotmail.com
CC:
On 30/10/2009 5:04 PM, Greg Michaelson wrote:
I would like to increase the precision/accuracy of R. That is, I'm
dealing with numbers exceedingly close to 1, and I would like to
increase the number of significant digits used in computation. In
matlab, you can use vpa() to accomplish this.
I think he means something like packages Rmpfr and gmp provide.
However, careful numerical analysis (and working with x-1 or log(x)
for numbers very near one) means that such packages are rarely
necessary (although they can be very convenient).
It all depends on the calculations: vague
Pablo Alvarez-2 wrote:
Hello,
We (two mac users) have been attempting to install rgdal from
http://www.kyngchaos.com/software:frameworks;, given that it is not
available as a binary on the CRAN (binaries) of the Package Installer.
The GDAL library contains an impressive collection
Hi R users
I want to do a multiple linear regression with R
With a normal model (Y=Cste+A1*X1+A2*X2) the program would be
lm(Y ~X2+X2)
My model is Y=Cste+A1*X1+A2*X2 with the constraint A1=-A2
What is the program for such a model?
Best regards
--
View this message in context:
I had earlier downloaded what I thought to be the most recent version
of rgdal (rgdal_0.6-18.tgz). After checking to see that all of the
other dependencies in the configure files were met, I finally got
around to checking the version of the rgdal and found my problem. It
was different (and
Sorry there was a mistake:
Hi R users
I want to do a multiple linear regression with R
With a normal model (Y=Cste+A1*X1+A2*X2) the program would be
lm(Y~X1+X2)
My model is Y=Cste+A1*X1+A2*X2 with the constraint A1=-A2
What is the program for such a model?
Best regards
--
View this
1 - 100 of 108 matches
Mail list logo