Get full P value in R when it's less than 2.225074e-308 (.Machine&double.xmin)











up vote
0
down vote

favorite












I have a question in regards to getting very very tiny p-value to print fully in R.



I ran a survival analysis using the code:



x <- coxph(Surv(my_data$Start,my_data$time_variable,my_data$event) ~ 
(my_data$variable1))
s <- summary(x)
print(s)


In the summary results as the P value was so small it just displayed as "< 2e-16"



Normally, if this happens, I just print out the p-value separately to get the actual number using the below command:



coef(summary(x))[1,5]


However, this time it is just printing "0", I think this is because the p-value is less than ".Machine$double.xmin" which is 2.225074e-308.
I have also tried changing this to see if I can print a smaller number but if I change .Machine$double.xmin to less than 2.225074e-308 it just automatically changes the number to 0.



Is there a way to get the actual number for a p-value when it's so small or is this a limitation of R that cannot be gotten around?



Any help would be greatly appreciated! Many thanks!










share|improve this question






















  • Why do you need the exact value of such a low p-value? The only way to calculate it is to work on the logarithmic scale, which is supported by R's distribution functions. Likely one of these functions is used internally to calculate the p-value here.
    – Roland
    Nov 21 at 16:46










  • I'm jealous this is a problem I never have... :)
    – Ian Wesley
    Nov 21 at 17:10















up vote
0
down vote

favorite












I have a question in regards to getting very very tiny p-value to print fully in R.



I ran a survival analysis using the code:



x <- coxph(Surv(my_data$Start,my_data$time_variable,my_data$event) ~ 
(my_data$variable1))
s <- summary(x)
print(s)


In the summary results as the P value was so small it just displayed as "< 2e-16"



Normally, if this happens, I just print out the p-value separately to get the actual number using the below command:



coef(summary(x))[1,5]


However, this time it is just printing "0", I think this is because the p-value is less than ".Machine$double.xmin" which is 2.225074e-308.
I have also tried changing this to see if I can print a smaller number but if I change .Machine$double.xmin to less than 2.225074e-308 it just automatically changes the number to 0.



Is there a way to get the actual number for a p-value when it's so small or is this a limitation of R that cannot be gotten around?



Any help would be greatly appreciated! Many thanks!










share|improve this question






















  • Why do you need the exact value of such a low p-value? The only way to calculate it is to work on the logarithmic scale, which is supported by R's distribution functions. Likely one of these functions is used internally to calculate the p-value here.
    – Roland
    Nov 21 at 16:46










  • I'm jealous this is a problem I never have... :)
    – Ian Wesley
    Nov 21 at 17:10













up vote
0
down vote

favorite









up vote
0
down vote

favorite











I have a question in regards to getting very very tiny p-value to print fully in R.



I ran a survival analysis using the code:



x <- coxph(Surv(my_data$Start,my_data$time_variable,my_data$event) ~ 
(my_data$variable1))
s <- summary(x)
print(s)


In the summary results as the P value was so small it just displayed as "< 2e-16"



Normally, if this happens, I just print out the p-value separately to get the actual number using the below command:



coef(summary(x))[1,5]


However, this time it is just printing "0", I think this is because the p-value is less than ".Machine$double.xmin" which is 2.225074e-308.
I have also tried changing this to see if I can print a smaller number but if I change .Machine$double.xmin to less than 2.225074e-308 it just automatically changes the number to 0.



Is there a way to get the actual number for a p-value when it's so small or is this a limitation of R that cannot be gotten around?



Any help would be greatly appreciated! Many thanks!










share|improve this question













I have a question in regards to getting very very tiny p-value to print fully in R.



I ran a survival analysis using the code:



x <- coxph(Surv(my_data$Start,my_data$time_variable,my_data$event) ~ 
(my_data$variable1))
s <- summary(x)
print(s)


In the summary results as the P value was so small it just displayed as "< 2e-16"



Normally, if this happens, I just print out the p-value separately to get the actual number using the below command:



coef(summary(x))[1,5]


However, this time it is just printing "0", I think this is because the p-value is less than ".Machine$double.xmin" which is 2.225074e-308.
I have also tried changing this to see if I can print a smaller number but if I change .Machine$double.xmin to less than 2.225074e-308 it just automatically changes the number to 0.



Is there a way to get the actual number for a p-value when it's so small or is this a limitation of R that cannot be gotten around?



Any help would be greatly appreciated! Many thanks!







r p-value






share|improve this question













share|improve this question











share|improve this question




share|improve this question










asked Nov 21 at 14:31









Caragh

255




255












  • Why do you need the exact value of such a low p-value? The only way to calculate it is to work on the logarithmic scale, which is supported by R's distribution functions. Likely one of these functions is used internally to calculate the p-value here.
    – Roland
    Nov 21 at 16:46










  • I'm jealous this is a problem I never have... :)
    – Ian Wesley
    Nov 21 at 17:10


















  • Why do you need the exact value of such a low p-value? The only way to calculate it is to work on the logarithmic scale, which is supported by R's distribution functions. Likely one of these functions is used internally to calculate the p-value here.
    – Roland
    Nov 21 at 16:46










  • I'm jealous this is a problem I never have... :)
    – Ian Wesley
    Nov 21 at 17:10
















Why do you need the exact value of such a low p-value? The only way to calculate it is to work on the logarithmic scale, which is supported by R's distribution functions. Likely one of these functions is used internally to calculate the p-value here.
– Roland
Nov 21 at 16:46




Why do you need the exact value of such a low p-value? The only way to calculate it is to work on the logarithmic scale, which is supported by R's distribution functions. Likely one of these functions is used internally to calculate the p-value here.
– Roland
Nov 21 at 16:46












I'm jealous this is a problem I never have... :)
– Ian Wesley
Nov 21 at 17:10




I'm jealous this is a problem I never have... :)
– Ian Wesley
Nov 21 at 17:10

















active

oldest

votes











Your Answer






StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














 

draft saved


draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53414335%2fget-full-p-value-in-r-when-its-less-than-2-225074e-308-machinedouble-xmin%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown






























active

oldest

votes













active

oldest

votes









active

oldest

votes






active

oldest

votes
















 

draft saved


draft discarded



















































 


draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53414335%2fget-full-p-value-in-r-when-its-less-than-2-225074e-308-machinedouble-xmin%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

A CLEAN and SIMPLE way to add appendices to Table of Contents and bookmarks

Calculate evaluation metrics using cross_val_predict sklearn

Insert data from modal to MySQL (multiple modal on website)