Re: [R] Updating to nlme 3.1-62 failing from source (OS X)
This is an internal compiler/assembler error - please check that your tools are consistent and if so report this as a MacOS X error. (It is the sort of thing which happens if gcc targetted for one assembler is used with another, for example on Solaris with GNU vs Sun assemblers.) Given that the C code being compiled is unchanged since 3.1-60 which was successfully compiled by the maintainers, this does point to a problem local to you. On Sun, 31 Jul 2005, Michael Kubovy wrote: R Version 2.1.1 (2005-06-20); R Cocoa GUI 1.12 (1622); Mac OS X 10.4.2 (8C46) For nlme, the R package installer for CRAN (binaries) gives 3.1-60 as Repository Version as well as Installed Version. The CRAN repository for MacOS X has not be updated for a couple of weeks. For CRAN (sources) it gives 3.1-62 as Repository Version and 3.1-60 as Installed Version. When I try to update nlme, I get * Installing *source* package 'nlme' ... ** libs gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c corStruct.c -o corStruct.o /var/tmp//ccOTBSDa.s:6088:** Removing '/Library/Frameworks/ R.framework/Versions/2.1.1/Resources/library/nlme' ** Restoring previous '/Library/Frameworks/R.framework/Versions/2.1.1/ Resources/library/nlme' The downloaded packages are in /private/tmp/RtmpL1j8Sa/downloaded_packages Unknown pseudo-op: .p2align /var/tmp//ccOTBSDa.s:6088:Rest of line ignored. 1st junk character valued 50 (2). /var/tmp//ccOTBSDa.s:6414:Unknown pseudo-op: .p2align /var/tmp//ccOTBSDa.s:6414:Rest of line ignored. 1st junk character valued 50 (2). /var/tmp//ccOTBSDa.s:6608:Unknown pseudo-op: .p2align /var/tmp//ccOTBSDa.s:6608:Rest of line ignored. 1st junk character valued 50 (2). /var/tmp//ccOTBSDa.s:7111:Unknown pseudo-op: .subsections_via_symbols make: *** [corStruct.o] Error 1 ERROR: compilation failed for package 'nlme' -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ 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] Conditional piece-wise dependent regression
Hi, after reading some R docs, I couldnt figure out how can I find the solution for the following problem, therefore I would ask this friendly list for an advice. Were making a least square approximation for an experiment described by the following model: T is the time, Y is some measured value. From time=0 till time=U: Y = b + p*T From time=U and on (some effect added): Y = b + p*T + q*(T-U) From time=V and on (some additional effect added): Y = b + p*T + q*(T-U) + r*(T-V) Measured: Yi, Ti pairs. Wanted: b, p, q, r. b and p are the same for all time ranges; q is the same for time=U and on. Thanks, Arie. __ 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] Conditional piece-wise dependent regression
Since you want least squares, I think you could use lm() here, i.e., U - 50 V - 100 Time - 1:150 dat - data.frame(y = rnorm(150), Time, f1 = as.numeric(Time U), f2 = as.numeric(Time V)) ### m - lm(y ~ Time + I(Time - U):f1 + I(Time - V):f2, data = dat) # check also the design matrix model.matrix(m) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.be/biostat/ http://www.student.kuleuven.be/~m0390867/dimitris.htm - Original Message - From: Arie [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Monday, August 01, 2005 10:17 AM Subject: [R] Conditional piece-wise dependent regression Hi, after reading some R docs, I couldn't figure out how can I find the solution for the following problem, therefore I would ask this friendly list for an advice. We're making a least square approximation for an experiment described by the following model: T is the time, Y is some measured value. From time=0 till time=U: Y = b + p*T From time=U and on (some effect added): Y = b + p*T + q*(T-U) From time=V and on (some additional effect added): Y = b + p*T + q*(T-U) + r*(T-V) Measured: Yi, Ti pairs. Wanted: b, p, q, r. b and p are the same for all time ranges; q is the same for time=U and on. Thanks, Arie. __ 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-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] fortran
Dear Friends I am trying to call a fortran subroutine in R, but i can not. How can i use the subroutines in R? Regrads Javad __ [[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
Re: [R] fortran
See ?dyn.load ?.Fortran and manual Writing R Extensions! Then ask again if you have additional questions! Hope this helps, Ales Ziberna - Original Message - From: javad Ashjari [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Sent: Monday, August 01, 2005 12:21 PM Subject: [R] fortran Dear Friends I am trying to call a fortran subroutine in R, but i can not. How can i use the subroutines in R? Regrads Javad __ [[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 __ 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] How to hiding code for a package
I always thought that ability to see and study the code of every package, was a great thing about R and other open source environments. So I hope there are no good ways of hiding the code of packages. Jarek \=== Jarek Tuszynski, PhD. o / \ Science Applications International Corporation \__,| (703) 676-4192 \ [EMAIL PROTECTED] `\ -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gary Wong Sent: Saturday, July 30, 2005 4:21 AM To: r-help@stat.math.ethz.ch Subject: [R] How to hiding code for a package Hey everyone, I have made a package and wish to release it but before then I have a problem. I have a few functions in this package written in R that I wish to hide such that after installation, someone can use say the function foo(parameters = ) but cannot do foo. Typing foo should not show the source code or at least not all of it. Is there a way to do this ? I have searched the mailing list and used google, and have found something like [R] Hiding internal package functions for the doc. pkg-internal.Rd but this seems different since it seems that the keyword internal just hides the function from showing in the index and hides documentation, not the function itself. Can someone help? Thanks Gary __ 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-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] converting stata's by syntax to R
I am struggling with migrating some stata code to R. I have a data
frame containing, sometimes, repeat observations (rows) of the same
family. I want to keep only one observation per family, selecting
that observation according to some other variable. An example data
frame is:
# construct example data
fam - c(1,2,3,3,4,4,4)
wt - c(1,1,0.6,0.4,0.4,0.4,0.2)
keep - c(1,1,1,0,1,0,0)
dat - as.data.frame(cbind(fam,wt,keep))
dat
I want to keep the observation for which wt is a maximum, and where
this doesn't identify a unique observation, to keep just one anyway,
not caring which. Those observations are indicated above by keep==1.
(Note, keep - c(1,1,1,0,0,1,0) would be fine too, but not
c(1,1,1,0,0,0,1)).
The stata code I would use is
bys fam (wt): keep if _n==_N
This is my (long-winded) attempt in R:
# first keep those rows where wt=max_fam(wt)
maxwt - by(dat,dat$fam,function(x) max(x[,2]))
maxwt - sapply(maxwt,[[,1)
maxwt.dat - data.frame(maxwt=maxwt,fam=as.integer(names(maxwt)))
dat - merge(dat,maxwt.dat)
dat - dat[dat$wt==dat$maxwt,]
dat
Now I am stuck - I want to keep either row with fam==4, and have tried
playing around with combinations of sample and apply or by, but with
no success. I can only find an inefficient for-loop solution:
# identify those rows with 1 observation
more - by(dat,dat$fam,function(x) dim(x)[1])
more - sapply(more,[[,1)
more.dat - data.frame(more=more,fam=as.integer(names(more)))
dat - merge(dat,more.dat)
# sample from those for whom more1
result-dat[dat$more==1,]
for(f in unique(dat$fam[dat$more1])) {
rows - rownames(dat[dat$fam==f,])
result - rbind(result,dat[sample(rows,1),])
}
result
I am sure that for something so simple in stata to be so complicated
in R must indicate ignorance of R on my part, but searches of help
files and RSiteSearch hasn't led to any better solution.
Any suggestions would be most helpful! Thanks, C.
__
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] reading from posgresql
Dear all I don't know how to read data from posgresql to R. Does anybody can help me? Thanks a lot Adrián __ 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] reading from posgresql
I don't know how to read data from posgresql to R. The RdbiPgSQL package will read from Postgres. Its current incarnation is available via Bioconductor (www.bioconductor.org). __ 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] converting stata's by syntax to R
try
attach(dat)
dat-dat[order(fam,wt),]
#sort the data ,as the stata's byable command does
lis-by(dat,fam,function(x) x[length(x$fam),])
#equall your stata command ,but return a list.
do.call(rbind,lis)
#to make the list to be a matrix-like result.
fam wt keep
1 1 1.01
2 2 1.01
3 3 0.40
4 4 0.40
=== 2005-08-01 22:24:27 您在来信中写道：===
I am struggling with migrating some stata code to R. I have a data
frame containing, sometimes, repeat observations (rows) of the same
family. I want to keep only one observation per family, selecting
that observation according to some other variable. An example data
frame is:
# construct example data
fam - c(1,2,3,3,4,4,4)
wt - c(1,1,0.6,0.4,0.4,0.4,0.2)
keep - c(1,1,1,0,1,0,0)
dat - as.data.frame(cbind(fam,wt,keep))
dat
I want to keep the observation for which wt is a maximum, and where
this doesn't identify a unique observation, to keep just one anyway,
not caring which. Those observations are indicated above by keep==1.
(Note, keep - c(1,1,1,0,0,1,0) would be fine too, but not
c(1,1,1,0,0,0,1)).
The stata code I would use is
bys fam (wt): keep if _n==_N
This is my (long-winded) attempt in R:
# first keep those rows where wt=max_fam(wt)
maxwt - by(dat,dat$fam,function(x) max(x[,2]))
maxwt - sapply(maxwt,[[,1)
maxwt.dat - data.frame(maxwt=maxwt,fam=as.integer(names(maxwt)))
dat - merge(dat,maxwt.dat)
dat - dat[dat$wt==dat$maxwt,]
dat
Now I am stuck - I want to keep either row with fam==4, and have tried
playing around with combinations of sample and apply or by, but with
no success. I can only find an inefficient for-loop solution:
# identify those rows with 1 observation
more - by(dat,dat$fam,function(x) dim(x)[1])
more - sapply(more,[[,1)
more.dat - data.frame(more=more,fam=as.integer(names(more)))
dat - merge(dat,more.dat)
# sample from those for whom more1
result-dat[dat$more==1,]
for(f in unique(dat$fam[dat$more1])) {
rows - rownames(dat[dat$fam==f,])
result - rbind(result,dat[sample(rows,1),])
}
result
I am sure that for something so simple in stata to be so complicated
in R must indicate ignorance of R on my part, but searches of help
files and RSiteSearch hasn't led to any better solution.
Any suggestions would be most helpful! Thanks, C.
__
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
= = = = = = = = = = = = = = = = = = = =
2005-08-01
--
Deparment of Sociology
Fudan University
Blog:http://sociology.yculblog.com
__
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] converting stata's by syntax to R
Chris Wallace [EMAIL PROTECTED] writes:
I am struggling with migrating some stata code to R. I have a data
frame containing, sometimes, repeat observations (rows) of the same
family. I want to keep only one observation per family, selecting
that observation according to some other variable. An example data
frame is:
# construct example data
fam - c(1,2,3,3,4,4,4)
wt - c(1,1,0.6,0.4,0.4,0.4,0.2)
keep - c(1,1,1,0,1,0,0)
dat - as.data.frame(cbind(fam,wt,keep))
dat
I want to keep the observation for which wt is a maximum, and where
this doesn't identify a unique observation, to keep just one anyway,
not caring which. Those observations are indicated above by keep==1.
(Note, keep - c(1,1,1,0,0,1,0) would be fine too, but not
c(1,1,1,0,0,0,1)).
The stata code I would use is
bys fam (wt): keep if _n==_N
This is my (long-winded) attempt in R:
# first keep those rows where wt=max_fam(wt)
maxwt - by(dat,dat$fam,function(x) max(x[,2]))
maxwt - sapply(maxwt,[[,1)
maxwt.dat - data.frame(maxwt=maxwt,fam=as.integer(names(maxwt)))
dat - merge(dat,maxwt.dat)
dat - dat[dat$wt==dat$maxwt,]
dat
Now I am stuck - I want to keep either row with fam==4, and have tried
playing around with combinations of sample and apply or by, but with
no success. I can only find an inefficient for-loop solution:
# identify those rows with 1 observation
more - by(dat,dat$fam,function(x) dim(x)[1])
more - sapply(more,[[,1)
more.dat - data.frame(more=more,fam=as.integer(names(more)))
dat - merge(dat,more.dat)
# sample from those for whom more1
result-dat[dat$more==1,]
for(f in unique(dat$fam[dat$more1])) {
rows - rownames(dat[dat$fam==f,])
result - rbind(result,dat[sample(rows,1),])
}
result
I am sure that for something so simple in stata to be so complicated
in R must indicate ignorance of R on my part, but searches of help
files and RSiteSearch hasn't led to any better solution.
Any suggestions would be most helpful! Thanks, C.
How about
unsplit(lapply(split(dat,dat$fam),
function(x) seq(length=nrow(x)) == which.max(x$wt)),
dat$fam)
or
do.call(rbind, lapply(split(dat,dat$fam),
function(x) x[which.max(x$wt),]))
or (same thing, basically)
do.call(rbind, by(dat,dat$fam,function(x) x[which.max(x$wt),]))
--
O__ Peter Dalgaard Øster Farimagsgade 5, Entr.B
c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907
__
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] converting stata's by syntax to R
Here is one way this can be done
do.call(rbind, by(dat, list(dat$fam) ,function(x) {
+ if(NROW(x)1) return(x[which.max(x$wt),])
+ else return(x)}
+ ))
and it returns
fam wt keep
1 1 1.01
2 2 1.01
3 3 0.61
4 4 0.41
hth,
On Mon, 1 Aug 2005, Chris Wallace wrote:
I am struggling with migrating some stata code to R. I have a data
frame containing, sometimes, repeat observations (rows) of the same
family. I want to keep only one observation per family, selecting
that observation according to some other variable. An example data
frame is:
# construct example data
fam - c(1,2,3,3,4,4,4)
wt - c(1,1,0.6,0.4,0.4,0.4,0.2)
keep - c(1,1,1,0,1,0,0)
dat - as.data.frame(cbind(fam,wt,keep))
dat
I want to keep the observation for which wt is a maximum, and where
this doesn't identify a unique observation, to keep just one anyway,
not caring which. Those observations are indicated above by keep==1.
(Note, keep - c(1,1,1,0,0,1,0) would be fine too, but not
c(1,1,1,0,0,0,1)).
The stata code I would use is
bys fam (wt): keep if _n==_N
This is my (long-winded) attempt in R:
# first keep those rows where wt=max_fam(wt)
maxwt - by(dat,dat$fam,function(x) max(x[,2]))
maxwt - sapply(maxwt,[[,1)
maxwt.dat - data.frame(maxwt=maxwt,fam=as.integer(names(maxwt)))
dat - merge(dat,maxwt.dat)
dat - dat[dat$wt==dat$maxwt,]
dat
Now I am stuck - I want to keep either row with fam==4, and have tried
playing around with combinations of sample and apply or by, but with
no success. I can only find an inefficient for-loop solution:
# identify those rows with 1 observation
more - by(dat,dat$fam,function(x) dim(x)[1])
more - sapply(more,[[,1)
more.dat - data.frame(more=more,fam=as.integer(names(more)))
dat - merge(dat,more.dat)
# sample from those for whom more1
result-dat[dat$more==1,]
for(f in unique(dat$fam[dat$more1])) {
rows - rownames(dat[dat$fam==f,])
result - rbind(result,dat[sample(rows,1),])
}
result
I am sure that for something so simple in stata to be so complicated
in R must indicate ignorance of R on my part, but searches of help
files and RSiteSearch hasn't led to any better solution.
Any suggestions would be most helpful! Thanks, C.
__
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-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] Clayton-Oakes failure time model
Dear R-list, does an R-package exist in which the Clayton-Oakes failure time model is implemented? Thanks for help, Stefan. [[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
Re: [R] converting stata's by syntax to R
if you also need to create the `keep' vector, then you could try this
approach:
fam - c(1,2,3,3,4,4,4)
wt - c(1,1,0.6,0.4,0.4,0.4,0.2)
dat - data.frame(fam, wt)
###
keep - unlist( lapply(split(wt, fam), function(x){
ind - rep(FALSE, length(x))
ind[which.max(x)] - TRUE
ind
}) )
as.numeric(keep)
dat[keep, ]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
Web: http://www.med.kuleuven.be/biostat/
http://www.student.kuleuven.be/~m0390867/dimitris.htm
- Original Message -
From: Chris Wallace [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Sent: Monday, August 01, 2005 4:24 PM
Subject: [R] converting stata's by syntax to R
I am struggling with migrating some stata code to R. I have a data
frame containing, sometimes, repeat observations (rows) of the same
family. I want to keep only one observation per family, selecting
that observation according to some other variable. An example data
frame is:
# construct example data
fam - c(1,2,3,3,4,4,4)
wt - c(1,1,0.6,0.4,0.4,0.4,0.2)
keep - c(1,1,1,0,1,0,0)
dat - as.data.frame(cbind(fam,wt,keep))
dat
I want to keep the observation for which wt is a maximum, and where
this doesn't identify a unique observation, to keep just one anyway,
not caring which. Those observations are indicated above by
keep==1.
(Note, keep - c(1,1,1,0,0,1,0) would be fine too, but not
c(1,1,1,0,0,0,1)).
The stata code I would use is
bys fam (wt): keep if _n==_N
This is my (long-winded) attempt in R:
# first keep those rows where wt=max_fam(wt)
maxwt - by(dat,dat$fam,function(x) max(x[,2]))
maxwt - sapply(maxwt,[[,1)
maxwt.dat -
data.frame(maxwt=maxwt,fam=as.integer(names(maxwt)))
dat - merge(dat,maxwt.dat)
dat - dat[dat$wt==dat$maxwt,]
dat
Now I am stuck - I want to keep either row with fam==4, and have
tried
playing around with combinations of sample and apply or by, but with
no success. I can only find an inefficient for-loop solution:
# identify those rows with 1 observation
more - by(dat,dat$fam,function(x) dim(x)[1])
more - sapply(more,[[,1)
more.dat - data.frame(more=more,fam=as.integer(names(more)))
dat - merge(dat,more.dat)
# sample from those for whom more1
result-dat[dat$more==1,]
for(f in unique(dat$fam[dat$more1])) {
rows - rownames(dat[dat$fam==f,])
result - rbind(result,dat[sample(rows,1),])
}
result
I am sure that for something so simple in stata to be so complicated
in R must indicate ignorance of R on my part, but searches of help
files and RSiteSearch hasn't led to any better solution.
Any suggestions would be most helpful! Thanks, C.
__
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-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] Conditional piece-wise dependent regression
In case you want to estimate U and V as well, there's a segmented regression package (called segmented) on http://cran.r-project.org for this. Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Dimitris Rizopoulos Sent: Monday, August 01, 2005 4:40 AM To: Arie Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Conditional piece-wise dependent regression Since you want least squares, I think you could use lm() here, i.e., U - 50 V - 100 Time - 1:150 dat - data.frame(y = rnorm(150), Time, f1 = as.numeric(Time U), f2 = as.numeric(Time V)) ### m - lm(y ~ Time + I(Time - U):f1 + I(Time - V):f2, data = dat) # check also the design matrix model.matrix(m) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: http://www.med.kuleuven.be/biostat/ http://www.student.kuleuven.be/~m0390867/dimitris.htm - Original Message - From: Arie [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Monday, August 01, 2005 10:17 AM Subject: [R] Conditional piece-wise dependent regression Hi, after reading some R docs, I couldn't figure out how can I find the solution for the following problem, therefore I would ask this friendly list for an advice. We're making a least square approximation for an experiment described by the following model: T is the time, Y is some measured value. From time=0 till time=U: Y = b + p*T From time=U and on (some effect added): Y = b + p*T + q*(T-U) From time=V and on (some additional effect added): Y = b + p*T + q*(T-U) + r*(T-V) Measured: Yi, Ti pairs. Wanted: b, p, q, r. b and p are the same for all time ranges; q is the same for time=U and on. Thanks, Arie. __ 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-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-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] How to hiding code for a package
There's something else you could try - since you can't hide the code, obfuscate it. Hide the real thing in a large pile of useless, complicated, awfully formatted code that would stop anyone except the most desperate (including yourself, after a couple of weeks/months) from trying to understand it. The best solution would be to compile the code, but R is not there yet. -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Saturday, July 30, 2005 5:35 AM To: Gary Wong Cc: r-help@stat.math.ethz.ch Subject: Re: [R] How to hiding code for a package What you ask is impossible. For a function to be callable it has to be locatable and hence can be printed. One possibility is to have a namespace, and something like foo - function(...) foo_internal(...) where foo is exported but foo_internal is not. Then foo_internal is hidden from casual inspection, but it can be listed by cognescenti. Why do you want to do this? Anyhone can read the source code of your package, and any function which can be called can be deparsed, possibly after jumping through a few hoops. On Sat, 30 Jul 2005, Gary Wong wrote: Hey everyone, I have made a package and wish to release it but before then I have a problem. I have a few functions in this package written in R that I wish to hide such that after installation, someone can use say the function foo(parameters = ) but cannot do foo. Typing foo should not show the source code or at least not all of it. Is there a way to do this ? I have searched the mailing list and used google, and have found something like [R] Hiding internal package functions for the doc. pkg-internal.Rd but this seems different since it seems that the keyword internal just hides the function from showing in the index and hides documentation, not the function itself. Can someone help? Thanks -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ 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-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] How to hiding code for a package
Dear Gary Wong: Are you trying to defeat the GNU License? If so, why not make it available as an S-Plus library? There, you don't have to provide the source; binary files will suffice. spencer graves bogdan romocea wrote: There's something else you could try - since you can't hide the code, obfuscate it. Hide the real thing in a large pile of useless, complicated, awfully formatted code that would stop anyone except the most desperate (including yourself, after a couple of weeks/months) from trying to understand it. The best solution would be to compile the code, but R is not there yet. -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Saturday, July 30, 2005 5:35 AM To: Gary Wong Cc: r-help@stat.math.ethz.ch Subject: Re: [R] How to hiding code for a package What you ask is impossible. For a function to be callable it has to be locatable and hence can be printed. One possibility is to have a namespace, and something like foo - function(...) foo_internal(...) where foo is exported but foo_internal is not. Then foo_internal is hidden from casual inspection, but it can be listed by cognescenti. Why do you want to do this? Anyhone can read the source code of your package, and any function which can be called can be deparsed, possibly after jumping through a few hoops. On Sat, 30 Jul 2005, Gary Wong wrote: Hey everyone, I have made a package and wish to release it but before then I have a problem. I have a few functions in this package written in R that I wish to hide such that after installation, someone can use say the function foo(parameters = ) but cannot do foo. Typing foo should not show the source code or at least not all of it. Is there a way to do this ? I have searched the mailing list and used google, and have found something like [R] Hiding internal package functions for the doc. pkg-internal.Rd but this seems different since it seems that the keyword internal just hides the function from showing in the index and hides documentation, not the function itself. Can someone help? Thanks -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ 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-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 -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA [EMAIL PROTECTED] www.pdf.com http://www.pdf.com Tel: 408-938-4420 Fax: 408-280-7915 __ 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] converting stata's by syntax to R
On Mon, 1 Aug 2005, Chris Wallace wrote: I am struggling with migrating some stata code to R. I have a data frame containing, sometimes, repeat observations (rows) of the same family. I want to keep only one observation per family, selecting that observation according to some other variable. An example data frame is: # construct example data fam - c(1,2,3,3,4,4,4) wt - c(1,1,0.6,0.4,0.4,0.4,0.2) keep - c(1,1,1,0,1,0,0) dat - as.data.frame(cbind(fam,wt,keep)) dat I want to keep the observation for which wt is a maximum, and where this doesn't identify a unique observation, to keep just one anyway, not caring which. Those observations are indicated above by keep==1. (Note, keep - c(1,1,1,0,0,1,0) would be fine too, but not c(1,1,1,0,0,0,1)). The stata code I would use is bys fam (wt): keep if _n==_N A reasonably direct translation of the Stata code is index - order(fam, -wt) keep - !duplicated(fam[index]) dat - data.frame(fam=fam[index], wt=wt[index], keep=keep) which sorts wt into decreasing order within family, then keeps the first observation in each family. This is less general than solutions other people have given, but I'd expect it to be faster for large data sets. 'keep' ends up TRUE/FALSE rather than 1/0; if this is a problem use as.numeric() on it. -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
Re: [R] Compiling nlme 3.1-62 succeeds if OS X Developer Tools updated; however, compiling mgcv_1.3-4 fails
Prof. Ripley was right. The problem went away when I updated the Developer Tools. On Aug 1, 2005, at 2:48 AM, Prof Brian Ripley wrote: This is an internal compiler/assembler error - please check that your tools are consistent and if so report this as a MacOS X error. On Sun, 31 Jul 2005, Michael Kubovy wrote: When I try to update nlme, I get * Installing *source* package 'nlme' ... ** libs gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c corStruct.c -o corStruct.o /var/tmp//ccOTBSDa.s:6088:** Removing '/Library/Frameworks/ R.framework/Versions/2.1.1/Resources/library/nlme' ** Restoring previous '/Library/Frameworks/R.framework/Versions/ 2.1.1/ Resources/library/nlme' The downloaded packages are in /private/tmp/RtmpL1j8Sa/downloaded_packages Unknown pseudo-op: .p2align /var/tmp//ccOTBSDa.s:6088:Rest of line ignored. 1st junk character valued 50 (2). /var/tmp//ccOTBSDa.s:6414:Unknown pseudo-op: .p2align /var/tmp//ccOTBSDa.s:6414:Rest of line ignored. 1st junk character valued 50 (2). /var/tmp//ccOTBSDa.s:6608:Unknown pseudo-op: .p2align /var/tmp//ccOTBSDa.s:6608:Rest of line ignored. 1st junk character valued 50 (2). /var/tmp//ccOTBSDa.s:7111:Unknown pseudo-op: .subsections_via_symbols make: *** [corStruct.o] Error 1 ERROR: compilation failed for package 'nlme' However, the following problem occurred when I tried to compile mgcv: * Installing *source* package 'mgcv' ... ** libs gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c gcv.c -o gcv.o gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c magic.c -o magic.o gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c mat.c -o mat.o gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c matrix.c -o matrix.o gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c mgcv.c -o mgcv.o gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c qp.c -o qp.o gcc-3.3 -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/ include -I/usr/local/include -fno-common -g -O2 -c tprs.c -o tprs.o gcc-3.3 -bundle -flat_namespace -undefined suppress -L/usr/local/lib - o mgcv.so gcv.o magic.o mat.o matrix.o mgcv.o qp.o tprs.o -framework vecLib -L/usr/local/lib/gcc/powerpc-apple-darwin6.8/3.4.2 -lg2c - lSystem -framework R ld: warning ** Removing '/Library/Frameworks/R.framework/Versions/ 2.1.1/Resources/library/mgcv' ** Restoring previous '/Library/Frameworks/R.framework/Versions/2.1.1/ Resources/library/mgcv' multiple definitions of symbol _xerbla_ /Library/Frameworks/R.framework/R(print.lo) definition of _xerbla_ /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/ vecLib.framework/Versions/A/libBLAS.dylib(single module) definition of _xerbla_ ld: warning multiple definitions of symbol _signgam /Library/Frameworks/R.framework/R(lgamma.lo) definition of _signgam /usr/lib/libSystem.dylib(gamma9.o) definition of _signgam ld: gcv.o has external relocation entries in non-writable section (__TEXT,__text) for symbols: restFP saveFP make: *** [mgcv.so] Error 1 ERROR: compilation failed for package 'mgcv' _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751 Fax:+1-434-982-4766 WWW:http://www.people.virginia.edu/~mk9y/ __ 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] How to hiding code for a package
--- Spencer Graves [EMAIL PROTECTED] wrote: Dear Gary Wong: Are you trying to defeat the GNU License? If so, why not make it available as an S-Plus library? There, you don't have to provide the source; binary files will suffice. spencer graves bogdan romocea wrote: There's something else you could try - since you can't hide the code, obfuscate it. Hide the real thing in a large pile of useless, complicated, awfully formatted code that would stop anyone except the most desperate (including yourself, after a couple of weeks/months) from trying to understand it. The best solution would be to compile the code, but R is not there yet. -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Saturday, July 30, 2005 5:35 AM To: Gary Wong Cc: r-help@stat.math.ethz.ch Subject: Re: [R] How to hiding code for a package What you ask is impossible. For a function to be callable it has to be locatable and hence can be printed. One possibility is to have a namespace, and something like foo - function(...) foo_internal(...) where foo is exported but foo_internal is not. Then foo_internal is hidden from casual inspection, but it can be listed by cognescenti. Why do you want to do this? Anyhone can read the source code of your package, and any function which can be called can be deparsed, possibly after jumping through a few hoops. On Sat, 30 Jul 2005, Gary Wong wrote: Hey everyone, I have made a package and wish to release it but before then I have a problem. I have a few functions in this package written in R that I wish to hide such that after installation, someone can use say the function foo(parameters = ) but cannot do foo. Typing foo should not show the source code or at least not all of it. Is there a way to do this ? I have searched the mailing list and used google, and have found something like [R] Hiding internal package functions for the doc. pkg-internal.Rd but this seems different since it seems that the keyword internal just hides the function from showing in the index and hides documentation, not the function itself. Can someone help? Thanks -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ 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-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 -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA [EMAIL PROTECTED] www.pdf.com http://www.pdf.com Tel: 408-938-4420 Fax: 408-280-7915 __ 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-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] New functions supporting GIF file format in R
Sorry about late response but I was traveling and did not follow R-Help for 2 weeks. I think, I prefer to leave read/write.GIF functions in my package where I have a full control of IO formats. I like to store images as integer or real matrices, and movies as 3D arrays. I believe that 'pixmap' package (as well as 'rimage' package and possibly others) use specialized classes for images , and do not handle movies at all (correct me if I am wrong). If you would like to add capability reading and writing GIF files (or ENVI files - another multi-frame file format - see read/write.ENVI) to the 'pixmap' class, feel free to write a shell around my function. With input in correct range (0:255) and option scale='never', no processing will be done to your data on the way in or out. Also, I am not planning on making any changes to that function, unless bugs are found. Regards, Jarek \=== Jarek Tuszynski, PhD. o / \ Science Applications International Corporation \__,| (703) 676-4192 \ [EMAIL PROTECTED] `\ -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Martin Maechler Sent: Tuesday, July 19, 2005 2:53 AM To: Tuszynski, Jaroslaw W. Cc: ([EMAIL PROTECTED]) Subject: Re: [R] New functions supporting GIF file format in R JarekT == Tuszynski, Jaroslaw W [EMAIL PROTECTED] on Mon, 18 Jul 2005 16:00:43 -0400 writes: JarekT Hi, A minor announcement. I just added two functions JarekT for reading and writing GIF files to my caTools JarekT package. Input and output is in the form of standard JarekT R matrices or arrays, and standard R color-maps JarekT (palettes). The functions can read and write both JarekT regular GIF images, as well as, multi-frame animated JarekT GIFs. Most of the work is done in C level code JarekT (included), so functions do not use any external JarekT libraries. JarekT For more info and examples go to JarekT http://cran.r-project.org/doc/packages/caTools.pdf JarekT http://cran.r-project.org/doc/packages/caTools.pdf JarekT and click GIF. Wouldn't it make sense to donate these to the 'pixmap' package which is dedicated to such objects and has been in place for a very long time? Regards, Martin __ 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-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] NA when using read.csv
Hi, when I used: Dist=read.csv(test.csv,header=TRUE) to read data from CSV file. For some cells, R mistakenly put in as NA, while most of the cells still appears to be right, and there is no error message in R. I am pretty sure the csv file is all right, but just can't figure out what went wrong. Can someone share your thoughts with me? Thanks! Ed __ 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] NA when using read.csv
On 8/1/05, Haibo Huang [EMAIL PROTECTED] wrote: Hi, when I used: Dist=read.csv(test.csv,header=TRUE) to read data from CSV file. For some cells, R mistakenly put in as NA, while most of the cells still Its not likely that there are errors in this software so its probably not a mistake but some other cause. appears to be right, and there is no error message in R. I am pretty sure the csv file is all right, but just can't figure out what went wrong. Can someone share your thoughts with me? Thanks! Read the posting guide at the bottom of each post. In particular, cut down the input file to a few lines that still exhibit the problem in order to provide a reproducible example. As an aside note that 'header = TRUE' is the default in read.csv (but not on read.table) so you don't need to specify it. __ 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] canonical correlations on multiple sets
Hi! I was wondering whether anybody could help me with canonical correlation on multiple sets (3 or 4) of variates, e.g. using one of the methods suggested by Kettenring (Biometrika 58(3):433, 1971), starting from covariance or correlation matrices. If someone could point me to an efficient way to achieve this or maybe share a method, I would be forever grateful. Some additional information on the data: Submatrices of covariance matrix S: S11, S12, S13, (S14, ) S22, S23, (S24, ) S33 (, S34, S44), with dim(S) = p + q + q or p + p + q + q I would like to retrieve the canonical correlates and the canonical variates for all sets. Thanks in advance, Tsjerk __ 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] Rgdal windows binary warning message
Hi all, I just downloaded windows binaries of RGDAL (from http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.1/) and installed. I ran the example R_HOME\library\rgdal\R-ex\getPixmapGDAL.R and everything seemed to work as expected. However I got the following warning message: DLL attempted to change FPU control word from 8001f to 9001f The R-help (?dyn.load) describes the problem as: External code must not change the floating point control word, but many DLLs do so. Common changes are to set it to use 53 bit precision instead of R's default 64 bit precision, or to unmask some exceptions. 'dyn.load' detects such changes, and restores R's control word to its default value of hex 8001F. This may cause the DLL to malfunction; if so, it should be rewritten to save and restore the control word itself. If 'warn.FPU' is set to 'TRUE' using the 'options' function, a warning will be printed. (If the warning says that the control word was changed from some other value than 8001F, please report the circumstances to the Windows maintainers: that probably indicates an internal bug.) Does anyone know if this is a major problem for reading and manipulating image files using RDAL? And if so are there windows binaries which have patched it up? Thanks in advance Tony ** Tony Gill - PhD Candidate Centre for Remote Sensing Spatial Information Science School of Geography, Planning Architecture University of Queensland Brisbane, Queensland, AUSTRALIA, 4072 Ph: 61-7-3365-7027 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
[R] can we manage memory usage to increase speed?
Hi, Thanks for reading. I am running a process in R for microarray data analysis. RedHat Enterprise Linux 4, dual AMD CPU, 6G memory. However, the R process use only a total of 200M memory. And the CPU usage is total to ~110% for two. The program takes at least 2 weeks to run at the current speed. Is there some way we can increase the usage of CPUs and memories and speed up? Any suggestion is appreciated. Thanks again. Zhilin [[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
Re: [R] nnet package for OS X
Hello every one-- My sincere apologies if I have missed something obvious, but I need to construct some multinomial logit models. I remember using nnet for this in the past, but I am having trouble locating a OS X version of this (it's been a year or so since I had to use it). I found a Windows version, but when I try to install across platforms, it crashes R immediately--no error messages or anything helpful. Can any one direct me to a repository, or alternatively, to another package readily available Mac side that can do simple multinomial logit models? This is basic survey data--nothing fancy needed. All best, Lisa __ 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] nnet package for OS X
Hi Lisa, as far as I know, nnet is part of the VR bundle which is a recommended package... refer http://www.stats.bris.ac.uk/R/src/contrib/Descriptions/VR.html so just try library(nnet) and see what happens! cheers! Sean On 02/08/05, Lisa Schweitzer [EMAIL PROTECTED] wrote: Hello every one-- My sincere apologies if I have missed something obvious, but I need to construct some multinomial logit models. I remember using nnet for this in the past, but I am having trouble locating a OS X version of this (it's been a year or so since I had to use it). I found a Windows version, but when I try to install across platforms, it crashes R immediately--no error messages or anything helpful. Can any one direct me to a repository, or alternatively, to another package readily available Mac side that can do simple multinomial logit models? This is basic survey data--nothing fancy needed. All best, Lisa __ 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-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] NA when using read.csv
Hi Haibo, Are there any character # in the cells? By default, the things at the right side of the character # will be masked as comments. Wuming On 8/2/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 8/1/05, Haibo Huang [EMAIL PROTECTED] wrote: Hi, when I used: Dist=read.csv(test.csv,header=TRUE) to read data from CSV file. For some cells, R mistakenly put in as NA, while most of the cells still Its not likely that there are errors in this software so its probably not a mistake but some other cause. appears to be right, and there is no error message in R. I am pretty sure the csv file is all right, but just can't figure out what went wrong. Can someone share your thoughts with me? Thanks! Read the posting guide at the bottom of each post. In particular, cut down the input file to a few lines that still exhibit the problem in order to provide a reproducible example. As an aside note that 'header = TRUE' is the default in read.csv (but not on read.table) so you don't need to specify it. __ 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-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] Can R do regression with absolute error loss?
Dear R users, I wonder if there are some functions or packages in R that can perform regression under absolute error loss. I have searched the full manual of R, but can not find it. I really appreciate your help. Best regards, Ben Yang [[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
Re: [R] Can R do regression with absolute error loss?
Check out the quantreg package.
On 8/1/05, Chunyu Yang {msbbb062} [EMAIL PROTECTED] wrote:
Dear R users,
I wonder if there are some functions or packages in R that can perform
regression under absolute error loss. I have searched the full manual of
R, but can not find it. I really appreciate your help.
Best regards,
Ben Yang
[[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
__
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] (no subject)
I've been using R for a while under Mac OS X, which thanks to the R on OS X developers, is probably the best platform for learning R. I recently built a Linux box with a Pentium D processor, and am running an AMD64 port of Ubuntu with the SMP kernel. After setting up the basics on the box, I thought I'd install R as well. I've had a very hard time finding good documentation on some configure flags/options for compiling a 64-bit version of R on Linux. I thought this was very strange, considering Linux is apparently the R developer's platform of choice. Can anyone point me to some good build instructions for 64-bit R under Linux? Should I compile it against any additional libraries for better performance? Thanks in advance for any pointers. --Jake __ 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] linux compile options (64-bit)
Sorry--sent this before I added a subject line... Begin forwarded message: From: Jacob Michaelson [EMAIL PROTECTED] Date: August 1, 2005 8:51:58 PM MDT To: r-help@stat.math.ethz.ch I've been using R for a while under Mac OS X, which thanks to the R on OS X developers, is probably the best platform for learning R. I recently built a Linux box with a Pentium D processor, and am running an AMD64 port of Ubuntu with the SMP kernel. After setting up the basics on the box, I thought I'd install R as well. I've had a very hard time finding good documentation on some configure flags/options for compiling a 64-bit version of R on Linux. I thought this was very strange, considering Linux is apparently the R developer's platform of choice. Can anyone point me to some good build instructions for 64-bit R under Linux? Should I compile it against any additional libraries for better performance? Thanks in advance for any pointers. --Jake [[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
Re: [R] can we manage memory usage to increase speed?
Please refer to the following post. Ed --- Mike Lawrence [EMAIL PROTECTED] wrote: Date: Mon, 1 Aug 2005 00:19:06 -0300 From: Mike Lawrence [EMAIL PROTECTED] To: Briggs, Meredith M [EMAIL PROTECTED] CC: r-help@stat.math.ethz.ch Subject: Re: [R] How do you increase memeory? memory.limit(size = x) where x is the desired memory limit in MB. Quoting Briggs, Meredith M [EMAIL PROTECTED]: Hello Function memory.size() =435109888. How do I increase it by, say 30%? Thanks Meredith [[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 -- Mike Lawrence, BA(Hons) Research Assistant to Dr. Gail Eskes Dalhousie University QEII Health Sciences Centre (Psychiatry) [EMAIL PROTECTED] The road to Wisdom? Well, it's plain and simple to express: Err and err and err again, but less and less and less. - Piet Hein __ 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 --- Zhilin Liu [EMAIL PROTECTED] wrote: Hi, Thanks for reading. I am running a process in R for microarray data analysis. RedHat Enterprise Linux 4, dual AMD CPU, 6G memory. However, the R process use only a total of 200M memory. And the CPU usage is total to ~110% for two. The program takes at least 2 weeks to run at the current speed. Is there some way we can increase the usage of CPUs and memories and speed up? Any suggestion is appreciated. Thanks again. Zhilin [[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 __ 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
