> On Wed, Sep 4, 2013 at 1:45 PM, Daniel Hornung
wrote:
> > Hello,
> >
> > can it be that xyplot does not support the lwd argument?
> >
> > At least here, the following still shows thin lines, as opposed to the
> > regular
> > plot command:
> >
> > xyplot(Sepal.Length ~ Sepal.Width, data = iris
It might be, but with appropriate indexes a SQL engine (via sqldf or RODBC for
example) might be able to do it that way anyway.
---
Jeff NewmillerThe . . Go Live...
DCN:Basics
On Sep 4, 2013, at 10:39 PM, Euna Jeong wrote:
Hi,
I have questions about R2 used in pls (or multivariate analysis).
Is R2 same with the square of the PCC (Pearson Correlation
Coefficient)?
I found the following description from wiki (Coefficient of
determination)
Dear Colleagues,
I'm working on a Delphi study comparing perceptions of high school
technology teachers and university engineering educators about the
importance of concepts about engineering for HS students to learn as part of
their fundamental education. I'm actually doing this as part of my
I have been using R version 2.15.1 happly along side R Commander. I then
tried to go through a fresh install using the latest version of R, R 3.0.1.
The trouble started when I wanted to install Rcmdr, I kept getting an error
message about previous installs of R Commander. I went through and
deinsta
Hi,
I have questions about R2 used in pls (or multivariate analysis).
Is R2 same with the square of the PCC (Pearson Correlation Coefficient)?
I found the following description from wiki (Coefficient of determination)
Similarly, in linear least squares regression with a
Hi Spencer,
It really is not very hard, and I have never had issue with it:
http://www.oracle.com/technetwork/java/javase/downloads/jdk7-downloads-1880260.html
Just download the x86 and x64 versions for your OS and install. Worst
case, you need to add the directory to the PATH variable in Windo
Sometimes one has to really read the manual carefully.
"If non-trivial bounds are supplied, this
method will be selected, with a warning." (re L-BFGS-B)
Several of us have noted problems occasionally with this code.
You might want to look at the box constrained codes offered in optimx
pa
On 9/4/2013 6:09 PM, Ista Zahn wrote:
Hi Spencer,
Why don't you want to install 64bit Java?
That may be a reasonable approach.
I may have Java confused with something else, but I remember
hearing that it was difficult or unwise to try to install both 32- and
64-bit versions of
Hi Spencer,
Why don't you want to install 64bit Java?
On Wed, Sep 4, 2013 at 6:12 PM, Spencer Graves
wrote:
> Hello, All:
>
>
> What do you recommend for reading a Google Spreadsheet into R? I didn't find
> anything useful using library(sos); findFn('google spreadsheet').
>
>
> I can solve the p
Thank you everyone for your help.
> Date: Wed, 4 Sep 2013 20:00:02 -0400
> From: murdoch.dun...@gmail.com
> To: dcarl...@tamu.edu
> CC: pmassico...@hotmail.com; ruipbarra...@sapo.pt; r-help@r-project.org
> Subject: Re: [R] Histogram
>
> On 13-09-04 4:44 PM, David Carlson wrote:> Good question. I
On 13-09-04 5:56 PM, Noah Silverman wrote:
Hi,
Working with R, I often want to copy and paste some values somewhere else.
(Its not worth saving a CSV file for a dozen or so entries.) Or, I may want to
copy all the names of an object into some code.
Besides the other suggestions, the data e
On 13-09-04 4:44 PM, David Carlson wrote:> Good question. It turns out
that the manual page does not tell
> the whole story.
Do you really think the manual page would be improved if it went into as
much detail as you give below? It does say clearly that breaks is a
"suggestion only". I don't
Hi Edouard,
In terms of speed, your new solution may not be that much different from the
old one:
#large matrix
M<- 10
N<- 1e3
set.seed(249)
A<- matrix(sample(1:10,1e5*4,replace=TRUE),1e5,4)
B<- matrix(NA,ncol=ncol(A),nrow=0)
system.time({
set.seed(54)
for (i in 1 : N) B <- rbind(B, apply(
summary(lm(Canopy_Height~Ground_Elevation, data=young400_1)) #use
data= instead of attach!
Or even
mylm <- lm(Canopy_Height~Ground_Elevation, data=young400_1)
mylm
summary(mylm)
coefficients(mylm)
Most intro to R guides cover the basics of modeling; you might benefit
from reading one of them.
S
Hello,
I am trying to run a fixed effects panel regression on data containing 5
columns and 1,494 rows.
I read the data in as follows:
>drugsXX<-read.csv(file="C:\\Folder\\vX.X\\Drugs\\drugsXX_panel.csv",
head=TRUE, sep=",")
Verified it read in correctly and had a good data.frame:
>dim(drugsXX
Thank you everyone.
Try executing this:
replicate(100, length(hist(rnorm(100), nclass = 10)$counts))
I'm still not sure why the number of bins (classes) is not consistent.
Thank in advance.
> Date: Wed, 4 Sep 2013 20:27:36 +0100
> From: ruipbarra...@sapo.pt
> To: pmassico...@hotmail.com
> CC: r-h
So look at the examples found in
?Control
and give it a try.
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 8/30/13 9:57 AM, "BJN1417" wrote:
>so I have to create a for loop of the geometric sequence
>h(x,n)=1+x+x^2+x^3^4
Good question. It turns out that the manual page does not tell
the whole story. Looking at the source code for hist.default,
the function starts with the number of breaks suggested by
nclass.Sturges(), but then this number (or any other number of
breaks that you specify) is passed to pretty() along
Hi,
Working with R, I often want to copy and paste some values somewhere else.
(Its not worth saving a CSV file for a dozen or so entries.) Or, I may want to
copy all the names of an object into some code.
R, rather nicely, wraps output with an index number on the left side.
For example:
[1
You should get no lines at all, as you have not specified that lines be
drawn. Use the "type" argument to do so.
xyplot(rnorm(5) ~1:5,pch=4) ## points only
xyplot(rnorm(5) ~1:5,pch=4,type="b",lwd=4) ## points with thick lines
read ?panel.xyplot carefully (the default panel function for xyplot) f
Hi,
You could use ?cat()
For ex:
vec1<-1:100
cat(vec1)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55
56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81
82 83 8
Hello,
can it be that xyplot does not support the lwd argument?
At least here, the following still shows thin lines, as opposed to the regular
plot command:
xyplot(Sepal.Length ~ Sepal.Width, data = iris, pch=4, lwd=4)
Cheers,
Daniel
--
Max-Planck-Institute for Dynamics and Self-Organization
> print(1:100) [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 16 17 18 19 20 21 22 23 24 25 26
[27] 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
43 44 45 46 47 48 49 50 51 52
[53] 53 54 55 56 57 58 59 60 61 62 63 64 65 66
Hello, All:
What do you recommend for reading a Google Spreadsheet into R? I didn't
find anything useful using library(sos); findFn('google spreadsheet').
I can solve the problem by downloading the file either as *.ods or
*.xlsx format, then opening it and saving it as *.xls, then using
re
On 04/09/2013 22:56, Noah Silverman wrote:
Hi,
Working with R, I often want to copy and paste some values somewhere else.
(Its not worth saving a CSV file for a dozen or so entries.) Or, I may want to
copy all the names of an object into some code.
R, rather nicely, wraps output with an ind
Depending on the OS you are working with awk or gawk are great utilities
for stripping columns from files. Also if you use a spreadsheet it is
quite easy to drop a column.
On Sep 4, 2013 5:59 PM, "Noah Silverman" wrote:
> Hi,
>
> Working with R, I often want to copy and paste some values somewhe
HI,
It is better to provide a reproducible example using ?dput().
you can also check in this link.
http://r.789695.n4.nabble.com/names-attribute-must-be-the-same-length-as-the-vector-td4503946.html
library(plm)
#Using the example from ?plm()
data("Produc", package = "plm")
zz <- plm(log(gsp) ~
Hi,
replicate(100,length(hist(10,0:10)$counts))
# [1] 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
10
#[26] 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
10
#[51] 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
10
#
On Wed, Sep 4, 2013 at 4:02 PM, philippe massicotte
wrote:
> Thank you everyone.
> Try executing this:
> replicate(100, length(hist(rnorm(100), nclass = 10)$counts))
> I'm still not sure why the number of bins (classes) is not consistent.
It depends on the range of x. If you look at the definitio
On 04/09/2013 4:02 PM, philippe massicotte wrote:
Thank you everyone.
Try executing this:
replicate(100, length(hist(rnorm(100), nclass = 10)$counts))
I'm still not sure why the number of bins (classes) is not consistent.
R is behaving as documented. You suggested 10 bins, but it finds for
so
First of I am new to using R.
I have a dataset that I plotted using R, I created a scatter plot and used
abline to create the line, what I need is to find the equation of the line.
Below is the script I have used up until this point.
>young400_1<-read.csv("Z:\\SOFTEL\\North Key Largo
project\\Can
Hello,
See the arguments 'right' and 'include.lowest' of ?hist.
To give what you want, try instead
h1 <- hist(1:10, 10) # counts are 2, 1, 1, ...
h2 <- hist(1:10, breaks = 0:10) # all counts are 1
and see the difference between h1 and h2, components 'breaks' and 'counts'.
Hope this helps,
We can just ask hist():
? hist
. . .
breaks
one of:
a vector giving the breakpoints between histogram cells,
a function to compute the vector of breakpoints,
a single number giving the number of cells for the
histogram,
=
Hi everyone.
I'm currently translating some Matlab code into R. However, I realized that the
hsit function produce different results in both languages.
in Matlab, hist(1:10, 10) will produce 10 bins with a count of 1 in each, but
in R it will produce 9 classes with count of 2,1,1,1,1,1,1,1,1.
I'm
On 04/09/2013 18:57, Hui Du wrote:
Hi R-community:
I heard 'snow' package is a good tool to parallelize processes and speed them
up. I tried to use it but was not successful. Could someboy point where I was
wrong? Thanks.
I want to read a HUGE file to R and hope 'snow' helps me to speed it up
Jean,
Thank you for the suggestion. Actually the dataset is quite large so that
method might be unmanageable.
Holly
hollym...@gmail.com
On Sep 4, 2013, at 10:14 AM, "Adams, Jean" wrote:
> Holly,
>
> I don't know of a clever way to do this, but I can think of a brute force
> way, which m
Hi R-community:
I heard 'snow' package is a good tool to parallelize processes and speed them
up. I tried to use it but was not successful. Could someboy point where I was
wrong? Thanks.
I want to read a HUGE file to R and hope 'snow' helps me to speed it up. Here
are codes:
library(snow)
iFi
Hi Simon,
What you need are regular expressions.
The help for gsub says this, but in such a way that if you didn't know
that's what you were looking for, you wouldn't learn it there:
See the help pages on regular expression for details of the
different types of regular expressions.
Th
>>> El día 01/09/2013 a las 15:13, Helen Sawaya
escribió:
> Dear R experts,
> I have data from a questionnaire that I would like to factor analyse.
It is
> in a likert scale form (0-3). I would like to check first for
univariate and
> multivariate outliers but the most common ways of doing so as
HI,
May be this helps.
input<- readLines(textConnection("
Private + Shared = RAM used Program
84.0 KiB + 14.5 KiB = 98.5 KiB sleep
108.0 KiB + 11.5 KiB = 119.5 KiB klogd
124.0 KiB + 15.0 KiB = 139.0 KiB hidd
128.0 KiB + 12.5 KiB = 140.5 KiB gpm
116.0 KiB + 28.5
Hi,
This is confusing because the error message suggests that you were not able to
read the csv file. Then how did you removed the Inf values?. I guess this
should be a different dataset.
>From your previous email:
"
prices=read.zoo("C:\\Users\\Desktopawc_au.csv",header=TRUE,sep=",",f
Rui's suggestion and your reproducible example work for me. Can you
give us more information about what you did, including the output of
sessionInfo() ?
library(dismo)
ex<-structure(list(env1 = c(182, 163.33, 443.02, 1240.16), env2 = c(1134,
550, 2111, 2523), env3 = c(24.53, 24.93, 24.71, 21.05),
Actually you have two loops, the for() loop you created and the
loop that is hidden inside apply(). You can hide the first loop
with lapply() or sapply():
B <- do.call(rbind, lapply(1:N, function(x)
colSums(A[sample.int(nrow(A), M, replace=TRUE),])))
Or
B <- t(sapply(1:N, function(x) colSums(A[s
Hello,
Where does the function bioclim come from? What package? If it's from
package dismo, then you should try predict() with the arguments reversed:
pred <- predict(mod, ex)
Hope this helps,
Rui Barradas
Em 04-09-2013 13:48, Kristi Glover escreveu:
Dear R User,
Would you give me some hi
Hi Rui,
Thanks for the message. i used dismo package.
I changed the argument as you suggested, still did not work.
Error in UseMethod("predict") :
no applicable method for 'predict' applied to an object of class "data.frame"
Thanks
KG
> Date: Wed, 4 Sep 2013 14:57:45 +0100
> From: ruipbarra.
Thanks Sarah and Rui, It worked now. I don't know why it was not working
before. I again restarted R and then it worked.
I really appreciated your help.
Cheers,
KG
> Date: Wed, 4 Sep 2013 10:20:39 -0400
> Subject: Re: [R] would you give me hints in r?
> From: sarah.gos...@gmail.com
> To: kristi
> -Original Message-
> When I try to apply mean to a list, I get the answer :
>
> argument is not numeric or logical: returning NA
>
Example:
l4 <- list(1:4)
class(l4) #not numeric or logical ...
mean(l4) #same error
#a list is not a number, a logical (TRUE/FALSE) or a vector or arra
Holly,
I don't know of a clever way to do this, but I can think of a brute force
way, which might only be feasible if you have a small data set (as in your
example). You could permute every possible set of connections, then choose
from that collection only the ones that meet your criteria.
Using
Hello,
I'm not sure I understand, but if you want a ?regexp to only match
numbers before a %, try the following.
gsub("[0-9]+%", "[percentagevalue]", text1)
gsub("[0-9]+%", "[percentagevalue]", text2)
[0-9] matches any character in the range from 0 to 9, and the + means to
repeat that char
> This is not efficient.
For whom?
> (a) L-BFGS-B does not identify itself even though it was called overriding
> the method
> parameter in optim.
Would you prefer that the warning
> o <- optim(par=c(1,2), fn=function(x)-sum(abs(sin(x))), method="CG",
lower=c(-1,-1), upper= c(2,3))
Warnin
Hi,
gsub("#%", "[percentagevalue]", text1)
#[1] "this is a number 23%"
gsub("\\d+%$", "[percentagevalue]", text1)
#[1] "this is a number [percentagevalue]"
gsub("bla", "23", text2)
#[1] "this is not a number 23%"
A.K.
- Original Message -
From: Simon Pickert
To: r-help@r-proje
Dear all,
I met some problems trying to install the TRR package.
I runed the command : install.packages("TRR")
I've received the following message :
In getDependencies(pkgs, dependencies, available, lib) :
package ‘TRR’ is not available (for R version 2.14.1)
I'm in Linux Mint and it seems
Greetings,
I have recently been exploring the 'glmnet' package and subsequently
cv.glmnet. The basic code as follows:
model <- cv.glmnet(variables, group, family="multinomial", alpha=.5,
standardize=F)
I understand that cv.glmnet does k-fold cross-validation to return a value
of lambda. Howeve
The latest release (2013-05-16, Good Sport) R-3.0.1 so perhaps you need to
upgrade to 3.0.1?
John Kane
Kingston ON Canada
> -Original Message-
> From: lucien.blanden...@unine.ch
> Sent: Wed, 4 Sep 2013 15:05:03 +
> To: r-help@r-project.org
> Subject: [R] Problem with installing the
Dear all,
I met some problems trying to install the TRR package.
I runed the command : install.packages("TRR")
I've received the following message :
In getDependencies(pkgs, dependencies, available, lib) :
package ‘TRR’ is not available (for R version 2.14.1)
I'm in Linux Mint and it seems
Hi all,
what are the placeholders for string operations/modifications? Is there a
placeholder for numbers, which would allow me to easily replace all numbers in
a string? Something like
text1 <- c("this is a number 23%")
text2 <- c("this is not a number bla%")
newtext1 <- gsub(#%, [percentagev
Hi there,
I am trying to apply multinomial Logit for a panel data set. I have 5016
observations for 22 countries (id). Each country has 228 observations over
time domain.
Following table shows the part of my dataset including 10 variables.
id t X1X2X3X4X5 X6
Hello everybody,
Without any loop and any package,
I would like to return N products of M rows in a matrix A :
Today, I managed to do it with a loop :
B <- matrix(NA, ncol = ncol(A), nrow = 0)
for (i in 1 : N) B <- rbind(B, apply(A[sample(1 : nrow(A), M, replace = T),
], 2, prod))
Do you have
thanks Pascal. actually, right after I sent the post I realized that I
didn't send my sessionInfo (at the bottom of this message now). Just
to make sure, everything in my script works except for the line with
stemDocument, which was commented in my gist. Below is the specific
warning I am getting:
Hello,
I noticed the commented line and I get no error when I run this specific
line.
Regards,
Pascal
2013/9/4 Ricardo Pietrobon
> thanks Pascal. actually, right after I sent the post I realized that I
> didn't send my sessionInfo (at the bottom of this message now). Just
> to make sure, eve
To give a specific example, the simple code for my test suite is given at the bottom of
this message. A simpler (simple-minded maybe) approach than creating a new packge for
testing. I now run this on the survival package every time that I submit a new version to
CRAN. It takes a while, since
Dear R User,
Would you give me some hints on why I could not predict using data. format'
data.
Here is the example:
ex<-structure(list(env1 = c(182, 163.33, 443.02, 1240.16), env2 = c(1134,
550, 2111, 2523), env3 = c(24.53, 24.93, 24.71, 21.05), env4 = c(0.05,
0, 0, 0)), .Names = c("env1", "en
On Tue, Sep 3, 2013 at 2:51 AM, Christoph Scherber
wrote:
> Dear all,
>
> I´ve played around with the "airquality" dataset, trying to solve the matrix
> equations of a simple
> multiple regression by hand; however, my matrix multiplications don´t lead to
> the estimates returned
> by coef().
Thank you for the reply, I wount post any other link anymore. I did not
know it's not allowed.
-
- László-András Zsurzsa,-
- Msc. Infromatics, Technical University Mu
I'll skip over the courtesy implications of double posting/pointing to
stackoverflow.
The stackoverflow thread makes it look as if you need to learn more
Excel. Do you really not know what an Excel template is?
It sounds as if you want what Excel calls "conditional formatting" which
you can
Hi Krzysztof,
Have a look at the packages sem, lavaan and psych.
Regards,
José
Prof. José Iparraguirre
Chief Economist
Age UK
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of fronc...@vizja.pl
Sent: 03 September 2013 13:33
To: r-h
Hi,
I have tried the ideas with an actual data set but couldn't
pass the parsing phase. The name of the 'Program' varies. MiB and KiB are
both included.
I should have shown the real-time data set.
Private + Shared = RAM used Program
84.0 KiB + 14.5 KiB = 98.5 Ki
Hi Christoph,
ginv() computes the Moore-Penrose generalized inverse by way of a
singular value decomposition. Part of the calculation involves taking
the reciprocal of the non zero values. In practice, non zero is
really "within some precision tolerance of zero". Numerical precision
can bite yo
Here is a solution I applied using qAIC and package bbmle so I share it
for next ones. It is not really automatized as I need to read every
results of the drop() test an enter manually the less significant
variable but I guess a function can be created in this goal.
nullQ <- update (null, fami
http://stackoverflow.com/questions/18511249/excel-cell-coloring-using-xlsx
This is the initial post on stackoverflow. Please look at this maybe I'm
clearer here.
Thank you in advance,
-
- László-András Zsurzsa,
Hi Michael,
You do not need to create a self-contained example from the mass of
code where it is embedded, but given that optim() works in many cases,
to file a bug report, you do need to give _an_ example where it is
failing.
Here is an example where it works great:
> optim(1, fn = function(x)
It would take some effort to extract selfcontained code from the mass of code
wherein this optimization is embedded. Moreover I would have to obtain
permission from my employer to do so.
This is not efficient.
However some things are evident from the trace log which I have submitted:
(a) L-BFGS-
Well, abusing the system does not appear to be working, does it? The correct
approach is to be more targeted in your audience. Is your problem with Ubuntu?
Perhaps you need to post in the r-sig-debian mailing list. Is your issue
specifically with the openNLS package? Try a simplified parallel pr
Hi
see
?contrasts
?model.matrix
go through archives
and through chapter
11.1.1 Contrasts
from R-Intro document.
Anyway, when you go through this chapter you will probably benefit from reading
previous chapters too.
Regards
Petr
> -Original Message-
> From: r-help-boun...@r-projec
I posted on both the forums because i needed urgent help on the issue. And
I am still waiting for a relevant reply on how to solve the issue. If you
have any suggestions please do share.
On Wed, Sep 4, 2013 at 6:08 AM, David Winsemius wrote:
> This was also posted on StackOverflow. Crossposting
Dear Arun,
Thanks a ton for your help.
Regards,
Anupam
On Tue, Sep 3, 2013 at 7:12 PM, arun wrote:
>
>
> HI,
>
> ?list.files()
> list.files() #created 4 files in my working directory
> #[1] "A_hubs.txt""A_nonhubs.txt" "B_hubs.txt""B_nonhubs.txt"
>
> #If you want to do wilcox.test in
77 matches
Mail list logo
