Displaying 20 results from an estimated 30000 matches similar to: "loop"
2007 Oct 19
2
header
Hi everyone, I read the following data
> setwd("d:/")
> mydata <- read.table("mydata.txt", header=T)
> attach(mydata)
> mydata
name X20_1_TT_Z_e4 X20_2_TT_Z_e4
1 A1 102 22
2 A2 103 21
3 A3 107 20
> name
[1] A1 A2 A3
Levels: A1 A2 A3
> mode(name)
[1] "numeric"
Why name is
2005 Mar 23
4
sampling from a mixture distribution
Dear R users,
I would like to sample from a mixture distribution p1*f(x1)+p2*f(x2). I
usually sample variates from both distributions and weight them with their
respective probabilities, but someone told me that was wrong. What is the
correct way?
Vumani
2019 Sep 11
4
Fw: Calling a LAPACK subroutine from R
Berend,
I do not think this works with gfortran 7+. I am calling the BLAS
subroutine dgemv from Fortran code in my package eha, and the check
(with R-devel) gives:
gmlfun.f:223:1: warning: type of ?dgemv? does not match original
declaration [-Wlto-type-mismatch]
& score, ione)
^
/home/gobr0002/R/src/R-devel/include/R_ext/BLAS.h:107:1: note: type
mismatch in parameter 12
2019 Sep 11
4
Fw: Calling a LAPACK subroutine from R
Sorry for cross-posting, but I realized my question might be more appropriate for r-devel...
Thank you,
Giovanni
________________________________________
From: R-help <r-help-bounces at r-project.org> on behalf of Giovanni Petris <gpetris at uark.edu>
Sent: Tuesday, September 10, 2019 16:44
To: r-help at r-project.org
Subject: [R] Calling a LAPACK subroutine from R
Hello R-helpers!
2019 Sep 11
1
Fw: Calling a LAPACK subroutine from R
On 2019-09-11 22:16, Avraham Adler wrote:
> Can you write a small C function that calls LAPACK call that fro your
> Fortran code? Yes, an extra step but maybe less traumatic than rewriting
> parts of LAPACK directly.
Yes, I know how to do that, but I find it somewhat bizarre that it is
impossible to call a Fortran subroutine from Fortran. And rewriting
'dgemv' was simple:
2002 May 14
2
tapply and empty subsets
I am using tapply to compute means by group, as in
> tapply(y, z, mean)
1 2 3
21.00757 20.50031 NA
where y is a vector and z a factor with levels 1:3.
I would like to get a zero instead of the NA I get when a subset is
empty. The problem is that FUN is not applied to empty cells of the
ragged array. (Of course I am willing to define a function my.mean
that returns
2010 Aug 18
3
libreadline problem
Hello,
Recently updated Ubuntu to 10.04, I have installed the r-base and
r-base-dev packages from CRAN with
sudo apt-get install r-base
But now when I want to start R I get the following error:
gpetris at definetti:~$ R
/usr/lib/R/bin/exec/R: symbol lookup
error: /usr/local/lib/libreadline.so.6: undefined symbol: PC
Any clue??? Has anybody seen that before? I have tried to reinstall the
2003 Oct 29
1
restarting split.screen
Is there a way of `restarting' split.screen?
This is what I am getting:
> close.screen()
[1] 10 11 12 13
> close.screen(all=TRUE)
Error in par(args) : parameter "i" in "mfg" is out of range
> graphics.off()
> x11()
> close.screen()
[1] 10 11 12 13
> close.screen(all=TRUE)
Error in par(args) : parameter "i" in "mfg" is out of range
As
2007 Jan 18
4
Reading contingency tables
I am trying to read an ftable using read.ftable, but I get the
following error message:
> jobSatTable <- read.ftable("http://definetti.uark.edu/~gpetris/stat5333/jobSatisfaction.dat",skip=2)
Error in seek(file, where = 0) : no applicable method for "seek"
In addition: Warning messages:
1: no non-missing arguments to max; returning -Inf
2: no non-missing arguments to
2004 May 12
4
non-interactive call to R (running an R package as a stand-alone application)
Is there a way I can have R automatically execute the commands in a source file without ever having to use R interactively? If so, what arguments should I pass to the UNIX call to R? I need to do this to run several R jobs in parallel.
An alternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on
2004 Jan 20
4
matrix exponential: M^0
> -----Original Message-----
> From: Federico Calboli
> Sent: Tuesday, January 20, 2004 5:40 PM
> To: r-help
> Subject: [R] matrix exponential: M^0
>
> I would like to ask why the zeroeth power of a matrix gives me a matrix
> of ones rather than the identity matrix:
>
> > D<-rbind(c(0,0,0),c(0,0,0),c(0,0,0))
> > D<-as.matrix(D)
> > D
>
2008 Oct 09
2
Two math expressions in plot
Hello!
I am trying to put two math expressions in the title of a plot. As
you can see below, I can place correctly one expression at a time, but
not both. Ideally I would like to have them separated by a comma. Any
suggestions?
> k <- 1
> n.eff <- c(20, 30)
> ### this works
> plot(0,0, main = substitute(n == k, list(k = k)))
> ### this works
> plot(0,0, main =
2012 May 25
1
knitr customization
I am trying to transition from Sweave to knitr, but there are a few
things about customization of the appearence of R input and output that
I did not get yet. Maybe somebody on the list can help me.
In my Sweave presentations I used a slanted font for the R input and a
normal font for the output, both in a small font. I also indented
everything by an extra 2em. All this was achieved by the
2019 Sep 12
1
Calling a LAPACK subroutine from R
Followup:
I have checked my package nleqslv which uses dgemv only from Fortran, on Kubuntu 18.04 with the development version of R.
No errors or problems.
Berend
> On 12 Sep 2019, at 08:57, Berend Hasselman <bhh at xs4all.nl> wrote:
>
>
> I have tried what I proposed in a virtual Kubuntu 18.04 which uses gfortran 7.4.
> I used the latest development version of R.
>
2011 Mar 09
4
Help with read.csv
Hello,
I have a file that looks like this:
Date,Hour,DA_DMD,DMD,DA_RTP,RTP,,
1/1/2006,1,3393.9,3412,76.65,105.04,,
1/1/2006,2,3173.3,3202,69.20,67.67,,
1/1/2006,3,3040.0,3051,69.20,77.67,,
1/1/2006,4,2998.2,2979,67.32,69.10,,
1/1/2006,5,3005.8,2958,65.20,68.34,,
where the ',' is the separator and I tried to read it into R, but...
> y <- read.csv("Data/Data_tmp.csv",
2004 Feb 06
3
quantile function
I am trying to `cut' a continuous variable into contiguous classes
containing approximately an equal number of observations. I thought
quantile() was the appropriate function to use in order to find the
breakpoints, but I end up with classes of different sizes - see
example below. Does anybody have an explanation for that? And what is
the `recommended' way of computing what I am looking
2003 Oct 24
5
how to remove NaN columns ?
How can I remove columns with NaN entries ?
Here is my simple example:
> data <- read.csv("test.csv")
> xdata <- data[3:length(data)]
> xs <- lapply(xdata, function(x){(x - mean(x))/sqrt(var(x))})
> x <- data.frame(xs)
> x
C D E F
1 -0.7071068 NaN -0.7071068 -0.7071068
2 0.7071068 NaN 0.7071068 0.7071068
2012 May 09
1
Sweave, beamer and alert within code chunks
Hi all,
Using Beamer, in order to highlight a piece of R code I do something
like this - note the "\structure" and "\alert" commands:
\begin{semiverbatim}
> mleOut <- \structure{dlmMLE}(Nile,
+ parm = c(0.2, 120), # initial values for optimizer
+ lower = c(1e-7, 0)) \alert<2>{# V must be positive}
> mleOut$convergence
2005 Mar 22
1
Package vignette and build
Hello,
I am writing a package called 'DLM' containing a vignette.
The vignette contains a chunck with the function call 'library(DLM)'.
This worked fine with 'R CMD check DLM', but when it comes to building
the package with 'R CMD build DLM' I get the following error message:
* creating vignettes ... ERROR
Error: chunk 1
Error in library(DLM) : There is no
2010 Jan 21
2
What is the difference between S and R?
On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris <GPetris at uark.edu> wrote:
>
> I have found a good refernce to be "S Programming" by Venables and
> Ripley.
I'll take a look at this book.
Since S and R are not completely the same, there are delicate
differences between S and R, which an S book may be confusing if I
want to learn R. One example is Statistical Models