[R] greek letters, text, and values in labels
Hello,
I want to have a title that will look something like:
"Results for \theta=2.1", given that I have a variable theta=2.1, and
\theta should show on the screen like the greek letter.
I've tried a lot of things:
theta <- 2.1
plot(1:10, main=expression(paste("Results for", theta, "=", eval(theta
or using bquote
plot(1:10, main=paste("Results for ", bquote(theta == .(theta
or using substitute, etc. I could not make it work. This should be easy.
I would appreciate your help.
Thanks,
Adrian
__
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 commented, minimal, self-contained, reproducible code.
[R] how to avoid infinite loops when sourcing R files
Hello,
I have a bunch of R files in a directory and I want to source all of them
with something like lapply(files, source).
I have a main.R file
source("C:/Temp/set.parms.R")
parms <- set.parms()
do.calcs(parms)
cat("I'm done with main.R\n")
Then I have set.parms.R function
set.parms <- function(){
cat("I'm in set.parms.\n"); flush.console()
directory <- "C:/Temp/"
files <- dir(directory, "\\.[rR]$", full.name=T)
files <- files[-grep("set.parms.R", files)] # remove infinite loop
lapply(files, source) # source them all
cat("Exiting set.parms.\n"); flush.console()
}
And other functions f1, f2, f3, etc. in the same directory that also
source set.parms.R. For example:
f1 <- function(){
source("H:/user/R/RMG/Energy/VaR/Overnight/Test/set.parms.R")
cat("I add two numbers.\n"); flush.console()
}
Because of the source command in f1, I get into an infinite loop. This
must be a common situation but I don't know how to avoid it.
I need the source(set.parms) in f1, f2, f3, etc. because I want use a
different combination of them in other projects.
Thanks,
Adrian
__
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 commented, minimal, self-contained, reproducible code.
Re: [R] learning lattice graphics
Check the documentation link from http://cm.bell-labs.com/cm/ms/departments/sia/project/trellis/software.writing.html Adrian On 5/26/07, Tyler Smith <[EMAIL PROTECTED]> wrote: > > Hi, > > I've just produced my first lattice plot - the graphic is very > impressive, but I only partly understand how it works. I crimped from > examples in MASS and the help pages to get most of what I want done. > Howver, I find both MASS and the help pages are a little too terse for > my needs as a beginner. Can anyone suggest web or print resources for > lattice graphics that provide a more thorough introduction? I've been > through the FAQ and standard manuals, but didn't see what I needed > there. > > Thanks, > > Tyler > > __ > 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 commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ 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 commented, minimal, self-contained, reproducible code.
[R] ROracle error in Windows. Memory could not be read.
I've seen from earlier posts that other people had problems installing
ROracle under Windows. I run R-2.3.1.
I got the Windows binaries for ROracle from
http://stat.bell-labs.com/RS-DBI/download/index.html
Here is my session:
> require(ROracle)
Loading required package: ROracle
Loading required package: DBI
[1] TRUE
> drv <- dbDriver("Oracle")
> drv
> con <- dbConnect(drv, user="USER", password="PASS", dbname="TEST")
A window pops up, with the message: The instruction at "0x6260c621"
referenced memory at "0x4a280ae2". The memory could not be read. And R
freezes.
Something is working, because I get to establish an Oracle driver.
Any ideas would be much appreciated.
Thank you,
Adrian
__
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 commented, minimal, self-contained, reproducible code.
Re: [R] Suppress blanks/spaces in character
> c <- "abcdef ghi"
> gsub(" +"," ", c)
[1] "abc def ghi"
>
On Fri, 27 Oct 2006 [EMAIL PROTECTED] wrote:
> Hi all
>
> I'm have a character vector and would like to suppress the blanks if there
> are more than one after the other.
>
> Example:
>
> Character value is: "abc def ghi"
> The result should be: "abc def ghi"
>
> I know that it's possible to delete the leading blanks with the command
> "trim". But how can I delete blanks within a character?
>
> Thank you very much in advance for an answer
> Thomas
>
> __
> 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 commented, minimal, self-contained, reproducible code.
>
__
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 commented, minimal, self-contained, reproducible code.
[R] Installing package rpvm under Windows
Hello, I'm trying to install the rpvm package under Windows, but I am having problems. I have pvm3.4 installed properly. I've defined the system variables PVM_ROOT = C:\PROGRA~1\pvm3.4\ PVM_ARCH = win32 When I try to install, I get this: C:\R\Packages>Rcmd INSTALL rpvm_1.0.1.tar.gz -- Making package rpvm ** WARNING: this package has a configure script It probably needs manual configuration ** adding build stamp to DESCRIPTION making DLL ... gcc -shared -s -o rpvm.dll rpvm.def rpvm_core.o rpvm_ser.o utils.orpvm_res.o -Lc:/PROGRA~1/R/R-24~1.0/bin -lR rpvm_core.o:rpvm_core.c:(.text+0x44): undefined reference to `pvm_perror' rpvm_core.o:rpvm_core.c:(.text+0x49): undefined reference to `pvm_exit' rpvm_core.o:rpvm_core.c:(.text+0x197): undefined reference to `pvm_mytid' rpvm_core.o:rpvm_core.c:(.text+0x1c7): undefined reference to `pvm_parent' rpvm_core.o:rpvm_core.c:(.text+0x207): undefined reference to `pvm_exit' rpvm_core.o:rpvm_core.c:(.text+0x2a8): undefined reference to `pvm_pstat' rpvm_core.o:rpvm_core.c:(.text+0x352): undefined reference to `pvm_kill' rpvm_core.o:rpvm_core.c:(.text+0x3ba): undefined reference to `pvm_tasks' rpvm_core.o:rpvm_core.c:(.text+0x6bf): undefined reference to `pvm_spawn' rpvm_core.o:rpvm_core.c:(.text+0x7b7): undefined reference to `pvm_initsend' rpvm_core.o:rpvm_core.c:(.text+0x7f7): undefined reference to `pvm_mkbuf' rpvm_core.o:rpvm_core.c:(.text+0x837): undefined reference to `pvm_freebuf' rpvm_core.o:rpvm_core.c:(.text+0x867): undefined reference to `pvm_getsbuf' rpvm_core.o:rpvm_core.c:(.text+0x897): undefined reference to `pvm_getrbuf' rpvm_core.o:rpvm_core.c:(.text+0x8d7): undefined reference to `pvm_setsbuf' rpvm_core.o:rpvm_core.c:(.text+0x917): undefined reference to `pvm_setrbuf' rpvm_core.o:rpvm_core.c:(.text+0x97d): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0x9f2): undefined reference to `pvm_pkdouble' rpvm_core.o:rpvm_core.c:(.text+0xa43): undefined reference to `pvm_pkstr' rpvm_core.o:rpvm_core.c:(.text+0xa9e): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xad4): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xb2e): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xb5f): undefined reference to `pvm_pkdouble' rpvm_core.o:rpvm_core.c:(.text+0xbc1): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xbf8): undefined reference to `pvm_pkstr' rpvm_core.o:rpvm_core.c:(.text+0xcc5): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xcf5): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xdc5): undefined reference to `pvm_pkint' rpvm_core.o:rpvm_core.c:(.text+0xdf5): undefined reference to `pvm_pkdouble' rpvm_core.o:rpvm_core.c:(.text+0xe6e): undefined reference to `pvm_send' rpvm_core.o:rpvm_core.c:(.text+0xee4): undefined reference to `pvm_mcast' rpvm_core.o:rpvm_core.c:(.text+0xf43): undefined reference to `pvm_recv' rpvm_core.o:rpvm_core.c:(.text+0xf9e): undefined reference to `pvm_nrecv' rpvm_core.o:rpvm_core.c:(.text+0xffe): undefined reference to `pvm_probe' rpvm_core.o:rpvm_core.c:(.text+0x107c): undefined reference to `pvm_trecv' rpvm_core.o:rpvm_core.c:(.text+0x113d): undefined reference to `pvm_bufinfo' rpvm_core.o:rpvm_core.c:(.text+0x1217): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x12ac): undefined reference to `pvm_upkdouble' rpvm_core.o:rpvm_core.c:(.text+0x130b): undefined reference to `pvm_upkstr' rpvm_core.o:rpvm_core.c:(.text+0x1374): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x13b4): undefined reference to `pvm_upkstr' rpvm_core.o:rpvm_core.c:(.text+0x1421): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x146f): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x14c1): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x150f): undefined reference to `pvm_upkdouble' rpvm_core.o:rpvm_core.c:(.text+0x1561): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x15be): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x1611): undefined reference to `pvm_upkint' rpvm_core.o:rpvm_core.c:(.text+0x166e): undefined reference to `pvm_upkdouble' rpvm_core.o:rpvm_core.c:(.text+0x16c0): undefined reference to `pvm_config' rpvm_core.o:rpvm_core.c:(.text+0x1894): undefined reference to `pvm_start_pvmd' rpvm_core.o:rpvm_core.c:(.text+0x191c): undefined reference to `pvm_addhosts' rpvm_core.o:rpvm_core.c:(.text+0x19c7): undefined reference to `pvm_delhosts' rpvm_core.o:rpvm_core.c:(.text+0x1a27): undefined reference to `pvm_halt' rpvm_core.o:rpvm_core.c:(.text+0x1acc): undefined reference to `pvm_mstat' rpvm_core.o:rpvm_core.c:(.text+0x1b53): undefined reference to `pvm_joingroup' rpvm_core.o:rpvm_core.c:(.text+0x1ba3): undefined reference to `pvm_lvgro
[R] how to clear the last error?
Hello, I would like to clear the last error that I get with geterrmessage(). Not even rm(list=ls()) clears it. Can I set it to NULL somehow? Thank you, Adrian __ 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
[R] weird behavior of nsmall in format
>From the help page of format, nsmall should control the number of digits. > format(0.123456789, nsmall = 10) [1] "0.1234567890" > format(0.123456789, nsmall = 1) [1] "0.1234568" > format(0.123456789, nsmall = 2) [1] "0.1234568" > format(0.123456789, nsmall = 8) [1] "0.12345679" It adds zeros fine but for format(0.123456789, nsmall = 1) I want the result to be 0.1. I want to format numbers with a fixed number of digits. A combination of round and format + nsmall does the job. sprintf will do it too, but I wondered if the current implementation of nsmall in format is correct. Thanks, Adrian __ 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
Re: [R] the large dataset problem
Take a look at the package filehash. It allows you to work with large
objects in R (bigger than your RAM) by storing them on the disk. The
objects are represented as pointers in R and have a small footprint in
memory. You can load all of them in an environment and access them with the
$ operator. I think filehash is more general than R.huge. R.huge works very
well with numerical 2D data only.
Adrian Dragulescu
On 7/31/07, Eric Doviak <[EMAIL PROTECTED]> wrote:
>
>
> Just a note of thanks for all the help I have received. I haven't gotten a
> chance to implement any of your suggestions because I'm still trying to
> catalog all of them! Thank you so much!
>
> Just to recap (for my own benefit and to create a summary for others):
>
> Bruce Bernzweig suggested using the R.huge package.
>
> Ben Bolker pointed out that my original message wasn't clear and asked
> what I want to do with the data. At this point, just getting a dataset
> loaded would be wonderful, so I'm trying to trim variables (and if possible,
> I would also like to trim observations). He also provided an example of
> "vectorizing."
>
> Ted Harding suggested that I use AWK to process the data and provided the
> necessary code. He also tested his code on older hardware running GNU-Linux
> (or Unix?) and showed that AWK can process the data even when the computer
> has very little memory and processing power. Jim Holtman had similar success
> when he used Cygwin's UNIX utilities on a machine running MS Windows. They
> both used the following code:
>
> gawk 'BEGIN{FS=","}{print $(1) "," $(1000) "," $(1275) "," $(5678)}'
> < tempxx.txt > newdata.csv
>
> Fortunately, there is a version of GAWK for MS Windows. ... Not that I
> like MS Windows. It's just that I'm forced to use that 19th century
> operating system on the job. (After using Debian at home and happily running
> RKWard for my dissertation, returning to Windows World is downright
> depressing).
>
> Roland Rau suggested that I use a database with RSQLite and pointed out
> that RODBC can work with MS Access. He also pointed me to a sub-chapter in
> Venables and Ripley's _S Programming_ and "The Whole-Object View" pages in
> John Chamber's _Programming with Data_.
>
> Greg Snow recommended biglm for regression analysis with data that is
> too large to fit into memory.
>
> Last, but not least, Peter Dalgaard pointed out that there are options
> within R. He suggests using the colClasses= argument for when "reading" data
> and the what= argument for "scanning" data, so that you don't load more
> columns than necessary. He also provided the following script:
>
> dict <- readLines("
> ftp://www.sipp.census.gov/pub/sipp/2004/l04puw1d.txt";)
> D.lines <- grep("^D ", dict)
> vdict <- read.table(con <- textConnection(dict[D.lines])); close(con)
> head(vdict)
>
> I'll try these solutions and report back on my success.
>
> Thanks again!
> - Eric
>
> __
> 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 commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
__
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 commented, minimal, self-contained, reproducible code.
[R] autofill for a tkentry field
Hello all,
I'm experimenting with tlk/tk and I find it very useful. I would like to
have the text of a tkentry be automatically filled when the user starts
typing, according to a predetermined list. I've found a solution on the
tcl wiki pages, but I don't quite know how to translate it in R... Here is
the tcl code: http://wiki.tcl.tk/13267
And here is my code that does not work.
autofill <- function(){
cat("I was here!")
entry <- tclvalue(f1)
ind <- grep(entry, fruitList)
if (length(ind)>0){f1 <- tclVar(fruitList[ind[1]])}
}
fruitList <<- c("Apple", "Banana", "Peach", "Pear", "Plum", "Mango")
tkdestroy(tt)
tt <- tktoplevel()
pw <- tkframe(tt)
ppw <- tk2labelframe(pw, text="Fruits:")
f1 <<- tclVar("")
entry1.f <- tkentry(ppw, width=30, textvariable=f1, validate="key",
validatecommand=autofill)
tkgrid(entry1.f)
tkgrid(ppw)
tkgrid(pw)
I don't know how to pass arguments to the autofill function, so I can
manipulate the textvariable. Any help is much appreciated. If we get a
working solution, maybe we should add this example to
http://www.sciviews.org/_rgui/tcltk/index.html.
Thank you,
Adrian
__
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 commented, minimal, self-contained, reproducible code.
[R] R in a production environment
Recently there have been a lot of discussions on the list on the use of R in industry and comparisons with SAS & SPSS. These are very important questions and it might be worth having them all grouped together for reference. I would like to add to these topics a related one. Can we share our experience of using R in a production environment? I work in a financial company that is cautiously interested in using R in a production environment. I'm trying to convince them that R is serious software, modern and powerful. Their first question is who else out there is using R. I can guess from some contributed packages, or from the list postings, but it would be great if there would be a centralized place. We could have a standard set of questions and a contact person (maybe on the wiki site?) What do you think? I believe it could create a snowball effect to increase the visibility and use of R beyond academia and personal use. Thanks, Adrian __ 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 commented, minimal, self-contained, reproducible code.
[R] dumping the summary of lm to a text file
Hello list,
I have a linear regression
ctl<- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt<- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2,10,20, labels=c("Ctl","Trt"))
weight <- c(ctl, trt)
reg<- lm(weight ~ group)
sreg <- summary(reg)
and I would like to dump exactly what I see on the console with
print(sreg) to a text file. I've tried using cat but it did not work.
I've read the list questions but I did not see a solution.
Can you help with this?
Thank you,
Adrian Dragulescu
__
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
[R] problem with DCOM server in R 2.1.0
Hello list, I just installed R 2.1.0. When I try to run the "Server 01 - Basic test" I get the error "Fatal error: unable to open the base package". Any solution? Thank you, Adrian Dragulescu __ 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
[R] profile error on an nls object
Hello all, This is the error message that I get. > hyp.res <- nls(log(y)~log(pdf.hyperb(theta,X)), data=dataModel, + start=list(theta=thetaE0), + trace=TRUE) 45.54325 : 0.100 1.3862944 -4.5577142 0.0005503 3.728302 : 0.0583857346 0.4757772859 -4.9156128701 0.0005563154 1.584317 : 0.0194149477 0.3444648833 -4.9365149150 0.0004105426 1.569333 : 0.0139310639 0.3824648048 -4.9024001228 0.0004089738 1.569311 : 0.0137155342 0.3888648619 -4.8979817546 0.0004137501 1.569311 : 0.0136895846 0.3893564152 -4.8976182201 0.0004141057 1.569311 : 0.0136876315 0.3894059947 -4.8975821760 0.0004141343 > hyp.res.S <- summary(hyp.res) > hyp.res.S Formula: log(y) ~ log(pdf.hyperb(theta, X)) Parameters: Estimate Std. Error t value Pr(>|t|) theta1 0.0136876 0.0359964 0.3800.705 theta2 0.3894060 0.3079860 1.2640.211 theta3 -4.8975822 0.2219928 -22.062 <2e-16 *** theta4 0.0004141 0.0005457 0.7590.451 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 Residual standard error: 0.1542 on 66 degrees of freedom Correlation of Parameter Estimates: theta1theta2theta3 theta2 -0.02168 theta3 -0.02029 0.997736 theta4 -0.97182 -0.008054 -0.008952 > pr1 <- profile(hyp.res) 1.825584 : 0.3894059947 -4.8975821760 0.0004141343 1.58426 : 0.373691474 -4.909091289 0.000824045 1.583673 : 0.4176596873 -4.8774106487 0.0008176545 1.583670 : 0.4196944963 -4.8760375504 0.0008187918 1.583670 : 0.4199010211 -4.8758854269 0.0008188162 1.624899 : 0.449756713 -4.854643555 0.001215014 1.624743 : 0.46804752 -4.84185838 0.00122343 1.624741 : 0.470384534 -4.840195293 0.001224199 1.624741 : 0.470638282 -4.840013199 0.001224298 1.624741 : 0.470670500 -4.839990112 0.001224309 1.692158 : 0.522188258 -4.803565745 0.001635778 1.691853 : 0.540794581 -4.791027785 0.001650730 1.691847 : 0.544973564 -4.788090229 0.001652321 1.691847 : 0.545500818 -4.787718964 0.001652616 1.691847 : 0.545592388 -4.787654441 0.001652658 1.784749 : 0.622277872 -4.734086833 0.002091090 1.784039 : 0.642139831 -4.721442413 0.00211 1.784022 : 0.649929188 -4.716068239 0.002119126 1.784021 : 0.651094995 -4.715267692 0.002119963 1.784021 : 0.65136956 -4.71507850 0.00212012 1.784021 : 0.651420684 -4.715043256 0.002120153 1.901667 : 0.760981513 -4.639871097 0.002604136 1.899765 : 0.782870773 -4.627113544 0.002644289 1.899703 : 0.798044378 -4.616913842 0.002653139 1.899699 : 0.800930030 -4.614996517 0.002655646 1.899699 : 0.801815115 -4.614404602 0.002656286 1.899699 : 0.802033012 -4.614258879 0.002656455 1.899699 : 0.802090888 -4.614220164 0.002656499 2.042311 : 0.960722487 -4.508069592 0.003221186 2.036028 : 0.985442669 -4.495703574 0.003291724 2.035767 : 1.017849320 -4.474692317 0.003317203 2.035736 : 1.026482531 -4.469203277 0.003326355 2.035733 : 1.029937702 -4.466987738 0.003329627 2.035732 : 1.031114541 -4.466233213 0.003330776 2.035732 : 1.03153247 -4.46596514 0.00333118 2.035732 : 1.031679019 -4.465871141 0.003331322 2.035732 : 1.031730150 -4.465838341 0.003331371 1.583425 : 3.595503e-01 -4.918824e+00 1.793668e-05 1.583270 : 3.783622e-01 -4.905413e+00 2.175254e-05 1.58327 : 3.793739e-01 -4.904683e+00 2.134353e-05 1.58327 : 3.794562e-01 -4.904622e+00 2.133204e-05 1.624847 : 0.3695765499 -4.9116128267 -0.0003687005 1.624748 : 0.3877816645 -4.8985581022 -0.0003699797 1.624747 : 0.3890277914 -4.8976645481 -0.0003706884 1.624747 : 0.3891837724 -4.8975508910 -0.0003707385 1.693266 : 0.3989123147 -4.8904787597 -0.0007628478 1.693170 : 0.4172349434 -4.8774484377 -0.0007692982 1.693168 : 0.4193992605 -4.8759045536 -0.0007705839 1.693168 : 0.4197576966 -4.8756463144 -0.0007707548 1.693168 : 0.4198094076 -4.8756090629 -0.0007707816 1.788041 : 0.450653198 -4.853510937 -0.001173674 1.787884 : 0.469850125 -4.840178495 -0.001186483 1.787878 : 0.473941725 -4.837285088 -0.001188994 1.787878 : 0.474783800 -4.836687703 -0.001189518 1.787878 : 0.474960093 -4.836562526 -0.001189627 1.787878 : 0.474996400 -4.836536741 -0.001189650 1.908362 : 0.531011661 -4.796878005 -0.001614805 1.907987 : 0.552282889 -4.782714884 -0.001637055 1.907965 : 0.560264832 -4.777154923 -0.001642437 1.907964 : 0.562346209 -4.775709348 -0.001644026 1.907963 : 0.56295255 -4.77528733 -0.00164447 1.907963 : 0.563123513 -4.775168321 -0.001644596 1.907963 : 0.563172014 -4.775134556 -0.001644632 2.053048 : 0.653571857 -4.712183501 -0.002111091 2.051874 : 0.678913471 -4.696423681 -0.002150202 2.051785 : 0.69544095 -4.68518302 -0.00216328 2.051772 : 0.701166962 -4.681323343 -0.002168564 2.051770 : 0.703471045 -4.679765756 -0.002170603 2.051769 : 0.704360670 -4.679164327 -0.002171397 2.051769 : 0.704707691 -4.678929683 -0.002171706 2.051769 : 0.704842651 -4.678838422 -0.002171826 2.051769 : 0.704895042 -4.678802995 -0.00217187
Re: [R] cumulative density function
Because the result is a list, you can get your number like this: integrate(dnorm,0,2)[[1]] You don't have an expression for the probability function? You can do a nonlinear regression to a distribution you belive in, and then integrate that. Or, you can calculate the experimental cdf by simple summation. Best, Adrian __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] memory problems with lm
Hello list, I've seen the recent discussions documenting problems with lm. I have encountered the following problem. I use WinXP Pro with service pack 1, and R 1.9.0, on a XEON 2GHz, with 1GB of RAM. > eff.fro std.dev mean NSTRDSP 7.403749e-01 1.215686e-01 CPFGEP 9.056763e+00 1.815686e+00 WSWOLF 4.703588e+05 1.112832e+05 NPILGRIM 1.017640e+06 2.134335e+05 WSNMILE 1.367312e+07 1.892021e+06 WSHIDESL 1.830811e+07 1.892021e+06 > reg <- lm(log(mean) ~ log(std.dev), data=eff.fro) Error in model.matrix.default(mt, mf, contrasts) : cannot allocate vector of length 1074790452 > log(eff.fro$mean) [1] -2.1072763 0.5964635 11.6198339 12.2710808 14.4531561 [6] 14.4531561 > reg <- lm(log(eff.fro$mean) ~ log(eff.fro$std.dev)) Error: cannot allocate vector of size 3360077 Kb > lef <- log(eff.fro) > lef std.dev mean NSTRDSP -0.3005986 -2.1072763 CPFGEP2.2035117 0.5964635 WSWOLF 13.0612512 11.6198339 NPILGRIM 13.8329973 12.2710808 WSNMILE 16.4309427 14.4531561 WSHIDESL 16.7228546 14.4531561 > lef <- log(eff.fro) > reg <- lm(lef$mean ~ lef$std.dev) Here the my computer completely crashed. A window poped-up and said memory problem at address ..., and if I want to debug. I ran the same code one more time, and it worked but it did not work how I wanted (where is the slope?): > reg <- lm(lef$mean ~ lef$std.dev) > reg Call: lm(formula = lef$mean ~ lef$std.dev) Coefficients: (Intercept) 8.548 > > summary(reg) Call: lm(formula = lef$mean ~ lef$std.dev) Residuals: 1 2 3 4 5 6 -10.655 -7.951 3.072 3.723 5.905 5.905 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept)8.548 2.9992.85 0.0358 * --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 Residual standard error: 7.346 on 5 degrees of freedom I ran again: reg <- lm(log(mean) ~ log(std.dev), data=eff.fro) and I get the pop-up: The instruction at "0x6b4c45a5" referenced memory at "0x0032374a". The memory could not be "read". Click OK to terminate the program. Any ideas? Thank you, Adrian __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] memory problems with lm
If I enforce the variables to be numeric it works fine. > str(eff.fro) `data.frame': 6 obs. of 2 variables: $ std.dev: num [, 1:6] 7.40e-01 9.06e+00 4.70e+05 1.02e+06 1.37e+07 ... ..- attr(*, "dimnames")=List of 1 .. ..$ : chr "NSTRDSP" "CPFGEP" "WSWOLF" "NPILGRIM" ... $ mean : num 1.22e-01 1.82e+00 1.11e+05 2.13e+05 1.89e+06 ... > gc() used (Mb) gc trigger (Mb) Ncells 578941 15.51166886 31.2 Vcells 589444 4.52377385 18.2 > eff.fro std.dev mean NSTRDSP 7.403749e-01 1.215686e-01 CPFGEP 9.056763e+00 1.815686e+00 WSWOLF 4.703588e+05 1.112832e+05 NPILGRIM 1.017640e+06 2.134335e+05 WSNMILE 1.367312e+07 1.892021e+06 WSHIDESL 1.830811e+07 1.892021e+06 > reg <- lm(log(as.numeric(mean)) ~ log(as.numeric(std.dev)), data=eff.fro) > reg Call: lm(formula = log(as.numeric(mean)) ~ log(as.numeric(std.dev)), data = eff.fro) Coefficients: (Intercept) log(as.numeric(std.dev)) -1.63680.9864 Adrian On Thu, 29 Apr 2004, Liaw, Andy wrote: > Can you show us the output of str(eff.fro)? Do you have other things in the > global environment or the search path that's taking up memory? What does > gc() say? > > Andy > > > From: Adrian Dragulescu > > > > Hello list, > > > > I've seen the recent discussions documenting problems with lm. > > > > I have encountered the following problem. I use WinXP Pro with > > service pack 1, and R 1.9.0, on a XEON 2GHz, with 1GB of RAM. > > > > > eff.fro > > std.dev mean > > NSTRDSP 7.403749e-01 1.215686e-01 > > CPFGEP 9.056763e+00 1.815686e+00 > > WSWOLF 4.703588e+05 1.112832e+05 > > NPILGRIM 1.017640e+06 2.134335e+05 > > WSNMILE 1.367312e+07 1.892021e+06 > > WSHIDESL 1.830811e+07 1.892021e+06 > > > reg <- lm(log(mean) ~ log(std.dev), data=eff.fro) > > Error in model.matrix.default(mt, mf, contrasts) : > > cannot allocate vector of length 1074790452 > > > log(eff.fro$mean) > > [1] -2.1072763 0.5964635 11.6198339 12.2710808 14.4531561 > > [6] 14.4531561 > > > reg <- lm(log(eff.fro$mean) ~ log(eff.fro$std.dev)) > > Error: cannot allocate vector of size 3360077 Kb > > > lef <- log(eff.fro) > > > lef > > std.dev mean > > NSTRDSP -0.3005986 -2.1072763 > > CPFGEP2.2035117 0.5964635 > > WSWOLF 13.0612512 11.6198339 > > NPILGRIM 13.8329973 12.2710808 > > WSNMILE 16.4309427 14.4531561 > > WSHIDESL 16.7228546 14.4531561 > > > lef <- log(eff.fro) > > > reg <- lm(lef$mean ~ lef$std.dev) > > > > Here the my computer completely crashed. A window poped-up and said > > memory problem at address ..., and if I want to debug. > > > > I ran the same code one more time, and it worked but it did not work > > how I wanted (where is the slope?): > > > reg <- lm(lef$mean ~ lef$std.dev) > > > reg > > > > Call: > > lm(formula = lef$mean ~ lef$std.dev) > > > > Coefficients: > > (Intercept) > > 8.548 > > > > > > > > summary(reg) > > > > Call: > > lm(formula = lef$mean ~ lef$std.dev) > > > > Residuals: > > 1 2 3 4 5 6 > > -10.655 -7.951 3.072 3.723 5.905 5.905 > > > > Coefficients: > > Estimate Std. Error t value Pr(>|t|) > > (Intercept)8.548 2.9992.85 0.0358 * > > --- > > Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 > > > > Residual standard error: 7.346 on 5 degrees of freedom > > > > I ran again: > > > > reg <- lm(log(mean) ~ log(std.dev), data=eff.fro) > > > > and I get the pop-up: > > The instruction at "0x6b4c45a5" referenced memory at "0x0032374a". > > The memory could not be "read". Click OK to terminate the program. > > > > > > Any ideas? Thank you, > > Adrian > > > > __ > > [EMAIL PROTECTED] mailing list > > https://www.stat.math.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide! > > http://www.R-project.org/posting-guide.html > > > > > > > -- > Notice: This e-mail message, together with any attachments, contains > information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New > Jersey, USA 08889), and/or its affiliates (which may be known outside the > United States as Merck Frosst, Merck Sharp & Dohme or MSD and in Japan as > Banyu) that may be confidential, proprietary copyrighted and/or legally > privileged. It is intended solely for the use of the individual or entity > named on this message. If you are not the intended recipient, and have > received this message in error, please notify us immediately by reply e-mail > and then delete it from your system. > -- > __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Problem Connecting to Oracle with R from Windows XP
You also need the ROracle package.
On 8/15/07, Song, Alex <[EMAIL PROTECTED]> wrote:
>
> Hello,
>
>
>
> I installed RGui 2.5.1 and package DBI on Windows XP and tried to connect
> to Oracle database which is on a Linux server. When I tried to use
> dbDriver("Oracle"), I got an error as follows:
>
>
>
> > drv <- dbDriver("Oracle")
>
> Error in do.call(as.character(drvName), list(...)) :
>
> could not find function "Oracle"
>
> >
>
>
>
> Could anyone tell me how to connect to Oracle with R from Windows XP? Do I
> need to configure any environment variables? Do I need to configure ODBC? Do
> I need to install any other packages? I have Oracle client 10g installed on
> my computer and I can connect to the Oracle database using other client
> software like Toad, SQLPlus, or SQL Developer.
>
>
>
> Thank you very much.
>
>
>
> Alex
>
>
>
> Alex Song
>
> Data Management Specialist / Spécialiste en gestion des données
>
> National Forest Inventory / Inventaire forestier national
>
> Pacific Forestry Centre / Centre de foresterie du Pacifique
> Canadian Forest Service / Service canadien des forêts
> Natural Resources Canada / Ressources naturelles Canada
>
> Government of Canada / Gouvernement du Canada
>
> 506 West Burnside Road / 506 chemin Burnside ouest
>
> Victoria, BC V8Z 1M5
>
> Phone (250) 363-3342
>
> Facs: (250) 363-0775
>
>
>
> Email: [EMAIL PROTECTED]
>
>
>
>
> [[alternative HTML version deleted]]
>
>
> __
> 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 commented, minimal, self-contained, reproducible code.
>
>
[[alternative HTML version deleted]]
__
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 commented, minimal, self-contained, reproducible code.
[R] filehash package: error in unserialize(con)
This question is related to the filehash package only.
I would like to use the filehash package because of it's ability to deal
with large amounts of data. I've set up a filehash database on a network
drive and populated it with several GB of data. Somehow, I must have
corrupted the files because I cannot access the objects anymore. Here is
what I do:
> require(filehash)
> filehashOption(defaultType = "RDS") # need this, dbInit will fail
> db <- dbInit("//nas-msw-07/app/Prices/Historical")
> hPP <- db2env(db)# put all the data in an enviroment
> ("commod xyz ex" %in% ls(hPP))
[1] TRUE
> dbExists(db, "commod xyz ex")
[1] TRUE
> dbFetch(db, "commod xyz ex")
Error in unserialize(con) : error reading from connection
> hPP[["commod xyz ex"]]
Error in unserialize(con) : error reading from connection
>
How do I get to touch the data again? It's really slow to recreate the
filehash again, and it seems to be very fragile. Any suggestions much
appreciated.
Thank you,
Adrian Dragulescu
__
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 commented, minimal, self-contained, reproducible code.
