Dear List
Forgive me for the slightly off topic query
Force.com and Salesforce have many (1009) apps at
http://sites.force.com/appexchange/home for cloud computing for
businesses, but very few forecasting and statistical simulation apps.
Example of Monte Carlo based app is here
On Wed, 2010-11-10 at 13:38 +0530, Ajay Ohri wrote:
Dear List
Forgive me for the slightly off topic query
In general OT is OK if it relates to R, but why should we forgive you
for posting to three R lists at once! This would seem to be OT for
R-Devel and so far OT for R-Packages it makes me
No I am not on the payroll for these people.
or any relation of any kind ;) I just like cloud computing thats all
Appreciate your kind consideration Mr Simpson- I was not sure of the
list so mailed all 3- .
The subject line has been modified-hopefully there should be a
technical answer to a
Hi Nikos.
There is quite a bit going on here, both in the code and in your
terminology.
You should really consider reading An introduction to R that comes with
your R installation.
A few pointers though:
* in R speak, you have nowhere declared 2 global matrices: it is not
completely clear why
Dear R users,
Here is my problem:
I have an array with at least four dimensions:
dim(myArray)
[1] 20 17 3 6
I'd like to apply a function to each occurrence of the matrix (3x6)
defined by the last two dimensions. This interpolation function always
return a matrix of the same dimensions
Hi Ivan,
thanks for the hint. Now it works as expected.
Sincerely,
/steffen
Am 10.11.2010 10:10, schrieb Ivan Calandra:
Hi,
The xlim and ylim arguments should be given the extremes of the range,
not a range:
plot(x=NULL,
y=NULL,
xlim=c(1, 10), ## with c()
ylim=v(1e-9, 1e-3), ## with c()
Sou usuário do R há bastante tempo e,
gostaria de participar dessa lista de dis-
cussão sobre o programa. Na verdade
sou apenas um iniciante e preciso ap-
prender muito para poder discutir.
Obrigado.
Doki.
[[alternative HTML version deleted]]
The other possibilities are:
(1) you're missing a necessary interaction term
(2) one of the variables affecting output just isn't in your data set.
(3) you need to transform 'hosp_days' or 'age' -- are those the only
two continuous variables? Might be worth trying to plot 'em versus #
of
You could use the following to achieve your objective. To start with
?readLines
?strsplit
?for
?ifelse
As you try, you may receive more specific answers for the issues you come up
with.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Date: Wed, 10 Nov 2010 12:47:10 +0200
From: hka...@gmail.com
To: R-help@r-project.org
Subject: [R] Installing R and an editor on a USB drive
Hi,
I have adviced my students to install R and an editor on a USB drive for
working in the
Hey,
I don't know whether this helps, but an editor for Windows with code
highlighting for any programming language (including R) that can be
installed on a pen drive (USB stick) is Notepad++
(http://notepad-plus-plus.org/) and with the NppToR application
(http://npptor.sourceforge.net/)
On Nov 10, 2010, at 3:47 AM, Nick Sabbe wrote:
Hi Nikos.
There is quite a bit going on here, both in the code and in your
terminology.
You should really consider reading An introduction to R that comes
with
your R installation.
A few pointers though:
* in R speak, you have nowhere
Hello Nasrin,
Please attach how each of your files look like (their first few rows), so we
could understand why the rbind doesn't work.
In general, be sure to keep the r-help e-mail also corresponded so others
might help if I don't know the answer (or if they answer before me).
Best,
Tal
Hi Sachin,
I guess there are several different possibilities that are more or less handy
depending on your data:
- lists were mentioned already, and I think they are the most natural
representation of ragged arrays. Also very flexible, e.g. you can introduce more
dimensions. But they can
On 10/11/2010 12:12 AM, Michael Bedward wrote:
Hello Sachin,
You have a ragged array and you can easily store this as a list of vectors...
x- list(c(0,0,1,1), c(1,3,5), 4, c(7, -1, 8, 9, 10, 6))
The only gotcha with this is that you will then need to use double
brackets for the first index
On Nov 10, 2010, at 5:13 AM, Jim Lemon wrote:
On 11/10/2010 07:55 PM, Santosh Srinivas wrote:
Basic question ... checked the help page but the only answer was to
use
paste!
Is there any way to format as %?
sprintf(%1.0f %%,10)
[1] 10 %
sprintf(%1.0f%%,10)
[1] 10%
degree =
On 11/10/2010 07:55 PM, Santosh Srinivas wrote:
Basic question ... checked the help page but the only answer was to use
paste!
Is there any way to format as %?
degree = c(0.20,0.5)
degree
[1] 0.2 0.5
print(degree)
[1] 0.2 0.5
Hi Santosh,
If you want to scale proportions to percentages,
Dear Yves,
You may not need to do more than set the dim attribute correctly:
dim (test) - c (dim (myArray) [c (3 : 4, 1 : 2)]
or
dim (test) - c (dim (myArray) [c (4 : 3, 1 : 2)]
Claudia
--
Claudia Beleites
Dipartimento dei Materiali e delle Risorse Naturali
Università degli Studi di Trieste
Dear mailing list readers!
Using R and the plot function I stumbled over this little issue:
plot(x=NULL,
y=NULL,
xlim=range(1:10),
ylim=range(1e-9:1),
log=y)
produces an plot empty plot, where the y-axis is in the specified range.
Changing ylim to
ylim=range(1e-9:1e-3)
Basic question ... checked the help page but the only answer was to use
paste!
Is there any way to format as %?
degree = c(0.20,0.5)
degree
[1] 0.2 0.5
print(degree)
[1] 0.2 0.5
--
Thanks
Hello,
I have a tab limited text document with multiple lines as mentioned below,
#FILE FORMAT
#Book booknameauthor publisher pages
#CD namecontent
Here is a start:
# read the input file
input - readLines('/tempxx.txt')
# process the file starting at each Book
result - lapply(which(grepl(^Book, input)), function(.line){
+ contents - NULL # initialize
+ name - strsplit(input[.line], '\t')[[1]][2] # book name
+ # process
Hi,
I am running the examples in page 70 of the ff package document, but it
failed with the following error
cat(let's create some ff objects\n)
let's create some ff objects
n - 8e3
a - ff(sample(n, n, TRUE), vmode=integer, length=n, filename=d:/tmp/a.ff)
b - ff(sample(255, n, TRUE),
Steffen,
Consider the input you have given to range:
1e-9:1
[1] 1e-09 1e+00
versus
1e-9:1e-3
[1] 1e-09
hth, Ingmar
On Wed, Nov 10, 2010 at 10:00 AM, Steffen Uhlig
steffen.uh...@htw-saarland.de wrote:
Dear mailing list readers!
Using R and the plot function I stumbled over this little
Hi,
The xlim and ylim arguments should be given the extremes of the range,
not a range:
plot(x=NULL,
y=NULL,
xlim=c(1, 10), ## with c()
ylim=v(1e-9, 1e-3), ## with c()
log=y)
HTH,
Ivan
Le 11/10/2010 10:00, Steffen Uhlig a écrit :
Dear mailing list readers!
Using R and
Date: Tue, 9 Nov 2010 18:25:18 -0800
From: djmu...@gmail.com
To: paul.rhee...@up.ac.za
CC: r-help@r-project.org
Subject: Re: [R] repeatedrepeated measures in ANOVA or mixed model
Hi:
This sounds like a 'doubly repeated measures problem'. Are
From: santosh.srini...@gmail.com
To: karthick.laksh...@gmail.com; r-help@r-project.org
Date: Wed, 10 Nov 2010 16:00:26 +0530
Subject: Re: [R] Parsing txt file
You could use the following to achieve your objective. To start with
?readLines
Hi,
I have adviced my students to install R and an editor on a USB drive for
working in the computer class. With R everything works fine following these
instructions: http://personal.bgsu.edu/~mrizzo/Rmisc/usbR.htm.
But several editors (e.g., Tinn-R and WinEdt) require administrator rights.
I
Nelson,
Esta lista é em inglês.
Para participar acesse https://stat.ethz.ch/mailman/listinfo/r-help
2010/11/10 Nelson Doki mnd...@gmail.com
Sou usuário do R há bastante tempo e,
gostaria de participar dessa lista de dis-
cussão sobre o programa. Na verdade
sou apenas um iniciante e
Tinn R does not require administrator privileges. I've installed it on my USB
stick and it works fine.
--
View this message in context:
http://r.789695.n4.nabble.com/Installing-R-and-an-editor-on-a-USB-drive-tp3035839p3036030.html
Sent from the R help mailing list archive at Nabble.com.
thanks Gabor
There is no way i can come up with your code. I guess if i am using loop in
R, I am doing it wrong.
--
View this message in context:
http://r.789695.n4.nabble.com/Help-make-this-simpler-count-business-day-tp3034621p3036042.html
Sent from the R help mailing list archive at
Is it possible to send data from an executing R script, to external R script
files, on linux in this case, get that external R script to execute, and pass
results back to the central script? If so what commands should I be looking
at? I've googled for R and external, stuff like that but no
Hello,
I have a short question about the prcomp function. First I cite the
associated help page (help(prcomp)):
Value:
...
SDEV the standard deviations of the principal components (i.e., the square
roots of the eigenvalues of the covariance/correlation matrix, though the
calculation is actually
Hi Dennis,
Yes indeed. But when I apply panel.text then the text is not nicely centered
in each bar.
KR,
Ashraf
On Wed, Nov 10, 2010 at 3:31 AM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
Did you mean
panel.text(x,y/2,label = round(y,3),cex=1)
??
HTH,
Dennis
On Tue, Nov 9, 2010
Unfortunately, bootcov is not meant to operate on fit objects produced by
fit.mult.impute. bootcov gets there too late in the process and does not
know how to penalize for imputation.
Frank
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in
On Tue, Nov 9, 2010 at 8:49 PM, Gabor Grothendieck
ggrothendi...@gmail.com wrote:
On Tue, Nov 9, 2010 at 11:49 AM, cameron raymond...@invesco.com wrote:
Help make this simpler – count business day
I am a beginner in R and this is my first post
Want to count the day in month. For example
As I recall (my memory may be faulty), Tinn-R installation required admin
privileges, but running does not and the installation does not include the
specific file path in any files.
To make Tinn-R available on a machine where I do not have admin privileges I
have
a) installed onto a machine
It would be possible to call other R scripts using calls to ?system or
?system2 and ?Rscript, provided you formatted what you passed as [args]
and how your scripts handled [args].
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
My parallel code is running slower than my non-parallel code! Can someone
pls advise what am I doing wrong here?
t and tTA are simple matrices of equal dimensions.
#NON PARALLEL CODE
nCols=ncol(t)
nRows=nrow(t)
tTA = matrix(nrow=nRows,ncol=nCols)
require(TTR)
system.time(
for (i in 1:nCols) {
I think PCA decomposes matrix A according to A'A, not to COV (A).
But if A is centered then A'A = (n + 1) COV (A).
So for non-centered A, you want to look at A'A instead:
crossprod(A) %*% evec[,1] / (nrow (A) - 1) - eval [1] * evec [,1]
[,1]
[1,] 0.000e+00
[2,] 0.000e+00
[3,]
Hi Frank,
Good to know. Thanks again for your help!
Kim
On Wed, Nov 10, 2010 at 8:53 AM, Frank Harrell f.harr...@vanderbilt.eduwrote:
Unfortunately, bootcov is not meant to operate on fit objects produced by
fit.mult.impute. bootcov gets there too late in the process and does not
know
On Wed, Nov 10, 2010 at 10:47 AM, Hannu Kahra hka...@gmail.com wrote:
(setq inferior-R-program-name G:/r-2.12.0/bin/i386/rterm.exe)
since R is installed on the G drive. Everything works if R is on G, but when
changing the computer, R is usually on another drive and Emacs cannot find
it. Is
Can you provide a little more information. What operating system are
you using? Have you monitored the CPU and memory utilizations of the
processes? Do you have enough physical memory; e.g., are you paging?
How big are the matrices that you are processing; e.g., str(tTA) and
object.size(tTA).
When I use
Plot(series1)
lines(series2)
The graph will use the y axis scaling for series 1 so some of series
2 is
cut off. How do I control the y axis scaling?
Plot them together as a multivariate series as in the examples
already
pointed to.
...or try
plot(series1, type=l,
From: santosh.srini...@gmail.com
To: r-help@r-project.org
Date: Wed, 10 Nov 2010 19:37:29 +0530
[[elided Hotmail spam]]
My parallel code is running slower than my non-parallel code! Can someone
pls advise what am I doing wrong here?
I cited
I guess this is not OT for r-devel and r-packages -- it's simply spam.
Regards,
Yihui
--
Yihui Xie xieyi...@gmail.com
Phone: 515-294-2465 Web: http://yihui.name
Department of Statistics, Iowa State University
2215 Snedecor Hall, Ames, IA
On Wed, Nov 10, 2010 at 2:19 AM, Gavin Simpson
Thank you Dave for the reply. I used solstudio12.2, and followed the R admin
installation instruction for Solaris. When I run R as the compiled result,
there is a line Platform: sparc-sun-solaris2.10 (32-bit).
Jun
CC=cc -xc99
CPPFLAGS=-I/opt/csw/include
CFLAGS=-O -xlibmieee
F77=f95
I used solstudio12.2 to build 32-bit R. How do you get around the Matrix
segfault?
-Original Message-
From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk]
Sent: Wednesday, November 10, 2010 12:46 AM
To: Zhang,Jun
Cc: 'r-help@r-project.org'
Subject: Re: [R] R-2.12.0 problem on Sparc
zhiji19 zhiji19 at gmail.com writes:
I need to write R code for lim(I+X/n)^n with repeat loop. (note: limit is
n from 1 to infinity, I is identity matrix, X is square symmetric matrix)
Can anyone please provide help with my question.
Here's a more efficient solution. If you *need* to use a
I am very sorry. It's just that it's not my computer that the problem
occured on... It's at work. So I am going to see what distribution it was
when I get there. But I have no chance to know how R was installed... since
it wasn't me who installed it.
But could you please tell me , at least, where
Hello.
I have a general question regarding to clustering of association rules.
According to http://cran.r-project.org/web/packages/arules/vignettes/arules.pdf
4.7 Distance based clustering transactions and associations there is
possibility for creating clusters of association rules.
I do not
you are right. face palm to forehead. sorry for the spam
Websites-
http://decisionstats.com
http://dudeofdata.com
Linkedin- www.linkedin.com/in/ajayohri
On Wed, Nov 10, 2010 at 8:08 PM, Yihui Xie x...@yihui.name wrote:
I guess this is not OT for r-devel and r-packages -- it's simply spam.
For a Pareto distribution, even a truncated one, the inverse CDF method
should be straightforward to implement.
Giovanni Petris
On Tue, 2010-11-09 at 10:50 -0600, cassie jones wrote:
Dear all,
I am trying to simulate from truncated Pareto distribution. I know there is
a package called
Hi, does anyone know if there is a way to easily extract p-values from
the regsubsets() function?
Thanks,
James Stegen
--
James C. Stegen
NSF Postdoctoral Fellow in Bioinformatics
University of North Carolina
Chapel Hill, NC
919-962-8795
ste...@email.unc.edu
Hi, does anyone know if there is a way to extract the variance of each
variable explained in a structural equation model when using the sem()
function?
Thanks,
James Stegen
--
James C. Stegen
NSF Postdoctoral Fellow in Bioinformatics
University of North Carolina
Chapel Hill, NC
919-962-8795
Hello,
Some code using RGoogleDocs, which had been working smoothly since the
summer, just stopped working. I know that it worked on November 3rd, but it
doesn't work today. I've confirmed that the login and password still work
when I log in manually. I've confirmed that the URL gives the same
Hello Jim, hello all,
Thanks very much for the inputs, I used the code and it solved my
problem
special thanks to Jim Holtman for the code.
Regards,
karthick
--
View this message in context:
http://r.789695.n4.nabble.com/Parsing-txt-file-tp3035749p3036284.html
Sent from the R help
Hi R users
I have a very simple function. As the return of that function I want the
answer with 5 decimal places, but it hasn't worked with sprintf, nor with
format, nor print. This is how I used sprintf:
cuant-function(r,n){
d-seq(-1,1,by=0.001)
.
.
.
(SOME CALCULATIONS)
.
.
.
On 10/11/2010 10:51 AM, Liliana Pacheco wrote:
Hi R users
I have a very simple function. As the return of that function I want the
answer with 5 decimal places, but it hasn't worked with sprintf, nor with
format, nor print. This is how I used sprintf:
cuant-function(r,n){
d-seq(-1,1,by=0.001)
.
On Nov 10, 2010, at 9:51 AM, Liliana Pacheco wrote:
Hi R users
I have a very simple function. As the return of that function I want the
answer with 5 decimal places, but it hasn't worked with sprintf, nor with
format, nor print. This is how I used sprintf:
cuant-function(r,n){
I am getting the following error when using odfWeave
Error in xmlEventParse(infile, handlers = handlers, trim = FALSE, state =
state) :
File content_1.xml does not exist
This appears to be the same issue detailed in
Xiaobo,
You indeed need external 'zip' and 'unzip' utlities in the path, citing from
ffsave's help: using an external zip utility, e.g. for windows in Rtools on
[http://www.murdoch-sutherland.com/Rtools/];.
Please note that the mentioned utilities have a 4 GB limit for the zip file,
AFAIK. I
Dear James,
Model-implied variances and covariances of the non-error variables (observed
and latent) are given by (I - A)^-1 P [(I - A)^-1]', where the A and P
matrices are from the RAM formulation of the model, and are in the object
returned by sem(), and I is an identity matrix. Your model
Hello List,
I'm trying to find a convenient way of performing a weighted convolution on
multiple n-dimensional kernel density estimates (n-d pdfs) produced by the
np package function npudens()/npudist(). Are there any functions in R that
will take a list of functions and a list of weights and
Dear R community members:
The R Cookbook from O'Reilly is now available on-line in a Rough Cuts
version.
http://oreilly.com/catalog/9780596809164/
The book is a collection of recipes for R covering a variety of topics,
including: getting started, data structures, input and output,
Hi Mark,
Thanks for your response. After some detective work I figured out the answer
to my question.
The models
lm.split=lm(Y~B*V+B*N+V*N)
lmer.split=lmer(Y~V+N+V:N+(1|B)+(1|B:V)+(1|B:N))
contain exactly the same terms. The difference is that blocking factor (B)
is fixed in first model
Dear all
I'm getting a strange error when trying to use rtags() to generate
tags in Emacs format.
r...@liv-laptop:/usr/lib/R# R CMD rtags -o /usr/local/build/ETAGS
--no-Rd --no-c library/
Tagging R files under library/; writing to /usr/local/build/ETAGS
(overwriting)...
[..]
Error: '\.' is an
Hi.. I am using R as a user, another group built everything I am using. But I
need to install some packages and my only internet access is via a web-proxy.
I have hunted though the docs and I'm not finding anything on setting a proxy
in the package install command. Any help would be really
Hi Sue,
Try to download the packages through here
http://cran.r-project.org/web/packages/
http://cran.r-project.org/web/packages/then r/packages/install packages
from local zip files
if you are running a script do not forget to cal them
i hope that helps,
N
On Wed, Nov 10, 2010 at 12:26 PM,
R community,
I am creating a bivariate return level plot by adding calculated return
period values as lines onto an existing plot using the following code with
the points representing the return periods.
plot(H2,D2,pch=+,axes=TRUE)
points(H.10,D.10, type=l,col=blue)
points(H.20,D.20,
http://r.789695.n4.nabble.com/file/n3036574/plot.jpeg
--
View this message in context:
http://r.789695.n4.nabble.com/Plot-Axes-tp3036571p3036574.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Greetings Supreme Council of R Masters,
Like toddler, I have gotten my head stuck in the banisters of R ... again.
Let it be know I am still a neophyte in the R-community forum world, so
please don't flame me too bad.
I have two sets of data, each with a set of timestamps. I would like to
A few years back, I wrote some code to plot maps with the maptools
package. Now I am trying to reproduce my results, only to find out
that maptools has been updated and my code no longer works.
I've been able to fix the first part of it by forcing
spb -
Hello,
I have a basic question. Sorry if it is so evident
I have the following data file :
http://ekumen.homelinux.net/mydata.txt
I need to model Y~X-1 (simple linear regression through the origin) with
these data :
load(file=mydata.txt)
X=k[,1]
Y=k[,2]
aa=lm(Y~X-1)
dev.new()
On Wed, Nov 10, 2010 at 12:57 PM, Ian Craig ian.jh...@gmail.com wrote:
Greetings Supreme Council of R Masters,
Nice. :)
I have two sets of data, each with a set of timestamps. I would like to
somehow merge the datasets based on the timestamps and an individual
identifier. That is there are
apply the xlim/ylim in the initial plot
plot(..., xlim=range(H2, H.10, H.20, H.50, H.100), ylim=range(D2,
D.10, D.20, D.50, D.100))
On Wed, Nov 10, 2010 at 12:50 PM, dpender d.pen...@civil.gla.ac.uk wrote:
R community,
I am creating a bivariate return level plot by adding calculated return
Let it be know I am still a neophyte in the R-community forum world, so
please don't flame me too bad.
I have two sets of data, each with a set of timestamps. I would like to
somehow merge the datasets based on the timestamps and an individual
identifier. That is there are several
I would like to build a model in R to simulate the seed dispersal by one plant.
The plant produced 5 seeds and the probability of falling inside the eight
closest space was 0.8 and in the next space 0.2 and in the rest space 0:
0
0
0
0
0
0
0.2
0.2
0.2
0.2
0.2
0
0.2
0.8
0.8
0.8
0.2
Hi!
I have a set of 49 pairwise comparisons that I have done. From this I
would like to plot either histograms or the density plots of the values
I get. Now, I can plot one histogram per comparison, but I have problems
getting the output I want. When plotting like I normally would do:
Hi all,
I'm trying to apply loess regression to my data and then use the fitted
model to get the *standardized/studentized residuals. I understood that for
linear regression (lm) there are functions to do that:*
*
*
fit1 = lm(y~x)
stdres.fit1 = rstandard(fit1)
studres.fit1 = rstudent(fit1)
I was
My problem is that I have a data set for every day of measurement in a
seperate file and I want to plot one parameter of the data for all the days
in one graph. I tried to use for loop but only the last data remains in the
program memory, I don`t know how to plot each day`s data continusly after
Hi, does anyone know if there is a way to easily extract p-values from the
regsubsets() function?
Thanks,
James Stegen
p.s. this is a reposting due to me not putting in a useful subject heading...
--
James C. Stegen
NSF Postdoctoral Fellow in Bioinformatics
University of North Carolina
Chapel
Hi Harlan
I just tried to connect to Google Docs and I had ostensibly the same problem.
However, the password was actually different from what I had specified.
After resetting it with GoogleDocs, the getGoogleDocsConnection() worked
fine. So I don't doubt that the login and password are
Am Mittwoch, 10. November 2010, 19:22:38 schrieb Nasrin Pak:
My problem is that I have a data set for every day of measurement in a
seperate file and I want to plot one parameter of the data for all the days
in one graph. I tried to use for loop but only the last data remains in the
program
Hi all,
I've read the emails of Dan, Deepayan and Sundar about adding error bars to
the lattice plots (
https://stat.ethz.ch/pipermail/r-help/2006-October/114883.html), but I still
have the problem when I want to adding error bars to barchart. I tried both
the solution of Deepayan and Sundar but
On 10/11/2010 12:13 PM, Liviu Andronic wrote:
Dear all
I'm getting a strange error when trying to use rtags() to generate
tags in Emacs format.
r...@liv-laptop:/usr/lib/R# R CMD rtags -o /usr/local/build/ETAGS
--no-Rd --no-c library/
Tagging R files under library/; writing to
On 10/11/2010 12:26 PM, Caskey, Susan wrote:
Hi.. I am using R as a user, another group built everything I am using. But I
need to install some packages and my only internet access is via a web-proxy.
I have hunted though the docs and I'm not finding anything on setting a proxy
in the
I know the default plot() method expects dat$x and dat$y to exist and be
equal if you simply try
plot(dat)
Try looking at
methods(plot)
to see if there is a specified method from maptools,
then check documentation on ?plot.methodfrommaptools
Otherwise, check spb for $x and $y objects.
Several weeks ago my father told me one funny or not thing. I suggested him a
tool for resolving a trouble. To my great surprise it helped him and father
thanked me a lot and presented a gift. I suppose that this tool might be
good solution for solving varied troubles with excel files -
If you want to read in all the files and then set the range so you can
print a parameter from each one on a single chart, there is some
information in the archives about how to do this. A brief outline is
below (definitely untested)
allFiles - lapply(fileList, read.table, header=TRUE, ... other
Karin -
An example would have been nice.
Perhaps this will be helpful:
somedat = data.frame(x=sample(1:7,1000,replace=TRUE),
+ y=sample(1:7,1000,replace=TRUE),
+ z=rnorm(1000))
somedat$grp = interaction(somedat$x,somedat$y)
Now we need to
Hi,
we have a function that calculates a value for each day of the year.
For example an iterative function like
actvalue_0 := 0
actvalue_ {t+1} := actvalue_t + f(param_1,param_2,actvalue_t)
and for some days we have measurements. Now we want to choose param_1
and param_2 so that the function
Hi:
I have some old scripts from when I used to use Systat and have the
NRM(value1,value2) that I
need to use with R. Does anyone know the R equivalency to this function? Thanks
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish Wildlife Service
California,
So here is a sample of the datasets. I have modified them for the purposes
of circulating to solve this problem. I have already added the ID and key
as seen in the code below. Any ideas?
gpsdata
gpsARC Protocol TrackUTCDate UTCTime LocalDate LocalTimey
x Altitude_m
Aleksandr -
What happens when you use
library(maptools)
spb = readShapePoly('/home/sasha/Documents/maps/spb.shp')
- Phil Spector
Statistical Computing Facility
Department
Hi Ian,
Did you see the second part of my request, for a simple
reproducible example? At the very least, we need your
merge and sort code, since it sounds like that's where
the problem truly lies.
Sarah
On Wed, Nov 10, 2010 at 2:02 PM, Ian Craig ian.jh...@gmail.com wrote:
So here is a sample
Hi Oliver,
As a warning, I may be missing something too. I did not see something
explicit in base R or MASS. In a quick scan of the fourth edition of
the MASS book, I did not read anything that it is
illogical/unreasonable to try to find standardized residuals (but my
knowledge of local
P-values have little meaning in this context. Nor do regression coefficient
estimates and especially standard errors (they are biased).
Frank
-
Frank Harrell
Department of Biostatistics, Vanderbilt University
--
View this message in context:
Hi all,
I defined the following
#
myhist=function(x){
hist(x,xlab=,main=)
h=hist(x)
xfit=seq(min(x),max(x),length=100)
yfit=dnorm(xfit,mean(x),sd=sd(x))
yfit=yfit*diff(h$mids[1:2])*length(x)
lines(xfit, yfit, col=blue,
What hypothesis do you expect the p-values to be testing?
What you would get from regsubsets is unlikely to test any hypothesis of
interest.
If you really feel the need for p-values then the safest approach is probably
the function SnowsCorrectlySizedButOtherwiseUselessTestOfAnything in the
1 - 100 of 144 matches
Mail list logo
