Re: [R] Installing a tar.gz file

2016-03-06 Thread Uwe Ligges



On 06.03.2016 17:49, sowmya s wrote:

Hi:
I am trying to install a tar.gz file and I get the error message,
Warning in install.packages :
   installation of package
‘/Users/subraman/Documents/BCA_Data/TDL_Data/tdllicor_0.1-21.tar.gz’ had
non-zero exit status



We need at least the full output to be able to guess what went wrong.

Best,
Uwe Ligges




I have tried to gunzip the file and use the .tar file to uninstall and that
does not work either. Previous users have suggested that they had no
problems installing the file up until at least R2.15.1.

Could someone please offer some advice ?
Thank you.
Sowmya

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.



__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] stargazer summary statistics by group

2016-03-06 Thread sbihorel

Thanks for the suggestions. I will look into these packages.

On 3/6/2016 9:08 PM, Mitchell Maltenfort wrote:

https://cran.r-project.org/web/packages/tableone/index.html might help

On Sunday, March 6, 2016, Jeff Newmiller > wrote:


Maybe what you really want is the tables package.
--
Sent from my phone. Please excuse my brevity.

On March 6, 2016 4:45:58 PM PST, sbihorel
> wrote:
>Hi Boris,
>
>Sorry, but not really. The example that comes closest is "Flip the
>table
>axes" but this is not right either.
>
>In the design that I need, the year, month, day, etc... variables
would
>
>each get a block of rows with statistics (mean, sd, median, min, max)
>provided for each level of another variable, these levels being
>reported
>in columns.
>
>Basically, I need to look at descriptive statistics that are
>stratified.
>Does this make sense?
>
>Sebastien
>
>On 3/6/2016 4:04 PM, Boris Steipe wrote:
>> Does this entry on the Stargazer Cheatsheet come close enough
to what
>you want?
>>

>http://jakeruss.com/cheatsheets/stargazer.html#the-default-summary-statistics-table
>>
>> B.
>>
>> On Mar 6, 2016, at 3:34 PM, sbihorel
>> wrote:
>>
>>> Hi,
>>>
>>> I saw a post on this topic on stackoverflow a while ago. It
does not
>seem to have got any reply... Just trying my luck here.
>>>
>>> Is there any way to use stargazer to create a table of descriptive
>statistics by group such as the one below?
>>>
>>> Thanks
>>>
>>> #-|-Stat--|--A--|--B--|--Overall-|
>>> # Variable 1  | Stat1 | | |  |
>>> #-|---|-|-|--|
>>> # | Stat2 | | |  |
>>> #-|---|-|-|--|
>>> # | ...   | | |  |
>>> #-|---|-|-|--|
>>> # Variable 2  | Stat1 | | |  |
>>> #-|---|-|-|--|
>>> # | Stat2 | | |  |
>>> #-|---|-|-|--|
>>> # | ...   | | |  |
>>> #-|---|-|-|--|
>>> # Variable ...| Stat1 | | |  |
>>> #-|---|-|-|--|
>>> # | Stat2 | | |  |
>>> #-|---|-|-|--|
>>> # | ...   | | |  |
>>> #-|---|-|-|--|
>>>
>>> stat1 and stat2 could mean, sd, median, etc...
>>>
>>> __
>>> R-help@r-project.org  mailing list -- To
UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>
>--
>Sebastien Bihorel
>Associate Director, Pharmacometrics
>Buffalo Office: +1-716-633-3463 ext. 323 | Website
>
>
>
>__
>R-help@r-project.org  mailing list -- To
UNSUBSCRIBE and more, see
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.

[[alternative HTML version deleted]]

__
R-help@r-project.org  mailing list -- To UNSUBSCRIBE
and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.



--
Sent from Gmail Mobile


--
Sebastien Bihorel
Associate Director, Pharmacometrics
Buffalo Office: +1-716-633-3463 ext. 323 | Website 




__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] ACF values with confidence limits + Plot Extaction

2016-03-06 Thread Jean-Claude Arbaut
a <- as.ts(rnorm(20))
png("acf.png")
a.acf <- acf(a)
dev.off()

# to see what is available
names(a.acf)
unclass(a.acf)


2016-03-06 20:06 GMT+01:00 Preetam Pal :
> Thanks, Michael. Appreciate it.
> But suppose I go for the plot, how to extract it from R ... say, I want to
> save it as a .png file.
> Regards,
> Preetam
>
> On Sun, Mar 6, 2016 at 3:39 PM, Michael Dewey 
> wrote:
>
>> According to the documentation acf returns what you want. It also says
>> that it returns it invisibly if plot = TRUE which I imagine is what you are
>> doing.
>>
>> So try
>>
>> res <- acf(insert_parameters_here, plot = FALSE)
>>
>> and then look at res
>>
>>
>>
>> On 06/03/2016 07:40, Preetam Pal wrote:
>>
>>> Hi R-users,
>>>
>>> I have a time series of residuals and I want to get the ACF
>>> (autocorrelation) values till lag = 12, along with the 12 upper/lower
>>> confidence limits. I understand that acf(residual) would give me the plot,
>>> but I will also need the actual values as an array etc. Plus, I'll have to
>>> extract the plot from R as well. Is there a way to achieve these two?
>>>
>>> Appreciate your help.
>>>
>>> Regards,
>>> Preetam
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> __
>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>>
>> --
>> Michael
>> http://www.dewey.myzen.co.uk/home.html
>>
>
>
>
> --
> Preetam Pal
> (+91)-9432212774
> M-Stat 2nd Year, Room No. N-114
> Statistics Division,   C.V.Raman
> Hall
> Indian Statistical Institute, B.H.O.S.
> Kolkata.
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Installing a tar.gz file

2016-03-06 Thread sowmya s
Hi:
I am trying to install a tar.gz file and I get the error message,
Warning in install.packages :
  installation of package
‘/Users/subraman/Documents/BCA_Data/TDL_Data/tdllicor_0.1-21.tar.gz’ had
non-zero exit status

I have tried to gunzip the file and use the .tar file to uninstall and that
does not work either. Previous users have suggested that they had no
problems installing the file up until at least R2.15.1.

Could someone please offer some advice ?
Thank you.
Sowmya

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] stargazer summary statistics by group

2016-03-06 Thread Mitchell Maltenfort
https://cran.r-project.org/web/packages/tableone/index.html might help

On Sunday, March 6, 2016, Jeff Newmiller  wrote:

> Maybe what you really want is the tables package.
> --
> Sent from my phone. Please excuse my brevity.
>
> On March 6, 2016 4:45:58 PM PST, sbihorel <
> sebastien.biho...@cognigencorp.com > wrote:
> >Hi Boris,
> >
> >Sorry, but not really. The example that comes closest is "Flip the
> >table
> >axes" but this is not right either.
> >
> >In the design that I need, the year, month, day, etc... variables would
> >
> >each get a block of rows with statistics (mean, sd, median, min, max)
> >provided for each level of another variable, these levels being
> >reported
> >in columns.
> >
> >Basically, I need to look at descriptive statistics that are
> >stratified.
> >Does this make sense?
> >
> >Sebastien
> >
> >On 3/6/2016 4:04 PM, Boris Steipe wrote:
> >> Does this entry on the Stargazer Cheatsheet come close enough to what
> >you want?
> >>
> >
> http://jakeruss.com/cheatsheets/stargazer.html#the-default-summary-statistics-table
> >>
> >> B.
> >>
> >> On Mar 6, 2016, at 3:34 PM, sbihorel
> >> wrote:
> >>
> >>> Hi,
> >>>
> >>> I saw a post on this topic on stackoverflow a while ago. It does not
> >seem to have got any reply... Just trying my luck here.
> >>>
> >>> Is there any way to use stargazer to create a table of descriptive
> >statistics by group such as the one below?
> >>>
> >>> Thanks
> >>>
> >>> #-|-Stat--|--A--|--B--|--Overall-|
> >>> # Variable 1  | Stat1 | | |  |
> >>> #-|---|-|-|--|
> >>> # | Stat2 | | |  |
> >>> #-|---|-|-|--|
> >>> # | ...   | | |  |
> >>> #-|---|-|-|--|
> >>> # Variable 2  | Stat1 | | |  |
> >>> #-|---|-|-|--|
> >>> # | Stat2 | | |  |
> >>> #-|---|-|-|--|
> >>> # | ...   | | |  |
> >>> #-|---|-|-|--|
> >>> # Variable ...| Stat1 | | |  |
> >>> #-|---|-|-|--|
> >>> # | Stat2 | | |  |
> >>> #-|---|-|-|--|
> >>> # | ...   | | |  |
> >>> #-|---|-|-|--|
> >>>
> >>> stat1 and stat2 could mean, sd, median, etc...
> >>>
> >>> __
> >>> R-help@r-project.org  mailing list -- To UNSUBSCRIBE
> and more, see
> >>> https://stat.ethz.ch/mailman/listinfo/r-help
> >>> PLEASE do read the posting guide
> >http://www.R-project.org/posting-guide.html
> >>> and provide commented, minimal, self-contained, reproducible code.
> >
> >--
> >Sebastien Bihorel
> >Associate Director, Pharmacometrics
> >Buffalo Office: +1-716-633-3463 ext. 323 | Website
> >
> >
> >
> >__
> >R-help@r-project.org  mailing list -- To UNSUBSCRIBE and
> more, see
> >https://stat.ethz.ch/mailman/listinfo/r-help
> >PLEASE do read the posting guide
> >http://www.R-project.org/posting-guide.html
> >and provide commented, minimal, self-contained, reproducible code.
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org  mailing list -- To UNSUBSCRIBE and
> more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>


-- 
Sent from Gmail Mobile

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] stargazer summary statistics by group

2016-03-06 Thread Jeff Newmiller
Maybe what you really want is the tables package. 
-- 
Sent from my phone. Please excuse my brevity.

On March 6, 2016 4:45:58 PM PST, sbihorel  
wrote:
>Hi Boris,
>
>Sorry, but not really. The example that comes closest is "Flip the
>table 
>axes" but this is not right either.
>
>In the design that I need, the year, month, day, etc... variables would
>
>each get a block of rows with statistics (mean, sd, median, min, max) 
>provided for each level of another variable, these levels being
>reported 
>in columns.
>
>Basically, I need to look at descriptive statistics that are
>stratified. 
>Does this make sense?
>
>Sebastien
>
>On 3/6/2016 4:04 PM, Boris Steipe wrote:
>> Does this entry on the Stargazer Cheatsheet come close enough to what
>you want?
>>
>http://jakeruss.com/cheatsheets/stargazer.html#the-default-summary-statistics-table
>>
>> B.
>>
>> On Mar 6, 2016, at 3:34 PM, sbihorel
> wrote:
>>
>>> Hi,
>>>
>>> I saw a post on this topic on stackoverflow a while ago. It does not
>seem to have got any reply... Just trying my luck here.
>>>
>>> Is there any way to use stargazer to create a table of descriptive
>statistics by group such as the one below?
>>>
>>> Thanks
>>>
>>> #-|-Stat--|--A--|--B--|--Overall-|
>>> # Variable 1  | Stat1 | | |  |
>>> #-|---|-|-|--|
>>> # | Stat2 | | |  |
>>> #-|---|-|-|--|
>>> # | ...   | | |  |
>>> #-|---|-|-|--|
>>> # Variable 2  | Stat1 | | |  |
>>> #-|---|-|-|--|
>>> # | Stat2 | | |  |
>>> #-|---|-|-|--|
>>> # | ...   | | |  |
>>> #-|---|-|-|--|
>>> # Variable ...| Stat1 | | |  |
>>> #-|---|-|-|--|
>>> # | Stat2 | | |  |
>>> #-|---|-|-|--|
>>> # | ...   | | |  |
>>> #-|---|-|-|--|
>>>
>>> stat1 and stat2 could mean, sd, median, etc...
>>>
>>> __
>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>
>-- 
>Sebastien Bihorel
>Associate Director, Pharmacometrics
>Buffalo Office: +1-716-633-3463 ext. 323 | Website 
>
>
>
>__
>R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] stargazer summary statistics by group

2016-03-06 Thread sbihorel

Hi Boris,

Sorry, but not really. The example that comes closest is "Flip the table 
axes" but this is not right either.


In the design that I need, the year, month, day, etc... variables would 
each get a block of rows with statistics (mean, sd, median, min, max) 
provided for each level of another variable, these levels being reported 
in columns.


Basically, I need to look at descriptive statistics that are stratified. 
Does this make sense?


Sebastien

On 3/6/2016 4:04 PM, Boris Steipe wrote:

Does this entry on the Stargazer Cheatsheet come close enough to what you want?

http://jakeruss.com/cheatsheets/stargazer.html#the-default-summary-statistics-table

B.

On Mar 6, 2016, at 3:34 PM, sbihorel  wrote:


Hi,

I saw a post on this topic on stackoverflow a while ago. It does not seem to 
have got any reply... Just trying my luck here.

Is there any way to use stargazer to create a table of descriptive statistics 
by group such as the one below?

Thanks

#-|-Stat--|--A--|--B--|--Overall-|
# Variable 1  | Stat1 | | |  |
#-|---|-|-|--|
# | Stat2 | | |  |
#-|---|-|-|--|
# | ...   | | |  |
#-|---|-|-|--|
# Variable 2  | Stat1 | | |  |
#-|---|-|-|--|
# | Stat2 | | |  |
#-|---|-|-|--|
# | ...   | | |  |
#-|---|-|-|--|
# Variable ...| Stat1 | | |  |
#-|---|-|-|--|
# | Stat2 | | |  |
#-|---|-|-|--|
# | ...   | | |  |
#-|---|-|-|--|

stat1 and stat2 could mean, sd, median, etc...

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


--
Sebastien Bihorel
Associate Director, Pharmacometrics
Buffalo Office: +1-716-633-3463 ext. 323 | Website 




__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] ACF values with confidence limits + Plot Extaction

2016-03-06 Thread Bert Gunter
My word!  Did you try ?png?

More to the point, any good R tutorial should provide such info. Pls make
some reasonable efforts on your own before posting here.

Cheers,
Bert


On Sunday, March 6, 2016, Preetam Pal  wrote:

> Thanks, Michael. Appreciate it.
> But suppose I go for the plot, how to extract it from R ... say, I want to
> save it as a .png file.
> Regards,
> Preetam
>
> On Sun, Mar 6, 2016 at 3:39 PM, Michael Dewey  >
> wrote:
>
> > According to the documentation acf returns what you want. It also says
> > that it returns it invisibly if plot = TRUE which I imagine is what you
> are
> > doing.
> >
> > So try
> >
> > res <- acf(insert_parameters_here, plot = FALSE)
> >
> > and then look at res
> >
> >
> >
> > On 06/03/2016 07:40, Preetam Pal wrote:
> >
> >> Hi R-users,
> >>
> >> I have a time series of residuals and I want to get the ACF
> >> (autocorrelation) values till lag = 12, along with the 12 upper/lower
> >> confidence limits. I understand that acf(residual) would give me the
> plot,
> >> but I will also need the actual values as an array etc. Plus, I'll have
> to
> >> extract the plot from R as well. Is there a way to achieve these two?
> >>
> >> Appreciate your help.
> >>
> >> Regards,
> >> Preetam
> >>
> >> [[alternative HTML version deleted]]
> >>
> >> __
> >> R-help@r-project.org  mailing list -- To UNSUBSCRIBE and
> more, see
> >> https://stat.ethz.ch/mailman/listinfo/r-help
> >> PLEASE do read the posting guide
> >> http://www.R-project.org/posting-guide.html
> >> and provide commented, minimal, self-contained, reproducible code.
> >>
> >>
> > --
> > Michael
> > http://www.dewey.myzen.co.uk/home.html
> >
>
>
>
> --
> Preetam Pal
> (+91)-9432212774
> M-Stat 2nd Year, Room No. N-114
> Statistics Division,   C.V.Raman
> Hall
> Indian Statistical Institute, B.H.O.S.
> Kolkata.
>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org  mailing list -- To UNSUBSCRIBE and
> more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>


-- 
Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] package FD

2016-03-06 Thread Jim Lemon
The values in a$x do look numeric. What do you get from:

class(a$x)

If the result is "factor", as it was for your ft$trait3 variable (and
I hope that a$x is the same variable with a different name), then at
least one of those values must have been read in as non-numeric. The
possible reasons for this are many as Jeff noted and may be a
non-printing character that has crept into your original data file. If
your data set is as small as your example, I would start by loading
the original data file into a hex editor and looking for a character
that shouldn't be there. I once had to write a program in C that
scanned very large files of customer data to find just such
troublemakers and tell me where they were.

Jim

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] stargazer summary statistics by group

2016-03-06 Thread Boris Steipe
Does this entry on the Stargazer Cheatsheet come close enough to what you want?
   
http://jakeruss.com/cheatsheets/stargazer.html#the-default-summary-statistics-table

B.

On Mar 6, 2016, at 3:34 PM, sbihorel  wrote:

> Hi,
> 
> I saw a post on this topic on stackoverflow a while ago. It does not seem to 
> have got any reply... Just trying my luck here.
> 
> Is there any way to use stargazer to create a table of descriptive statistics 
> by group such as the one below?
> 
> Thanks
> 
> #-|-Stat--|--A--|--B--|--Overall-|
> # Variable 1  | Stat1 | | |  |
> #-|---|-|-|--|
> # | Stat2 | | |  |
> #-|---|-|-|--|
> # | ...   | | |  |
> #-|---|-|-|--|
> # Variable 2  | Stat1 | | |  |
> #-|---|-|-|--|
> # | Stat2 | | |  |
> #-|---|-|-|--|
> # | ...   | | |  |
> #-|---|-|-|--|
> # Variable ...| Stat1 | | |  |
> #-|---|-|-|--|
> # | Stat2 | | |  |
> #-|---|-|-|--|
> # | ...   | | |  |
> #-|---|-|-|--|
> 
> stat1 and stat2 could mean, sd, median, etc...
> 
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] stargazer summary statistics by group

2016-03-06 Thread sbihorel

Hi,

I saw a post on this topic on stackoverflow a while ago. It does not 
seem to have got any reply... Just trying my luck here.


Is there any way to use stargazer to create a table of descriptive 
statistics by group such as the one below?


Thanks

#-|-Stat--|--A--|--B--|--Overall-|
# Variable 1  | Stat1 | | |  |
#-|---|-|-|--|
# | Stat2 | | |  |
#-|---|-|-|--|
# | ...   | | |  |
#-|---|-|-|--|
# Variable 2  | Stat1 | | |  |
#-|---|-|-|--|
# | Stat2 | | |  |
#-|---|-|-|--|
# | ...   | | |  |
#-|---|-|-|--|
# Variable ...| Stat1 | | |  |
#-|---|-|-|--|
# | Stat2 | | |  |
#-|---|-|-|--|
# | ...   | | |  |
#-|---|-|-|--|

stat1 and stat2 could mean, sd, median, etc...

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] ACF values with confidence limits + Plot Extaction

2016-03-06 Thread Preetam Pal
Thank you very much,Jean-Claude and Michael.

1> @Michael, your suggestion "plot = FALSE" only returns the estimated ACF
values, not the confidence limits at different lags. May be I am missing
something here.Do you know any way around for this?
2>@Jean-Claude, got it, thanks.

Regards,
Preetam

On Mon, Mar 7, 2016 at 1:03 AM, Jean-Claude Arbaut 
wrote:

> a <- as.ts(rnorm(20))
> png("acf.png")
> a.acf <- acf(a)
> dev.off()
>
> # to see what is available
> names(a.acf)
> unclass(a.acf)
>
>
> 2016-03-06 20:06 GMT+01:00 Preetam Pal :
> > Thanks, Michael. Appreciate it.
> > But suppose I go for the plot, how to extract it from R ... say, I want
> to
> > save it as a .png file.
> > Regards,
> > Preetam
> >
> > On Sun, Mar 6, 2016 at 3:39 PM, Michael Dewey 
> > wrote:
> >
> >> According to the documentation acf returns what you want. It also says
> >> that it returns it invisibly if plot = TRUE which I imagine is what you
> are
> >> doing.
> >>
> >> So try
> >>
> >> res <- acf(insert_parameters_here, plot = FALSE)
> >>
> >> and then look at res
> >>
> >>
> >>
> >> On 06/03/2016 07:40, Preetam Pal wrote:
> >>
> >>> Hi R-users,
> >>>
> >>> I have a time series of residuals and I want to get the ACF
> >>> (autocorrelation) values till lag = 12, along with the 12 upper/lower
> >>> confidence limits. I understand that acf(residual) would give me the
> plot,
> >>> but I will also need the actual values as an array etc. Plus, I'll
> have to
> >>> extract the plot from R as well. Is there a way to achieve these two?
> >>>
> >>> Appreciate your help.
> >>>
> >>> Regards,
> >>> Preetam
> >>>
> >>> [[alternative HTML version deleted]]
> >>>
> >>> __
> >>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> >>> https://stat.ethz.ch/mailman/listinfo/r-help
> >>> PLEASE do read the posting guide
> >>> http://www.R-project.org/posting-guide.html
> >>> and provide commented, minimal, self-contained, reproducible code.
> >>>
> >>>
> >> --
> >> Michael
> >> http://www.dewey.myzen.co.uk/home.html
> >>
> >
> >
> >
> > --
> > Preetam Pal
> > (+91)-9432212774
> > M-Stat 2nd Year, Room No.
> N-114
> > Statistics Division,   C.V.Raman
> > Hall
> > Indian Statistical Institute, B.H.O.S.
> > Kolkata.
> >
> > [[alternative HTML version deleted]]
> >
> > __
> > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
>



-- 
Preetam Pal
(+91)-9432212774
M-Stat 2nd Year, Room No. N-114
Statistics Division,   C.V.Raman
Hall
Indian Statistical Institute, B.H.O.S.
Kolkata.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] ACF values with confidence limits + Plot Extaction

2016-03-06 Thread Preetam Pal
Thanks, Michael. Appreciate it.
But suppose I go for the plot, how to extract it from R ... say, I want to
save it as a .png file.
Regards,
Preetam

On Sun, Mar 6, 2016 at 3:39 PM, Michael Dewey 
wrote:

> According to the documentation acf returns what you want. It also says
> that it returns it invisibly if plot = TRUE which I imagine is what you are
> doing.
>
> So try
>
> res <- acf(insert_parameters_here, plot = FALSE)
>
> and then look at res
>
>
>
> On 06/03/2016 07:40, Preetam Pal wrote:
>
>> Hi R-users,
>>
>> I have a time series of residuals and I want to get the ACF
>> (autocorrelation) values till lag = 12, along with the 12 upper/lower
>> confidence limits. I understand that acf(residual) would give me the plot,
>> but I will also need the actual values as an array etc. Plus, I'll have to
>> extract the plot from R as well. Is there a way to achieve these two?
>>
>> Appreciate your help.
>>
>> Regards,
>> Preetam
>>
>> [[alternative HTML version deleted]]
>>
>> __
>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
> --
> Michael
> http://www.dewey.myzen.co.uk/home.html
>



-- 
Preetam Pal
(+91)-9432212774
M-Stat 2nd Year, Room No. N-114
Statistics Division,   C.V.Raman
Hall
Indian Statistical Institute, B.H.O.S.
Kolkata.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] package FD

2016-03-06 Thread Jeff Newmiller
That code doesn't show what a$x IS, just what comes out after you force things. 
The fact that you felt compelled to apply those functions just makes it seem 
more likely that Jim is onto something. The output of

str( a$x )

would show what kind of data it is, and

dput( a$x )

would let us put the data column into our copy of R and guess about why it 
might not already be numeric, and

dput( a )

followed by your actual code that uses `a` would give us a possibility to 
reproduce your actual error.  At this point the cause of the problem seems 
likely to be in the original input data or the code that you are using to 
import and manipulate that data, which is where most questions like this seem 
to end up. 
-- 
Sent from my phone. Please excuse my br
evity.

On March 6, 2016 3:57:04 AM PST, Fabio Monteiro  
wrote:
>Hey Jim
>
>they are all numeric as you can see
>
>as.numeric(as.character(a$x))
> [1]  20.0  50.0   7.9  25.0  20.0  20.0  15.0  30.0  48.0  75.0  75.0
> 25.0 300.0
>[14] 103.0  20.0  45.0  15.0  20.0  50.0   6.0  18.0  59.0  70.0  80.0
>100.0  40.0
>[27]  15.0  30.0  40.0  60.0   9.0  11.0  27.5  75.0  60.0  70.0  70.0
> 60.0  20.0
>[40]  21.0  50.0  35.0  46.0
>
>Fábio
>
>2016-03-05 9:38 GMT+00:00 Jim Lemon :
>
>> Hi Fabio,
>> What has probably happened is that ft$trait3 looks like numbers but
>> when it was read in at least one value could not be read as a number.
>> The default behavior in R is to transform the variable into a factor:
>>
>> testcase<-read.table(text="1 2 3 4
>>  1 2 3 4
>>  1 2 B 4")
>> > testcase
>>  V1 V2 V3 V4
>> 1  1  2  3  4
>> 2  1  2  3  4
>> 3  1  2  B  4
>> > sapply(testcase,class)
>>   V1V2V3V4
>> "integer" "integer"  "factor" "integer"
>>
>> So testcase$V3 was read in as a factor. There are a couple of things
>> you can do. First, try to convert the factor to numeric and live with
>> the NA values that will be generated:
>>
>> as.numeric(as.character(testcase$V3))
>> [1]  3  3 NA
>> Warning message:
>> NAs introduced by coercion
>>
>> If the number of non0numeric values in the original data is small and
>> you can find them:
>>
>> testcase$V3[3]<-3
>> as.numeric(as.character(testcase$V3))
>> [1] 3 3 3
>>
>> correct the original data and read it in again. I don't really know
>> enough to answer your second question.
>>
>> JIm
>>
>>
>> On Sat, Mar 5, 2016 at 2:58 AM, Fabio Monteiro
>>  wrote:
>> > I still have another question. apart from that one
>> >
>> > in the dbFD function in FD package there is one option which is if
>FRic
>> > should be standardized or not. What does that mean? Why should our
>> shouldn't
>> > I have the FRic Standardized?
>> >
>> > Thank you Jim
>> >
>> > 2016-03-04 14:52 GMT+00:00 Fabio Monteiro
>> >:
>> >>
>> >> class(ft$trait3)
>> >> [1] "factor
>> >>
>> >>
>> >> Yes is a factor. And now?
>> >>
>> >> Thank you
>> >>
>> >> Kind Regards
>> >> Fábio
>> >>
>> >> 2016-03-04 7:15 GMT+00:00 Jim Lemon :
>> >>>
>> >>> Hi Fabio,
>> >>> You should write:
>> >>>
>> >>> class(...)
>> >>>
>> >>> where ... is the same as what you would type to have the variable
>> >>> displayed on the console. Looking at your earlier message, it
>might
>> >>> be:
>> >>>
>> >>> x$trait3
>> >>>
>> >>> so try:
>> >>>
>> >>> class(x$trait3)
>> >>>
>> >>> Jim
>> >>>
>> >>>
>> >>> On Fri, Mar 4, 2016 at 11:30 AM, Fabio Monteiro
>> >>>  wrote:
>> >>> > i just called trait3 to my variable.
>> >>> >
>> >>> > Is this what i'm suppose to wright? class(trait3), or class
>> >>> > (my_trait3_variable?
>> >>> >
>> >>> > both give error
>> >>> >
>> >>> > 2016-03-03 23:42 GMT+00:00 Jim Lemon :
>> >>> >>
>> >>> >> Hi Fabio,
>> >>> >> It is possible that your remaining "numeric" variable is a
>factor.
>> >>> >> What
>> >>> >> does:
>> >>> >>
>> >>> >> class(my_numeric_variable)
>> >>> >>
>> >>> >> say? (where you substitute the name of your "numeric"
>variable)
>> >>> >>
>> >>> >> Jim
>> >>> >>
>> >>> >>
>> >>> >> On Fri, Mar 4, 2016 at 2:25 AM, Fabio Monteiro
>> >>> >>  wrote:
>> >>> >> > Hello, my name is Fábio and I'm a Marine Ecology student in
>> >>> >> > Portugal.
>> >>> >> >
>> >>> >> > I'm currently using the FD package for my work and yesterday
>one
>> >>> >> > message
>> >>> >> > appeared that I wasn't expecting and I really need your help
>to
>> try
>> >>> >> > to
>> >>> >> > figure out what's happening.
>> >>> >> > I'm using the dbFD function and the following message
>appeared:
>> >>> >> >
>> >>> >> > FRic: Only categorical and/or ordinal trait(s) present in
>'x'.
>> FRic
>> >>> >> > was
>> >>> >> >  measured as the number of unique trait combinations, NOT as
>the
>> >>> >> > convex
>> >>> >> > hull volume.
>> >>> >> > FDiv: Cannot be computed when only categorical and/or
>ordinal
>> >>> >> > trait(s)
>> >>> >> > present in 'x'.
>> 

Re: [R] package FD

2016-03-06 Thread Fabio Monteiro
Hey Jim

they are all numeric as you can see

as.numeric(as.character(a$x))
 [1]  20.0  50.0   7.9  25.0  20.0  20.0  15.0  30.0  48.0  75.0  75.0
 25.0 300.0
[14] 103.0  20.0  45.0  15.0  20.0  50.0   6.0  18.0  59.0  70.0  80.0
100.0  40.0
[27]  15.0  30.0  40.0  60.0   9.0  11.0  27.5  75.0  60.0  70.0  70.0
 60.0  20.0
[40]  21.0  50.0  35.0  46.0

Fábio

2016-03-05 9:38 GMT+00:00 Jim Lemon :

> Hi Fabio,
> What has probably happened is that ft$trait3 looks like numbers but
> when it was read in at least one value could not be read as a number.
> The default behavior in R is to transform the variable into a factor:
>
> testcase<-read.table(text="1 2 3 4
>  1 2 3 4
>  1 2 B 4")
> > testcase
>  V1 V2 V3 V4
> 1  1  2  3  4
> 2  1  2  3  4
> 3  1  2  B  4
> > sapply(testcase,class)
>   V1V2V3V4
> "integer" "integer"  "factor" "integer"
>
> So testcase$V3 was read in as a factor. There are a couple of things
> you can do. First, try to convert the factor to numeric and live with
> the NA values that will be generated:
>
> as.numeric(as.character(testcase$V3))
> [1]  3  3 NA
> Warning message:
> NAs introduced by coercion
>
> If the number of non0numeric values in the original data is small and
> you can find them:
>
> testcase$V3[3]<-3
> as.numeric(as.character(testcase$V3))
> [1] 3 3 3
>
> correct the original data and read it in again. I don't really know
> enough to answer your second question.
>
> JIm
>
>
> On Sat, Mar 5, 2016 at 2:58 AM, Fabio Monteiro
>  wrote:
> > I still have another question. apart from that one
> >
> > in the dbFD function in FD package there is one option which is if FRic
> > should be standardized or not. What does that mean? Why should our
> shouldn't
> > I have the FRic Standardized?
> >
> > Thank you Jim
> >
> > 2016-03-04 14:52 GMT+00:00 Fabio Monteiro  >:
> >>
> >> class(ft$trait3)
> >> [1] "factor
> >>
> >>
> >> Yes is a factor. And now?
> >>
> >> Thank you
> >>
> >> Kind Regards
> >> Fábio
> >>
> >> 2016-03-04 7:15 GMT+00:00 Jim Lemon :
> >>>
> >>> Hi Fabio,
> >>> You should write:
> >>>
> >>> class(...)
> >>>
> >>> where ... is the same as what you would type to have the variable
> >>> displayed on the console. Looking at your earlier message, it might
> >>> be:
> >>>
> >>> x$trait3
> >>>
> >>> so try:
> >>>
> >>> class(x$trait3)
> >>>
> >>> Jim
> >>>
> >>>
> >>> On Fri, Mar 4, 2016 at 11:30 AM, Fabio Monteiro
> >>>  wrote:
> >>> > i just called trait3 to my variable.
> >>> >
> >>> > Is this what i'm suppose to wright? class(trait3), or class
> >>> > (my_trait3_variable?
> >>> >
> >>> > both give error
> >>> >
> >>> > 2016-03-03 23:42 GMT+00:00 Jim Lemon :
> >>> >>
> >>> >> Hi Fabio,
> >>> >> It is possible that your remaining "numeric" variable is a factor.
> >>> >> What
> >>> >> does:
> >>> >>
> >>> >> class(my_numeric_variable)
> >>> >>
> >>> >> say? (where you substitute the name of your "numeric" variable)
> >>> >>
> >>> >> Jim
> >>> >>
> >>> >>
> >>> >> On Fri, Mar 4, 2016 at 2:25 AM, Fabio Monteiro
> >>> >>  wrote:
> >>> >> > Hello, my name is Fábio and I'm a Marine Ecology student in
> >>> >> > Portugal.
> >>> >> >
> >>> >> > I'm currently using the FD package for my work and yesterday one
> >>> >> > message
> >>> >> > appeared that I wasn't expecting and I really need your help to
> try
> >>> >> > to
> >>> >> > figure out what's happening.
> >>> >> > I'm using the dbFD function and the following message appeared:
> >>> >> >
> >>> >> > FRic: Only categorical and/or ordinal trait(s) present in 'x'.
> FRic
> >>> >> > was
> >>> >> >  measured as the number of unique trait combinations, NOT as the
> >>> >> > convex
> >>> >> > hull volume.
> >>> >> > FDiv: Cannot be computed when only categorical and/or ordinal
> >>> >> > trait(s)
> >>> >> > present in 'x'.
> >>> >> >
> >>> >> > My data:
> >>> >> > x is a matrix with species vs functional traits
> >>> >> > a is a matrix with species vs sampling (in abundances)
> >>> >> >
> >>> >> > Previously I used the dbFD function and was working just fine.
> >>> >> > Yesterday
> >>> >> > I
> >>> >> > removed 2 traits and this message appeared.
> >>> >> >
> >>> >> > My traits now are 3 categorical traits and 1 numeric. The 2 trais
> >>> >> > that I
> >>> >> > removed were numeric traits as well. I really need to remove those
> >>> >> > trait,
> >>> >> > but I still need the FDiv to be calculated. Can you explain to me
> >>> >> > why is
> >>> >> > this error occurring? I need to know how the dbFD is measuring the
> >>> >> > indexes
> >>> >> > so I can understanding the error and if I can or can't continue to
> >>> >> > use
> >>> >> > this
> >>> >> > package (if it applies or not to my goals)
> >>> >> >
> >>> >> > Kind regards
> >>> >> >
> >>> >> > Fábio Monteiro
> >>> >> >
> >>> >> > [[alternative HTML version 

Re: [R-es] usar R a traves de la web

2016-03-06 Thread eric
Hola Jesus, muchas gracias, esta opcion se ve bastante bien, aunque no 
encuentro un modo de subir mis archivos de datos, pero debe poderse, no ?


Muchas gracias,

Eric.




On 03/05/2016 09:03 AM, Jesús Para Fernández wrote:

Buscas algo asi?

http://www.r-fiddle.org/#/



 > From: ericconchamu...@gmail.com
 > To: r-help-es@r-project.org
 > Date: Sat, 5 Mar 2016 09:01:42 -0300
 > Subject: [R-es] usar R a traves de la web
 >
 > Estimada comunidad, para mi trabajo uso latex y R normalmente, ahora
 > debo viajar sin mi portatil, pero tengo la opcion de llevar un pequeño
 > tablet (con android) ... para suplir latex he estado usando
 > www.overleaf.com y trabaja excelente, practicamente todos los paquetes
 > que uso estan disponibles ahi ... pero no he encontrado algo similar
para R.
 >
 > Saben ustedes si existe algun proyecto que permita usar R en la web ? ya
 > he mirado muchas opciones
 >
(https://nsaunders.wordpress.com/2009/11/30/a-brief-survey-of-r-web-interfaces/)

 > como las de Rstudio (tengo android, no puedo instalar R en el tablet),
 > tampoco tengo la opcion de instalar un RApache por ejemplo, pero ninguna
 > me sirve ... necesito una pagina en que simplemente entre, pegue el
 > codigo y me devuelva las salidas sin tener que instalar nada localmente.
 >
 > Muchas gracias,
 >
 > Eric.
 >
 >
 >
 > --
 > Forest Engineer
 > Master in Environmental and Natural Resource Economics
 > Ph.D. student in Sciences of Natural Resources at La Frontera University
 > Member in AguaDeTemu2030, citizen movement for Temuco with green city
 > standards for living
 >
 > Nota: Las tildes se han omitido para asegurar compatibilidad con algunos
 > lectores de correo.
 >
 > ___
 > R-help-es mailing list
 > R-help-es@r-project.org
 > https://stat.ethz.ch/mailman/listinfo/r-help-es


--
Forest Engineer
Master in Environmental and Natural Resource Economics
Ph.D. student in Sciences of Natural Resources at La Frontera University
Member in AguaDeTemu2030, citizen movement for Temuco with green city 
standards for living


Nota: Las tildes se han omitido para asegurar compatibilidad con algunos 
lectores de correo.


___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es


Re: [R] ACF values with confidence limits + Plot Extaction

2016-03-06 Thread Michael Dewey
According to the documentation acf returns what you want. It also says 
that it returns it invisibly if plot = TRUE which I imagine is what you 
are doing.


So try

res <- acf(insert_parameters_here, plot = FALSE)

and then look at res


On 06/03/2016 07:40, Preetam Pal wrote:

Hi R-users,

I have a time series of residuals and I want to get the ACF
(autocorrelation) values till lag = 12, along with the 12 upper/lower
confidence limits. I understand that acf(residual) would give me the plot,
but I will also need the actual values as an array etc. Plus, I'll have to
extract the plot from R as well. Is there a way to achieve these two?

Appreciate your help.

Regards,
Preetam

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.



--
Michael
http://www.dewey.myzen.co.uk/home.html

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.