Dear Arun,
Thanks a ton for your help.
Regards,
Anupam
On Tue, Sep 3, 2013 at 7:12 PM, arun smartpink...@yahoo.com wrote:
HI,
?list.files()
list.files() #created 4 files in my working directory
#[1] A_hubs.txtA_nonhubs.txt B_hubs.txtB_nonhubs.txt
#If you want to do
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 dwinsem...@comcast.netwrote:
This was also posted on
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:
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
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)
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)
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 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 you
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
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:
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
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
On Tue, Sep 3, 2013 at 2:51 AM, Christoph Scherber
christoph.scher...@agr.uni-goettingen.de 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
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, env2,
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,
Hello,
I noticed the commented line and I get no error when I run this specific
line.
Regards,
Pascal
2013/9/4 Ricardo Pietrobon pietr...@gmail.com
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
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 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 a
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
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(#%, [percentagevalue],
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 TRR
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. However,
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,
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 simon.pick...@t-online.de
To:
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))
Warning
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 character
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.
-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 array of
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:
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:
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
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)
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),
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:
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
El día 01/09/2013 a las 15:13, Helen Sawaya
helensaw...@hotmail.com
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
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.
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)
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 jvad...@usgs.gov wrote:
Holly,
I don't know of a clever way to do this, but I can think of a brute force
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
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.
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,
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,
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
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
On Wed, Sep 4, 2013 at 4:02 PM, philippe massicotte
pmassico...@hotmail.com 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
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
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) ~
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 noahsilver...@ucla.edu wrote:
Hi,
Working with R, I often want to copy and paste
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
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
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,
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
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
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) for
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:
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()
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 bjn1...@uncw.edu wrote:
so I have to create a for loop of the geometric sequence
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)
[1]
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.
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,
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
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
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. It turns
Hi Spencer,
Why don't you want to install 64bit Java?
On Wed, Sep 4, 2013 at 6:12 PM, Spencer Graves
spencer.gra...@structuremonitoring.com 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
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
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
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
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
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
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)
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...
73 matches
Mail list logo
