regards,
Katharina.
--
Dr Katharina Fritsch B.Sc. M.Sc. MRSC
Chemical Modeller, Chemical and Process Modelling
E.
katharina.frit...@nnl.co.uk
T.
+44 (0)1925 289387
@uknnl
National Nuclear Laboratory Limited, 5th Floor, Chadwick House,
Birchwood Park, Warrington, WA3 6AE, UK
www.nnl.co.uk
From
Hiya,
I have lost track of half the discussion so forgive me if that suggestion has
been made before, but would it be suitable to wrap your code in a shiny app and
disseminate it this way?
Kind regards,
Katharina.
--
Dr Katharina Fritsch B.Sc. M.Sc. MRSC
Chemical Modeller, Chemical and Process
list said
it would be alright' is going to convince your IT department of the
harmlessness of R?
Cheers,
Katharina.
--
Dr Katharina Fritsch B.Sc. M.Sc. MRSC
Chemical Modeller, Chemical and Process Modelling
E.
katharina.frit...@nnl.co.uk
T.
+44 (0)1925 289387
@uknnl
National Nuclear
after updating R from 3.3.3. to 3.4.0 i cannot import spss-data files
anymore. for the european social survey (europeansocialsurvey.org) i get
this warning:
re-encoding from CP1252
Fehler in levels<-(*tmp*, value = if (nl == nL) as.character(labels)
else paste0(labels, :
factor level [3] is du
Hi,
I want to obtain post hoc comparisons for a model with a three-way interaction
(FullModel5 <- lmer(CIVADmc ~ 1 + Factor1 * Factor2 * Factor3 + (1|ID) + (1 |
Item), data=Exp3_WTP_diffi).
I tried: "summary(glht(FullModel4, linfct = mcp(Factor1*Factor2*Factor3 =
"Tukey")), test = adjusted("holm
i came across the GREA add-in for RStudio on R-bloggers:
http://feedproxy.google.com/~r/RBloggers/~3/dU-3B7fQJ2c/?utm_source=feedburner&utm_medium=email
https://github.com/Stan125/GREA
and wanted to try it straight away - yet I failed already with the installation.
after typing in
> devtools::inst
Dear Community,
Just a short and simple question, but the code does not come to my mind.
I want to plot the Resids by quarterly Data. But there are many quarters.
Just rotate them, does not make the plot clearer/ easier to read.
Is there a possibility to just show every 4th value on the x-axi
Dear Community,
unfortunately I can´t give you an reproducable example, because I really do
not understand why this messages pops up.
I estimate an Fixed Effects Modell, controlling for HAC, because F-statistic
changes, I want to compute it, for the other model-specifications it works,
But for t
-Ursprüngliche Nachricht-
Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im
Auftrag von Katharina Mersmann
Gesendet: Dienstag, 1. Juli 2014 16:22
An: r-help@r-project.org
Betreff: [R] (PLM- package) Residual-Plotting and missing Values
Dear R-Community,
I tried
Dear R-Community,
I tried plotting the residuals of an FE-model estimated via plm .
And detected that there are no residuals in the plot for the last two
countries.
I guess this happens because for some countries values are missing and R
gives me the following for
> fixed.reg1.new$resid[1]
Dear Community,
I emphasize the use of graphical methods to examine residuals for a Panel
model and want to plot the Standard Deviation of residuals by group.
I used the following for plotting residuals by group:
>library(lattice)
>xyplot((residuals(fixed.reg1.1))~countrynr,data=data.plm
Hey R Users!
I´m trying to run a waldtest comparing an restricted and unrestricted model.
But sample size is different (because of NA´s), so I get an error message:
> caus2<-plm(CSmean~ lag(CSmean,1)+ lag(numfull_FCRlong,1)+
+GDP+linpol_GDPpCapita+fort_budget+fort_percentDeb
Hello R-Community,
I actually don´t know how to implement the following:
I have quartaly Panel- Data and want to run a Granger- Test. I´ve already
done a univariate autoregression with grangertest() from package lmtest.
Because the Granger test is designed to handle pairs of variables
Dear R-community,
I am totally lost and need help.
For a Visualization I need to plot two regressions in the same plot. The
intention is to provide a visual basis for a synthesized theoretical
discussion.
1.)fixed.1<-plm(CSmean~ FCRlong, data = data.plm, index =
c("countrynr","quartal")
Hey arun,
I am thankful for your help and assistance! With the example it works!
Have a nice day,
Katie
-Ursprüngliche Nachricht-
Von: arun [mailto:smartpink...@yahoo.com]
Gesendet: Freitag, 25. April 2014 18:01
An: r-help@r-project.org
Cc: Katharina Mersmann
Betreff: Re: [R] purtest
NA/NaN/Inf in 'x'
Thanks for your hints and suggestions!
Katie
-Ursprüngliche Nachricht-
Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im
Auftrag von Katharina Mersmann
Gesendet: Freitag, 25. April 2014 15:53
An: r-help@r-project.org
Betreff: [R]
Hello R-Community,
I have quarterly panel-data and not surprisingly missing values, in a few
variables, which are differently distributed around the panel.
Now I want to run different unit-root tests.
For ADF on the pooled data set, I chose CADF-package, which can determine the
number of
This post was nearly what I was searching for.
Im actually trying to reproduce my Stata results in R and don´t overcome the
problem of the NeweyWest Estimators. I have quarterly PanelData
In Stata i used:
newey y x, lag(4) force
In R this should be
reg1.2<-lm(y~x)
coeftest(reg1.2,
Dear R-Users,
I am fairly new to R and got in trouble by understanding how to run a GAM
using penalized regression splines with 4 degress of freedom (even by
reading the R Documentation).
I tried the following:
> gamreg1.2<-gam(num_FCRlong ~ s(GDP,df=4)+s(cupol_GDPpCapita,df=4)
+
Hello R-Users,
I have a Problem in understanding my results using the davies.test
I read the pdf-documentation on the segmented package and followed the
instructions.
1st performing the Davies.test :
> davies.test(reg1.2,~lagBYmean)
Davies' test for a change in the slop
Hi Simon,
many thanks for looking into this and making me understand the problem!
I'll adjust my factor levels right away...
Best, Katharina
On 3 February 2014 12:42, Simon Wood wrote:
> Hi Katharina,
>
> Thanks for sending this.
>
> The problem is that the prediction d
Hi Simon,
thank you for your reply, I really appreciate any help to understand
the problem here...
Unluckily the package upgrade didn't help with this issue.
An example reproducing the error, and a current sessionInfo() Output
can be found below.
Many thanks once again,
Katharin
lem or mistake within my model
declaration/approach?
I attached the sessionInfo() Output below and the xylemRohWeekXnn2011
dump can be downloaded here:
https://webdisk.ads.mwn.de/Handlers/AnonymousDownload.ashx?folder=1a7cbaa4&path=xylemRohWeekXnn2011.txt
I`m appreciating any
many thanks for the enlightenment - I guess I completely misunderstood the
na.rm option here...
On 18 June 2013 00:50, Gabor Grothendieck wrote:
> On Mon, Jun 17, 2013 at 6:23 PM, Katharina May
> wrote:
> > Dear R Users,
> >
> > I've got a strange problem, whi
13.00 14.33 15.00 15.00
Am I missing something here?
Any help is very much appreciated...
Thanks,
Katharina
###
Background Info:
###
sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-w64-mingw32/x64 (64-bit)
locale:
[1] LC_COLLAT
I am investigating various R packages that facilitate estimation of linear
Gaussian multivariate state space models. I stumbled across the MARSS
package (http://cran.r-project.org/web/packages/MARSS/index.html), which I
believe is very well done, but am finding one missing feature that I cannot
liv
On Jun 7, 2012, at 12:55 , Gabor Grothendieck wrote:
> On Wed, Jun 6, 2012 at 4:40 PM, May Katharina
> wrote:
>> On Jun 6, 2012, at 10:19 , Gabor Grothendieck wrote:
>>
>>> On Wed, Jun 6, 2012 at 3:55 PM, May Katharina
>>> wrote:
>>>> Hello,
>
On Jun 6, 2012, at 10:19 , Gabor Grothendieck wrote:
> On Wed, Jun 6, 2012 at 3:55 PM, May Katharina
> wrote:
>> Hello,
>>
>> I'm trying to use na.spline (package zoo) to fill some missing data in a
>> time series.
>> this works fine, however, if I appl
ion I'm new to time
series & zoo and that I'm terrible sorry if this is
an absolute dull question).
The data I use is rather large, but the error can be reconstructed in a small
example (s. below).
If of any interest: I use R 2.15.0 GUI 1.51 Leopard build 32-bit (6148).
Many thanks
%d.%m.%Y %H:%M"), 10 * 60)
which uses the align.time function from the library xts BUT as far as I can
see only supports either all timestamps being either rounded up or down
(not both, depending on the minutes 0-9 being either > or < 5 line the way
I
fine factor levels as na directly in an svydesign-object?
best and thanks
katharina
dr. katharina manderscheid
soziologisches seminar
universität luzern
kasernenplatz 3
6000 luzern 7
tel. ++41 41 228 4657
Von: Thomas Lumley [tlum...@uw.edu]
Gesendet: Don
i defined the na's as such.
any ideas???
thanks a lot,
katharina
--
dr. katharina manderscheid
soziologisches seminar
universität luzern
kasernenplatz 3
6000 luzern 7
tel. ++41 41 228 4657
web: http://www.unilu.ch/deu/dr.-kath
Thank you very much everybody for your quick reply. My problem was that I
could not install packages.
Solution of the problem: It was a Windows 7 problem. I can install the
packages if I run R as administrator.
Thanks again
Katharina
[[alternative HTML version deleted
n")
Warnung in install.packages("foreign") :
'lib = "C:/PROGRA~1/R/R-212~1.1/library" ist nicht schreibbar
Fehler in install.packages("foreign") : kann keinee Pakete installieren
> install.packages("Rcmdr")
Warnung in install.pac
Hi,
I am trying to manipulate a gls regression model output to adjust for use of
two-stage least squares. Basically, I want to estimate a model, then feed in
a new set of residuals, then re-calculate all of the model output (i.e. the
standard errors of the estimators, etc.). I have found some docu
I am using the gls function to fit a two-stage least squares model with
first order autoregressive error terms. Since there is no automated
adjustment for the use of two-stage least squares in this package, I am
trying to manually replicate standard errors of the coefficient estimates in
order to a
hi all,
maybe my last post was to short and cryptic.
i conducted a correspondence analysis with ca(my.data) and would like to have
the x and y axis labelled with the eigenvalues of the dimensions? the argument
xlab/ylab does not work...
thanks for suggestions,
kat
hi all,
how can i get correspondence analysis plots (package ca, plot.ca()) with the
inertia on the two dimensions of the map?
thanks,
kat
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
hi,
i have a spss portable file which i cannot open in R.
i tried
library(Hmisc)
File.spss <- spss.get("File.por", use.value.labels=T)
and received:
Fehler in read.spss("File.por", :
Datei »File.por« ist nicht in einem unterstützten SPSS-Format
(in engl.: error, not supported spss format)
hi lexi
indeed, for a week or so i can also plot the diagnostics - same script, no
change in data... maybe an update of a package, i don't know.
katharina
--
dr. katharina manderscheid
soziologisches seminar
universität luzern
kasernenplatz 3
An: Manderscheid Katharina
Cc: r-help@r-project.org
Betreff: Re: [R] package ca/rgl for ubuntu?
On 06/12/2010 1:45 PM, Manderscheid Katharina wrote:
> hi,
> for some unknown reason i cannot install the package ca in R running in a
> ubuntu mint system. i keep getting the following erro
7; hatte Exit-Status ungleich 0
2: In install.packages("ca") :
Installation des Pakets 'ca' hatte Exit-Status ungleich 0
thanks a lot for your help!
katharina
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
=25, C=50, B=27,
A=100),
while also leaving the NA in the dataframe. Can I transform the whole dataframe
at once or do I need to transform each column separately? How do I correctly
specify the levels and labels of the new dataframe?
Thank you so much,
Katharina
[[alternative HTML
quot;M"] <- "1841-2250"
income[hinctnta=="F"] <- "2251-2710"
income[hinctnta=="S"] <- "2711-3190"
income[hinctnta=="K"] <- "3191-3770"
income[hinctnta=="P"] <- "3771-4550"
income[hinctnta==&quo
:
> reg.satlife <- lm(satlife ~ as.numeric(income)+ freunde + gesund+ satwirt)
the problem appears when calling
> par(mfrow=c(2,2))
> plot(reg.satlife)
the problem persists also when running the model without the factors income and
gesund.
do you need more info
this problem seems to only exist in R 2.12.0 but not in R 2.11.1.
any ideas? a bug?
--
dr. katharina manderscheid
soziologisches seminar
universität luzern
kasernenplatz 3
6000 luzern 7
tel. ++41 41 228 4657
web: http://www.unilu.ch/deu/dr
working with
an attached dataset.
can anyone help, please??
thanks a lot,
kat
--
dr. katharina manderscheid
soziologisches seminar
universität luzern
kasernenplatz 3
6000 luzern 7
tel. ++41 41 228 4657
web: http://www.unilu.ch/deu/dr.-katharina-mander
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
I have imported a spss-data set with the read.spss command and received the
following warning
messages. Does this mean that something did not work correctly? What am I to
make out of these messages?
Thanks for your help, Katharina
Warning messages
I have imported a spss-data set with the read.spss command and received the
following warning
messages. Does this mean that something did not work correctly? What am I to
make out of these messages?
Thanks for your help, Katharina
Warning messages:
1: In read.spss("/Daten/ESS4DE.spss/ESS4D
doing an abroad internship, so
my Spanish is far from perfect and my tutorial is in English.
Is there any way that I´m not aware of to change the language environment of R?
Thanks a lot!
Katharina
Hi Andreas, thank you for replying so soon. You are right, your code works
perfectly for my example. Unfortunately, the real data contains inner NAs,
which I forgot to include. Sorry about that. Any other suggestions?
Thanks,
Katharina
Original-Nachricht
> Datum: Wed, 12
ing values were
replaced by NA.
Does that already impede creating a ts-object?
Is there another possibility to perform the cross-correlation?
Thanks in advance,
Katharina
--
___
Katharina Appel
Große Weinmeisterstraße 30
14469 Potsdam
Tel.: 0331-2370828
e-mail: kaj...@gmx.de
ubling you: I struggling with the data analysis of my
bachelor thesis and
just want to compare the coefficients of my linear regression to
published ones...
a method often used I guess, but still I cannot find any appropriate
documentations.
Thanks,
Katharina
> On Aug 8, 2009, Kathari
nderstanding what I have
to provide as arguments.
I would be more than happy if somebody can point out an example similar to the
comparison I have to do...
Thanks,
Katharina
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
If anybody might know how to calculate the standard error of the coefficients
(slope, intercept)
of a lmodel2 that would help me a lot...
katharina wrote:
>
> sorry for spaming, but IU just had an idea not sure if that may be a way
> of
> doing it:
>
> 1. calculate the sta
2009/7/20 Katharina May
> Hi *,
>
> is there a way to obtain some kind of p-value for a model fitted with RMA
> using the lmodel2 package?
> I know that p-values are discussed and criticized a lot and as you can
> image from my question I'm not
> very much of a statisti
ut of lm?
Sorry for bothering everybody with this, well, probably rather idiotic
question, but I don't know where to
continue from this point...
Thanks,
Katharina
Here the output of my lmodel2 regression:
Model II regression
Call: lmodel2(formula = log(AGB) ~ log(BM_roots), data =
think more creatively about how to make a
> more meaningful display to provide viewers interpretable information.
>
> Bert Gunter
> Genentech Nonclinical Biostatistics
>
>
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] O
ing to the par {graphics} help page customizable to up
to 8 different line types).
Does anybody have any idea how these 26 different lines can be made
distinguishable from each other without the use
of colors?
Thanks,
Katharina
__
R-he
thanks Deepayan, that works great!
2009/6/19 Deepayan Sarkar :
> On 6/18/09, Katharina May wrote:
>> Hi there,
>>
>> sorry for troubling everybody once again, I've got a problem rewriting
>> Sarkar's function for
>> rewriting the tick locations
,
Katharina
2009/6/18 Jim Lemon :
> Katharina May wrote:
>>
>> What I sort of looking for is a xyplot with both the x axis at the bottom
>> and y axis at the
>> left going through 0, but continuing in the positive and negative area
>> (forming a kind
em is that if I replace
2 with "e" (using paste()), I get the error message that the location
isn't a numeric value.
Is there any way to get this working somehow or do I have to take a
different approach?
Thanks, Katharina
Here my failing approach:
require(lattice)
data(Ear
http://upload.wikimedia.org/wikipedia/commons/thumb/2/2c/Erf_plot.svg/600px-Erf_plot.svg.png
I'm not sure if this possible at all or if I just have to add to reference
lines at x=0 and y =0 as a
workaround even though not fully satisfying...?
Thanks,
Katharina
2009/6/18
>
> Hi Katharina,
&
=0,
x-axis at the bottom through y=0).
And is there a way to have the top axis and right axis not drawn at all? I
managed to get rid of the ticks, but not any further...
Hope somebody can help me struggling with lattice,
Thanks,
Katharina
katharina wrote:
>
> Hi there,
>
to draw the
right and top axes at all an how to align the y axis
that it goes through x=0?
Thanks,
Katharina
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.o
he variables for which I group my
lattice plots. The color variable is of type character so this
shouldn't be the reason)
The object sub_agb_object can be downloaded here:
http://www.wzw.tum.de/oekophys/fileadmin/sub_agb_data.rData
Maybe somebody has any i
uot; "yellow" "violet"
> "black" "violet" "violet" "black"
> [12] "yellow" "violet" "yellow" "violet" "violet" "violet" "yellow"
> "violet" &
sorry, the attachment can be found here for download:
http://www.wzw.tum.de/waldinventur/fileadmin/coefficient2.RData
katharina wrote:
>
> Hi there,
>
> I trying to solve this problem for the whole day not going anywhere,
> so I really hope maybe somebody can help
> me
Maybe someone can explain me the color issue and the problem with the
NA values which results in a wrong labeling and to few
points being displayed? I'm new to R as you can guess and my code
isn't really elegant but I really cannot get faults within it...
Attached you can find the re
thanks, now I know what to look for
Sorry the missunderstanding
-Katharina
David Winsemius wrote:
>
>
> On Jun 1, 2009, at 2:27 PM, katharina wrote:
>
>>
>> I've calculated coefficient$intercept for each site I took the data
>> from.
>> I want to
graph/plot?
I hope I made myself more clear this time...?
Thanks, Katharina
jholtman wrote:
>
> ?text
>
> On Mon, Jun 1, 2009 at 11:07 AM, Katharina May
> > wrote:
>
>> Hi everyone,
>>
>> is it possible to write a certain value beside a poin
needles per site [kg]")
and I want to have the value of coefficient$site_no written beside
each point within the plot (to recognize them)
sorry for troubling you with this question, but I somehow can't find
any information on it
(maybe because I'm not quite sure what to se
r i:
select * from [biomass_data$] where site_no in (2, 4, 6, 7)
Here my approach (which doesn't work):
site_all_data <= sqlQuery(channel, "select * from [biomass_data$]
where site_no in (",paste(i,sep=","),") )
sorry for spaming so much today to the mailin
thanks to all your solutions, works out perfectly!
2009/5/19 Henrique Dallazuanna :
> Try this:
>
> rm(list=ls(patt="site[0-9]$"))
>
> On Tue, May 19, 2009 at 7:47 AM, Katharina May
> wrote:
>>
>> Hi,
>>
>> how can I use rm() on objects
t))[1]) ||
dim((get(object)))[1] == 0) rm(list=object)
Thanks a lot!
2009/5/19 Jim Lemon :
> Katharina May wrote:
>>
>> Hi,
>>
>> how can I remove all empty objects (which are NA or have zero rows)
>> from my workspace?
>>
>>
>
> Hi Katharina,
> T
thought paste() would create
something like that...?
Thanks,
Katharina
--
Time flies like an arrow, fruit flies like bananas.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting g
Hi,
how can I remove all empty objects (which are NA or have zero rows)
from my workspace?
Thanks,
Katharina
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
nnel,
paste("select * from [biomass_data$] where site_no = ",i,sep="")))
}
Thanks for all your solution proposals!
Katharina
On 21.04.2009, at 02:43, David Winsemius wrote:
On Apr 20, 2009, at 8:04 PM, Katharina May wrote:
Hi there,
I've got a database or rath
[biomass_data$] where site_no = ",i,sep="")))
}
I would be very, very glad if anybody sends me a clue about this, as
I'm a obviuos Newbie using
R...
Thanks,
Katharina
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
79 matches
Mail list logo
